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Chapter  1 


INTRODUCTION  AND  SUMMARY 


1-1.  Introduction 

This  Volume  II  of  this  final  technical  report  of  four  volumes 
documents  a  ray  tracing  radome  analysis  computer  program  written 
in  Fortran  IV  for  use  on  the  Cyber  70/74  computing  system  at  Georgia 
Institute  of  Technology  and  the  IBM  3033  computing  system  at  Johns 
Hopkins  University  Applied  Physics  Laboratory.  The  program  was  developed 
at  Georgia  Institute  of  Technology  over  the  past  four  years;  however, 
considerable  development  work  in  computer  aided  radome  analysis  has  taken 
place  here  prior  to  that  time  [1-7]. 

This  analysis  package  was  used  during  the  research  carried  out 
under  this  grant  to  analyze  the  antennas  and  radomes  using  the  fast  re¬ 
ceiving  formulation  as  described  in  Volume  I.  Its  documentation  was 
done  in  conjunction  with  the  on-going  radome  technology  program  at  JHU/ 
APL  under  the  cognizance  of  R.  C.  Mallalieu  (APL  Contract  601053) .  It 
is  intended  to  serve  as  part  of  a  technology  base  for  the  radome  tech¬ 
nical  community. 

The  report  is  organized  by  chapters  according  to  the  approximate 
order  in  which  the  subprograms  are  called,  and  each  chapter  describes 
one  subprogram.  Each  chapter  is  essentially  self-contained  since  it  is 
meant  to  serve  as  the  complete  documentation  on  a  single  subroutine . 
References  are  provided  at  the  end  of  each  chapter.  In  some  cases, 
figures  are  duplicated  in  different  chapters  for  completeness.  Each 
chapter  is  terminated  with  the  listing  of  the  subroutine. 
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Chapter  2  describes  the  main  program  and  instructions  for  its  use. 


Chapters  3  through  28  describe  the  thirty  four  subroutines  required  for 
execution,  including  those  for  producing  Calcomp  pattern  plots  and  three- 
dimensional  plots.  Appendices  A  through  D  present  computed  results  for 
four  test  cases  for  use  in  verifying  correct  operation  on  other  systems. 
These  results  were  obtained  on  the  Cyber  70/74  computing  system  at  Georgia 
Tech.  The  remaining  part  of  this  chapter  describes  background  of  the  pro¬ 
gram  development  and  summarizes  the  features  of  the  computer  analysis. 

This  report  comprises  Volume  XI  of  four  volumes.  Volume  I  describes 
the  salient  results  of  this  overall  investigation  to  determine  the  accu¬ 
racies  and  ranges  of  validity  of  various  analysis  methods.  Volume  III 
documents  the  additional  software  required  to  analyze  radomes  using  a 
surface  integration  method.  Volume  IV  presents  the  experimental  results 
obtained  and  is  meant  to  serve  as  a  data  base  for  other  investigators 
reeking  to  verify  the  accuracy  of  their  computer  codes. 

1-2.  Background 

Development  of  the  radome  analysis  computer  program  (RACP)  was 
initiated  in  1971  in  an  effort  to  include  the  effects  of  the  radome  on 
a  ground  mapping  radar  [1] .  A  three-dimension  geometry  and  vector  field 
formulation  were  used.  A  plane  wave  spectrum  (PWS)  representation  of 
the  radiation  from  the  antenna  greatly  facilitated  the  computations  since 
tin :  Fast  Fourier  Transform  ( FFT)  could  be  used.  The  program  was  used  to 
compute  power  patterns  on  the  ground  for  many  different  cases  of  antenna/ 
missile  orientations.  From  these  data,  the  effects  of  the  radome  on 
pattern  shape,  power  loss  and  VSWR  were  determined. 


Monopulse  tracking  antennas  were  next  introduced  into  the  computer 


analysis  to  evaluate  radome  materials  and  shapes  for  seeker  systems  in 
the  8-18  GHz  band  [2|.  Tangent  ogive  shapes  of  various  fineness  ratios 
were  analyzed.  Monolithic  and  multilayer  wall  structures  were  used. 
Algorithms  were  developed  to  compute  boresight  errors  from  the  sampled 
data  difference  patterns  in  two  orthogonal  planes.  A  modification  of  this 
program  was  also  used  to  conduct  a  trade-off  and  development  study  for 
the  Multipurpose  Missile  (MF'M) ,  later  known  as  ASALM  [3]. 

The  next  step  in  the  development  of  RACP  came  in  1977  witii  the 
introduction  of  a  conical  scan  tracking  antenna  into  the  analysis  [4]. 

This  antenna  necessitated  a  reformulation  of  the  analysis  from  the  trans¬ 
mitting  formulation  used  earlier  to  a  receiving  formulation.  The  big 
advantage  offered  by  the  latter  is  that  the  antenna  response  can  be 
calculated  for  only  one  direction  of  arrival  of  tile  target  return  (plane 
wave).  In  the  former,  the  FFT  automatically  computes  "responses"  for 
many  directions  of  arrival  and,  lienee,  is  computationally  slower. 
Subsequent  versions  of  the  program  have  used  the  same  receiving  for¬ 
mulation  with  monopul  so  and  other  types  of  antenna  models. 

The  computed  results  obtained  with  the  receiving  and  transmitting 
formulations  are  not  always  the  same  [5].  A  computed -aided  analysis 
which  utilizes  the  Huygens-Fresnel  principle  [6,  7]  is  generally  con¬ 
sidered  to  be  more  accurate  than  the  two  methods  already  mentioned,  but 
requires  considerably  more  computation  time  that  may  not  be  warranted 
in  all  cases.  A  research  program  is  now  underway  at  Georgia  Tech  whose 
objective  is  to  establish  the  accuracies  and  ranges  of  validity  of  these 
three  methods  of  radome  analysis  | ’>  1  . 
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I  -  ( . 


Descri  pt  ion  of  the  Analysis 


Tin-  current  version  of  the  ray  tracing  analysis  computer  proqram 
utilises  a  receiving  formulation  based  on  the  Lorentz  reciprocity  theorem 
1  >)•  A  plane  wavt'  <<f  selectable  linear  or  circular  pol  ari  ration  is 
assum'd  ihcidint  on  !  In'  outside  of  the  radome  and  is  represented  by  a 
svst  em  id  parallel  ravs.  There  is  one  iay  for  each  sample  data  point,  in 
tiie  antiima  aperture  inside  the  radome.  Kaeh  ray  is  (raced  from  the  point 
where  it  iim  inges  on  the  outside  surface  to  the  corresponding  aperture 
point.  The  electric  and  magnetic  fields  K  ,  II  associated  with  each  ray 

'l  ‘  1 

are  weighted  by  the  c!at  panel  transmission  coefficients  ,  T j j  as  d<- 

‘  .  nni:.-  d  1  v  the  unit  normal  n,  the  di  recti  on  of  propagation  k,  and  the 

dti  b  etrie  pioperlieK  of  the  radome  wall.  The  weighted  incident  f  ields 

K '  ,  11'  at  each  aperture  point  art*  then  used  in  the  fol  lowing  integral  to 

obtain  the  complex  voltage  response  V  of  the  antenna  as 

r 


V  -  C 
r 


[ 


( F,  X  II’  -  lb  X  11  ) 
T  a  -  l  -T 


dxdy 


(!' 


where  1-;  ,  n  are  the  aperture  fields  whi'ii  the  antenna  is  transmitting, 
i'  is  a  complex  constant,  and  i  s  the  unit  vector  normal  to  the  xy 
< . i j  ,  rl ure)  plane.  f  or  digital  computer  implementation,  the  integral  in 
i'.quat  ion  (1)  reduces  to  a  double  summation,  and  the  equal-area  elements 
dxdv  i»  come  /  xAv  and  can  be  absorbed  into  the  constant  P. 

In  it:;  present  form,  t  he  program  accomrnod.it  es  only  one  radome 
shale;  viz.,  the  tangent  ogive.  The  length,  diameter  and  fineness  ratio 
■  if.  ,  of  ■  •<  Hu  :  •  ,  all  variable  in  tin  input  data.  Monolithic  and  tnulti- 
i  .tv.  r  will  eeuf  i'liirat  i"ia.  an  be  analyzed;  however,  only  uniform  wall 
■  mini  it  ieiis  wiles'  (toi.  itie:;  do  not  vnv\  from  point  to  point  on  'he 


wall  can  be  handled.  Provisions  are  made  to  allow  for  a  metal  tip  on 
the  radome  whose  effect  is  aperture  blockage. 

The  geometry  subroutines  provide  for  three  separate  coordinate 
systems  and  the  point  and  vector  transformations  among  them.  A  reference 
coordinate  system  is  provided  to  orient  the  antenna/radome  combination 
with  respect  to  other  bodies.  The  coordinate  systems  for  the  antenna 
and  the  radome  comprise  the  other  two  systems.  Boresight  error  and 
pattern  computations  are  carried  out  and  expressed  in  the  antenna  coor¬ 
dinate  system. 

The  primary  outputs  of  the  program  art  boresight  error  (mrad.), 
boresight  error  slope  (deg. /deg.),  gain  loss,  and  when  selected,  prin¬ 
cipal  plane  patterns.  Outputs  include  both  printing  and  ['lotting 
(Calcomp) .  Plotting  options  allow  for  selection  of  aperture  fields 
with  and  without  the  radome.  A  feature  is  also  [provided  to  either  ob¬ 
tain  or  suppress  intermediate  calculated  results  for  debugging  purposes. 

Boresight  error  calculations  for  monopulse  antennas  are  carried 
out  by  setting  the  first  target  return  at  a  known  direction  within  a 
few  degrees  of  true  boresight .  The  rosi'onses  in  the  two  difference 
channels  and  the  sum  channel  art'  then  computed  and  stored.  Another  set 
of  responses  for  a  return  180°  away  from  the  first  is  computed  next. 

The  two  sets  of  data  are  then  used  to  construct  a  linear  tracking  model 
in  the  two  orthogonal  planes,  and  the  process  is  repeated  until  a  bore¬ 
sight  null  is  indicated.  The  true  direction  of  arrival  of  the  plane 
wave  at  this  point  represents  the  boresight  error  directly. 

The  current  subroutine  used  to  characterize  the  antenna  permits 
selection  of  various  polarizations  and  two  aperture  distributions.  A 
uniform,  circular  aperture  distribution  having  vertical,  horizontal  or 


ci  rcular  (UIR  or  RHC)  polarizations  is  one  combination.  Tin  second 
distribution  is  a  tapered  rectangular  distribution  having  vortical  polar- 
i. ’.at  ion  as  lound  in  flat  plate  antennas.  This,  basic  subiout  t  no  would  not 
be  difficult  to  modify  to  accommodate  other  distribut  ions,  such  as  t  <  c  - 
t angular  aperture  with  cosine  taper. 

Computation  time  is  independent  ol  radome  size  but  depends  on  the 
number  of  samples  used  in  the  aperture.  For  2 66  sample  points  (ib  X  16 
array) ,  the  time  to  compute  the  received  voltaqes  in  the  three  cl  auiels 
is  1 . h  seconds. 

Thi'  proqram  is  organized  as  a  main  program  and  a  number  of  sup- 
poiting  subroutines,  all  written  in  Fortran  IV.  The  complete  program, 
including  plotting  software,  contains  thirty  four  subroutines.  The  core 
storage  required  for  the  complete  program,  including  a  LI  library  and 
system  I/O  routines,  is  just  over  46,000  (decimal)  words.  Integer,  real 
and  complex  variables  and  arrays  are  utilized.  Single,  double  and  three- 
dimensional  data  arrays  are  present.  Only  single  precision  variables  and 
computations  are  required  with  the  l>0-bit  word  available  on  the  Cyber  70 
at  Georgia  Tech. 
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Chapter  2 


PROGRAM  FTFRACP 

2-1.  Purpose:  RTFRACP  is  a  Fortran  computer  program  used  to  analyze 

the  effects  iif  a  tangent  ogive  radome  on  the  performance  of  a  mono- 
pulse  aperture  antenna.  It  consists  of  a  main  program  ant  11  sub¬ 
routines.  It  user,  complex  arithmetic  and  requires  1712  1  octal 
words  of  core  memory  for  execut  ion  on  t  In  Cbe  Cyber  Vo  system  (( ...  lit 
words)  at  Georgia  Institute  of  Technology.  Execution  time  to  com¬ 
pute  boresight  error  on  the  Cyber  70  is  approximately  two  seconds 
per  look  direction  when  the  antenna  aperture  is  represented  by 
16  x  16  =  216  sample  data  points.  Execution  time  to  compute  trans¬ 
mitting  and  receiving  patterns  and  aperture  near  fields,  and  to 
compute  the  necessary  Cal comp  ommands  for  two-  and  three-dimensional 
plotting,  is  a[  proximate ly  35  seconds  for  one  look  direction. 

The  computer-aided  ri'oiv  analysis  us<.  s  a  i  •*«  -eiving  forrtuiat  ion 
based  on  the  Lorontz  reciprocity  theorem  as  described  earlier  [1,2).  The 
voltage  produced  at  the  terminals  of  a  linear  antenna  by  an  incident,  plane 
wave  is  given  by 

Vk)  =  III  x  Hr  "  £r  x  H^)  •  n  da  (1) 

where  E  ,  H  are  the-  fields  produced  on  the  surface  S  enclosing  the  antenna 
when  ttie  antenna  is  transmitting;  E  ,  H  are  the  incident  fields  produced 
on  S  by  the  incident  plane  wave  or  perturbations  thereof;  k  is  a  unit  vector 
which  points  from  the  antenna  toward  the  direction  from  which  the  i  lane 
wave  arrives;  arid  n  is  a  unit  vector  normal  to  the  surface  S  and  pointing 
tru 1  ,-kdb  C  ilO  .  Hi 


outward.  The  fields  E^,  are  taken  to  be  those  produced  in  the  planar 
ap.-rt.uri  when  the  antenna  is  transmitting  in  the-  absence  of  the  radome . 
The  geometrical  optics  approximation 


H 

-T 


x 


(2) 


i:.  used  to  generate  l  lie  magnetic  field  in  the  aperture  from  the  aperture 

illumination  specified  by  E  .  Rays  are  traced  from  each  sample  point  in 

T 

the  uprture  in  the  uirection  k  to  the  inner  radome  wall.  The  plane  wave 
to- Ids  associat'd  with  each  ray  are  weighted  with  the  flat  panel  insertion 
Voltage  transmission  coefficients  as  determined  by  the  radome  wall  configu¬ 
ration,  *  he  angle  of  incidence,  and  the  plane  of  incidence.  The  individual 
>nt  »  i  but  ions  an  summed  u(  as  indicated  in  Equation  (1). 

Tin  ;  a rami  t  .g  s  of  the  tangent  ogive  radome  are  indicated  in  Figure 

-1.  ’l'lii  out  id.  base  diameter  D  and  fineness  ratio  F  determine  thi 

os  os 

\it  .  i  i.  1.  ngt  !,  a.  cord  l  ng  to 


F  -I.  /D 

os  os  OS 


(  3) 


A  mi  I  u  r.  lot  ini,  hold-;  for  t  t ,,  inside  dimensions;  viz., 


F  =-  I,  /!■ 

1  S  1  S  1  s 


(4) 
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iirv.it  nr.-  of  the  outside  wall  R  is  given  by 

os 


R  E  I)  /  sin  (»  - 


Tan_1(2F  )) 

os 


(M 


Id 


and  the  dimension  B  is  given  by 


B  =  R  -  D  /2  (G) 

os  os 

The  placements  of  a  bulkhead  (bottom  disk)  and  metal  tip  (top  disk)  can 

be  specified  by  Z  _  and  Z  ,  respectively.  The'  thickness,  dielectric 

BOT  TOP 

constant,  and  loss  tangent  of  the  wall  may  also  be  specified  for  up  to 
N-S  layers.  The  radome  is  assumed  to  be  a  body  of  revolution  with  uni¬ 
form  wall  dimensions  independent  of  location.  The  dashed  cylindrical 
shape  of  a  diameter  D  in  Figure  2-1  was  used  earlier  to  simulate  a  laser - 
induced  defect  and  is  not  pertinent  hero. 

The  subroutine  which  generates  the  antenna  aperture  fields  ropre- 
sent:-;  two  tyt  os  of  antennas:  circular  aperture-  with  uniform  illumination 
and  any  one  of  four  polarizations  (vertical,  horizontal,  RHC,  LHC);  flat 
[  1  a  ft  antenna  with  tapered  illumination  and  vertical  polarization .  For 
"ith.'  r  antenna,  the  fields  are  computed  for  oil"  of  three  selected  channels 
.'.um,  azimuth  difference,  elevation  difference.  Inputs  include  the  number 

sum)  ha  rj  ,  N  and  the  apertun  diameter  D  ,/>.  in  wavelengths. 

X  Y  A1 

Thi  ant.  tma/ radome  orient  at  ion  is  q  > oified  according  to  the  para- 
•  i  (.  fined  in  i  itiur.,  The  atu)l<  ;  s>  I’.’ts  the  plane  of  scan  of 

t  : i •  to. Some  t  ii  with  ,.s|.  ct  t<>  the  antenna  coordinate  system:  ^  =  0° 

•  |.  at  t  t:  azimuth  ;  lain  ;  1  *  *<  s.  led;,  t  lie  eh-v.il  Ion  j  lane.  The 

anal.  *  i :  if.  til  in  tin  so  1 .  <  t •  d  (lane. 

in.  :  re  p  air  .'Mmi'it.  .  boresiil.t  enors  in  th.-  azimuth  and  el  vat  inn 

ill!..  •  t  I  •  .•  a,  t  ■  nil  I .  Tf ,.  f  a.  (|  urn  orientation  is  -,pe.  i  f  I .  .1  by  and  !  . 

P  !• 

7 1  .  tii  a  t  ,  i  ;•  i  r  •  t  •!»  n  i  i  Kino  w.iv.-l  is  made  to  arriv<  from  the  Ii  nit  mn 
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strns  Used  in  Radome  Analysis. 
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x  sin  0  +  y  sin 
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whore  8  is  the  initial  specified  offset  angle;  c.g.,  2°.  The  voltage 
os 

received  by  each  channel  is  computed  and  stored.  The  second  return  is 
made  to  arrive  from 


k  =  x  (-sin  0  )  +  y  (-sin  0  )  +  z  /l  -  2  sin^  0  (8) 

2  A  os  A  os  A  os 


and  the  voltages  are  again  computed.  The  data  from  these  two  points  are 
used  to  construct  a  linear  tracking  model  in  the  two  planes,  and  a  direc¬ 
tion  of  arrival  k  is  predicted  which  will  yield  null  indications  in  both 
planes.  The  process  is  repeated  until  a  desired  error  tolerance  is  satis¬ 
fied  or  a  maximum  number  of  iterations  is  exceeded.  Upon  completion,  the 
output  k  indicates  the  direction  from  which  the  plane  arrives  which  yields 
an  electrical  boresight  indication.  If  a  and  6  represent  the  boresight 
error  angles  in  the  azimuth  and  elevation  planes,  respectively,  then  they 

arc  related  to  the  direction  k  =  x  k  +  y  k  +  z  k  by 

Ax  Ay  A  z 

k 

x 

sin  a  =  >-  =  ==  (9) 

A  -  k 
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sin 


(10) 


z  x  y 


(11) 
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Options  dte  also  provided  whereby  principal  plane  patterns  as 
shown  in  Fiqure  2-3  and  additional  outputs  around  boresight  can  be  com¬ 
puted  and  printed.  These  options  are  useful  when  preparing  software  for 
a  new  type  of  antenna  and  to  ensure  correct  operation  whenever  curious 


results  are  obtained. 

2-2.  Usage:  Line  No. 

DATA  APIN/O .  /  47 

DATA  ZROTIN  >  > .  tH '/  49 

DATA  RADIUS/lEO/  52 

DATA  THETAA ,  PHI A,  AGAM3A/0 .0,90.0,0. 0/  53 

DATA  NX,  NY,  NXE ,  NYE  ,  NXY/16 , lb , 25b , 1 , 512/  56 

DATA  NREC,  NS,  MX,  MY/32,16,16,1/  57 

READ  (5,6)  TITLE  62 


READ  (5,*)  GRAF3D,  GRAFSA,  GRAFTR ,  GRAFRV,  SUTPRS,  IPENCD  65 
READ  (5,*)  NFINL,  NFHI ,  NTHE ,  DIAOS,  RA,  RR,  ZTOPIN,  FREQ. 
OSANG  67 

READ  (5,*)  LMAX,  DMRAD,  I OPT,  RAPMAX,  VAIRM,  I POL,  I  CASE, 


N,  IPWR 

76 

READ 

d,  *) 

DIN  ( I )  ,  ER ( I )  ,  TD ( I ) 

(1=1, N) 

108 

READ 

(6,*) 

FINR(I) 

( 1=1 , NFINE ) 

117 

READ 

(5,*) 

PHI  (I) 

(1=1,  NPHI ) 

120 

READ 

('->,*) 

THETA ( I ) 

(1=1 , NTHE) 

122 

2-3.  Arguments 

a.  Inputs.  Units  of  arguments  on  input  are  distances  in  inches, 
angles  in  degn-.-s,  and  fro  (tt-ivy  in  g  i  gahi  at:’. ,  un  bum  otherwise  noted. 
Units  of  arguments  passed  to  subroutines  are  centimeters,  radians,  and 
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£  =  elevation  component 

A 

a  =  azimuth  component 


gigahertz.  An  asterisk  is  used  to  denote  those  DATA  arguments  that  do 


not  normally  need  to  be  changed  by  the  user. 

APIN*  -  Height  of  a  cylindrical  base  section  of  the  tangent 

ogive  radome.  It  is  no  longer  included  in  the  ray 
tracing  algorithms  and  should  not  be  changed  from 
its  zero  value. 

ZBOTIN  -  Distance  from  base  of  tangent  ogive  radome  to 

missile  bulkhead  (Figure  2-1). 

_ -i  If  D 

RADIUS*  -  The  radius  R  used  in  the  far  field  factor  e  J  ,'R 
by  Subroutine  FAR.  Do  not  change. 

THETAA*  -  Angle  0  between  z-axis  and  the  position  vector  r 

a  —a 

to  the  antenna  origin.  This  angle  was  used  in 
earlier  work  to  locate  the  antenna  origin  in  the 
reference  system  using  spherical  coordinates 
(r  ,  0  ,  ).  Do  not  change .  S'  o  Chapter  7. 

c3  3  3 

PHIA*  -  Angle  if  between  the  projection  of  axis  onto 

the  xy-plane  and  the  x-axis.  Do  not  change. 

AGAM3A*  -  Angle  between  z^-axis  and  z-axis  in  Figure  2-2. 

Do  not  change. 

NX, NY  -  Integer  powers  of  two  equal  to  the  number  of  sample 

points  in  the  antenna  aperture;  e.g.,  lb,  32,  P-4 , 
ill'.  Changing  NX  and  NY  necess i tutos  compatible 
changer,  in  Dines  16-18. 

NXE.NYE  -  Integer  powers  of  two  which  specify  the  expanded 
number  of  sample  points  desired  when  computing 
the  transmitting  patterns  of  the  antenna  by 


i  nv  r 


Fourier  transforming  the  aperture  fields. 


Subroutine  JOYFFT  provides  this  capability  of 
increased  resolution  in  one  or  both  dimensions. 
Changes  in  NXE,  NYE  necessitate  compatible  changes 

in  Lines  16,  20,  22,  and  23.  Note  that  NXF,*NYESNX*NY 
and  either  NXE<NX  or  NYEXNY. 

NXY  -  Integer  power  of  two  used  by  Subroutine  JOYI-TT  for 

dimension  of  complex  working  array  XYFFT.  Note 
that  MX*NXFNXY  and  MY*NY'  NXY.  See  below  for  MX 
and  MY. 

NREC  -  Integer  power  of  two  equal  to  the  number  of  joints 

at  which  to  compute  the  receiving  pattern  in  either 

principal  plane.  The  received  voltage  is  computed 

at  points  t!  equally  spaced  in  sinO,  where  v  is 

the  angle  measured  from  the  z  -axis  as  indicated 

A 

in  Figure  2-3,  where  sin  0  =  -  K.MAX  +  ( I  —  1 )  *2*KMAX./ 

iiKEC,  and  where  KKAX  -  sin  0  <  1.0. 

max 

NF  -  Not  used.  It  was  originally  used  by  .Subroutine 

RE CBS .  Do  not  remove. 

MX, MY  -  Integer  powers  of  two  equal  to  the  magnification 

factors  desired  in  the  k  and  k  (E-j  lane)  direc- 

x  y 

tions,  respecti vely ,  of  the  transmitting  antenna 
I  attorns.  Note  that  the  restrictions  MX*NY  'NXY 
and  MY*NY'NXY  must  bo  observed.  The  data  cited 
above  indicated  increased  resolutions  in  the.  NX 
direction  of  MX- 16  and  no  magnification  ( MY  —  1 )  in 
t'h>  NY  direihion.  Consequently,  note  that 
N  XH  -  M  X  *  N  X  - '  g  •• . 


1H 


TITLE  -  A  Hollerith  string  of  up  to  72  characters  which 

describes  briefly  the  analysis  being  done.  A 
format  of  18A4  is  specified  and  should  work  for 
machines  with  word  length  greater  than  or  equal 
to  32  bits.  The  dimension  of  TITLE  (Line  31) 
should  be  at  least  18. 

GRAF 3D  -  A  logical  variable  used  to  control  the  plotting 

of  the  incident  fields  on  the  antenna  aperture. 

This  feature  has  been  removed  from  the  program, 
and  GRAF 3D  should  always  be  FALSE. 

GRAFSA  -  A  logical  variable  which  (if  TRUE)  controls  the 

plotting  of  the  transmitting  power  patterns  of 
the  antenna  as  follows:  E-plane  sum,  E-plane 
difference  equation  (A  ) ,  H-plane  sum,  and 

CjXj 

H-plane  difference  azimuth (A,  ).  The  radome 

AZ 

is  absent. 

GRAFTR  -  A  logical  variable  which  controls  the  plotting  of 

the  amplitude  and  phase  of  the  antenna  aperture 
fields  in  the  following  order: 

Exx'  eyj:'  exaei/  eyael,  exaax '  eyaaz‘ 

GRAFKV  -  A  logical  variable  which  controls  the  plotting  of 

the  receiving  patterns  of  the  antenna  with  radome 
in  the  same  order  as  specified  under  GRAFSA  above. 

SUPERS  -  A  logical  variable  which  controls  the  printing  of 

numerous  results  as  illustrated  in  the  test  data 
m  Section  2-6  below.  When  TRUE,  the  printing  of 
these  numerous  results  are  suppressed.  This  feature 


1 ') 


r punch 


NFTNh 


NPHI 


NrniK 


m  a 


HA 


HP 


is  tMiivi'ii  i  i'll  t  to  aid  in  debuqqinq  new  portion:; 
of  .software  prior  to  makinq  produotion  rums. 

An  inteqer  variable  which  so  loots  pen  and  paper 
for  Mm  Cal  romp.  This  variable  may  be  system 
dependent.  For  the  Cyber  70,  IPENCD=00  yields 
ballpoint  pen  and  11"  wide  plain  paper;  IPENOD-40 
yields  a  heavier  ink  pen  and  tile  same  paper. 
Tnteaer  variable  equal  to  the  number  of  fineness 
ratios  to  be  considered  lor  t  Ik  tanqent  ogiv< 
radome;  e.q.,  NFlNl’-l. 

Tnteqer  variable  equal  to  the  numlier  of  scan  plain 
e.q. ,  NPHI=2 . 

Inteqer  variable  equal  to  the  numlier  of  angles 
in  each  scan  i  lane  at  which  to  compute  boresiqht 
errors,  etc.  Note:  The  program  is  set  up  to 
iterate  on  fineness  ratio,  scan  plane,  and  scan 
anqle  as.  outer  loop,  middle  loop,  and  inner  loop, 
respect  ively.  Therefore,  for  each  of  NI-'INE  fine¬ 
ness  ratios,  the  analysis  will  be  done  for  NTHF 
sran  anqles  in  NPHI  different  scan  planes. 

Real  variable  equal  to  the  outside  base  diameter 
(in.)  of  the  radome.  See  Figure  2-1. 

Real  variable  equal  to  the  distance  (in.)  from 
the  qimbn]  point  to  the  antenna  aperture. 

Real  variable  equal  to  the  distance  (in.)  from 


the  qimbal  point  to  the  base  of  the  radome. 


ZTOPIN 


Real  variable  equal  to  the  distance  (in.)  from 


FREQ 


OSANG 


LMAX 


DMRAD 


I  OPT 


RAP  MAX 


VAIRM 


TPOL 


the  base  of  the  radome  to  the  face  of  a  metal 
ti]i  on  the  radome. 

Real  variable  equal  to  the  frequency  of  operation 
in  gigahertz. 

Real  variable  equal  to  the  offset  angle  in  degrees 
at  which  the  first  target  return  is  to  arrive  on 
the  antenna;  e.g.,  OSANG=3.0. 

Integer  variable  equal  to  the  maximum  number  of 
iterations  allowed  by  Subroutine  RECBS  in  com¬ 
puting  boresight  error;  e.g.,  LMAX=5. 

Real  variable  equal  to  the  tolerance  in  milliradians 
allowed  on  computing  boresight  error;  e.g.,  DMRAD=0.1. 
Integer  variable  which  selects  the  polarization 
of  the  incident  plane  wave  as  follows: 

1.  Linear,  elevation  component 

2.  Linear,  azimuth  component 

3.  Right  hand  circular 

4.  Left  hand  circular 

Real  variable  equal  to  the  maximum  radius  (in.) 
of  the  antenna  aperture.  See  Figure  3-1. 

Real  variable  equal  t  :>  the  maximum  amplitude  of 
sum  channel  received  voltage  without  radome.  Any 
real  value  can  be  entered  for  this  variable  since 
a  subsequent  program  modification  (Lines  326-328) 
causes  VAIRM  to  be  computed  automatically. 

Integer  variable  which  selects  the  polarization 
of  the  antenna  when  ICASE=1  according  to  the 
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same  code  as  used  above  for  IOPT. 


I  CASE 


N 


IPWR 


DIN , ER 


FINK 


PHI 


Integer  variable  which  selects  one  of  two  typos 
of  antenna  apertures  for  the  analysis:  ICASE-1 
or  2  selects  a  circular  aperture  with  uniform 
illumination;  ICASE=3  selects  a  flat  plate 
antenna  with  programmed  illumination.  See 
Subroutine  HACNf'  in  Chapter  3. 

Integer  variable  equal  to  the  number  of  layers 
(up  to  5)  in  the  radome  wall.  For  cases  where 
more  than  5  layers  are  required,  the  dimensional 
arrays  on  Line  37  must  be  changed  to  NN=N+1. 
Integer  variable  which  selects  the  component 
for  which  to  compute  the  transmitting  power 
patterns  as  follows: 

1.  Elevation  Component 

2.  Azimuth  Component 

3.  Total  {lower 

TD  -  Subscripted  real  variables  equal  to  the  thickness 
(in.),  dielectric  constant  (t  ) ,  and  loss  tangent 
(tan  6)  of  each  layer  of  the  radome  wall.  1=1 
corresponds  to  the  first  layer  and  is  the  layer 
on  tile  exit  side  of  tile  wall.  Layer  N  is  the 
first  layer  encountered  by  the  incident  plane 
wave.  See  Subroutine  WALL. 

Subscripted  real  variable  equal  to  NFINE  fine¬ 
ness  ratios. 

Subscripted  real  variable  equal  to  NPHI  angles 


(degrees)  which  specify  the  scan  planes. 


THETA  -  Subscripted  real  variable  equal  to  NTHM  anqles 

(degrees)  which  specify  the  scan  angles  in  the 
scan  plane. 

b.  Outputs.  The  parameters  of  analysis  which  are  computed  and 
and  outputted  by  the  program  depend  on  whether  SUPPRS  is  true.  In  what 
follows,  it  is  assumed  that  SUPPRS=FALSE  so  that  all  possible  outputs 
are  obtained.  Since  many  of  the  original  input  parameters  are  printed 
directly,  only  those  parameters  not  already  explained  above  will  be  in¬ 
cluded  below.  Additional  clarification  may  be  found  in  Section  2-b. 

TABLE  -  Logical  variable  which,  if  TRIE ,  causes  a  look-up 

table  to  be  used  in  computing  transmission  coeffi 
cients.  When  SUPPRS=FALSE,  an  abbreviated  table 
of  transmission  coefficients  of  the  radomc  wall 
is  printed  by  Subroutine  WALL  with  variables  as 
explained  immediately  below. 

AN  (ILL  -  a  1  variable  equal  to  the  angle  of  incidence 

(degrees)  of  the  plane  wave  on  a  plane  sheet  of 
infinite  extent  having  the  layered  configuration 
specified  for  the  radome  wall.  The  entries  in 
the  table  are  computed  at  2r>0  equal  increments 
in  sin  0.,  but  only  every  fifth  result  is  printed 

TPERI , TPARI -  Complex  variables  equal  to  the  voltage  insertion 

transmission  coefficients  of  the  sheet  for  the 

two  cases  of  E  ( >erpendi cular  to  the  plane  of 

incident e  (T  )  and  E.  parallel  to  the  plane  of 

incidence  (T|.)  .  In  the  printed  table,  the 

power  transmission  coefficients  | T  ‘  are 
1 -  I  ±  l 
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1 
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KXMAX  -  Rea  l  variable  equal  to  the  folding  wavenumber 

associated  with  satn|  l  i  n  < )  the  ai"  ttuiv  f  i  *  *  I  ds 
accord inq  t  o  KXMAX  l  .  A  2(  Ax/  )  .  where  Ax  i  s  t  lu¬ 
ll  i  st  ance  hot  we*  n  samples.  See  Subr<  <ut  ines  HAl'NI' 
and  l-'FTA. 


i > XW I .  -  Ri  al  variable  e  ,uaJ  to  A.\/.\. 

KXM,KYM  -  Real  variables  equal  to  the  foKILtiq  wavenumbers 

of  the  principal  plane  pattern:;  after  maqnifi- 
eation  for  increased  resolution.  KXM  KYMAX*f»'XE/ 
(MX*NX)  and  applies  to  the  11-plane. 

KYM- KYMAX*NYK/ (MY*NY)  and  applies  to  the  E-plane. 
Usually,  the  expanded  dimension  NXE  and  magnifiea- 
t ion  factor  MX  are  selected  so  that  KXM- KXMAX. 
Also,  NYE  and  MY  are  usually  si  looted  so  that 
KYM<  <KYMAX. 

MIN, MAX  Real  variables  equal  to  the  minimum  and  maximum 

values  of  the  amplitude  of  the  complex  arrays 
con t a i n i nq  the  aperture  fields  as  processed  by 
Subroutine  NORM! I  in  preparation  for  3D  plotting 
by  Subroutine  P1.T3P11. 
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ROS 


Rival  variable  equal  to  the  radius  of  curva- 


BOS 

FINOS 

FINIS 

The  following  variables 
putod  and  printed: 

ICUT 

I  COMP 


KMAX 

NREO 


turc-  of  the  outside  shape  of  the  tangent  ogive 
radome . 

Real  variable  equal  to  the  distance  B  in  inches 
defined  ir.  Figure  3- 1 . 

Real  variable  equal  to  the  fineness  ratio  of 
the  radome  as  based  on  the  outside  dimensions. 

Real  variable  equal  to  the  fineness  ratio  of 
the  radome  as  based  on  the  inside  dimensions, 
are  printed  when  the  receiving  patterns  are  corn- 

integer  variable  which  defines  the  E-plane  (ICUT=1) 
or  H-plane  (ICUT=2)  pattern.  See  Figure  2-3. 

Integer  variable  which  defines  the  field  comp'onent 
of  the  plane  wave  incident  on  the  receiving  antenna: 
IC0MP=1  for  elevation  component;  IC0MP=2  for  azimuth 
component. 

Real  variable  equal  to  the  sine  of  the  maximum 
angle  off  broadside  for  which  the  received  voltage 
is  computed. 

Integer  variable  (power  of  2)  equal  to  the  number 
of  points  at  which  the  receiving  [.'attorn  is  com¬ 
puted.  Tlu  pattern  is  computed  at  NREC  points 

spaced  equally  in  k  =sinO  according  to  Ak  = 
xy  xy 

2  KMAX/NREC. 


DK 


Real  variable  equal  to  2 *KMAX/NREC . 
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b.  Supporting  Subroutines.  Thirty  four  supporting  subroutines 


are  required  by 
below . 

(1) 

(2) 

(3) 

(-1) 

CO 

(•  ) 


RTFRACP.  The  purpose  of  each  one  is  briefly  described 

HACNF — Computes  complex  vector  aperture  electric  fields 
of  antenna  for  all  three  monopulse  channels  at  NX  x  NY 
sample  points. 

ORIENT- -Com;  utes  matrices  ROTATE  and  TRANSLate  used  for 
coordinate  transformations  by  Subroutines  POINT  and 
VECTOR . 

POINT — Transforms  a  point  P(xft,  y^ ,  z  )  in  antenna  system 

to  the  same  point  I'(x  ,  y  ,  z  )  in  radomc  coordinate  system, 

K  H  H 

and  vice  versa. 

VECTOR — Transforms  a  vector  from  radome  to  antenna  coor¬ 
dinate  system,  and  vice  versa. 

1NC'PW--Oom|  utes  the  rect.  angular  electric  field  corn:  orient  s 

of  a  r  lan«  w.«v<*  incident  from  the  direct  ion  k  in  int.nna 
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dielectric  wall  with  unit  inner  normal  n.  The  unit  vectors 
k,  n  are  used  to  resolve  the  incident  plane  wave  into  vector 
components  perpendicular  and  parallel  to  the  plane  of  inci¬ 
dence,  and  to  determine  the  angle  of  incidence.  KXMI T  colls 
Subroutines  WALL  and  AMPHS. 

('')  WALL — Computes  the  voltage  insertion  transmission  coeffi¬ 
cients  of  flat  panel  model  of  the  radome  wall  as  function 
of  the  sine  of  the  incidence  angle. 

(10)  AXR — Comjut.es  real  vector  cross  product  C  =  A  X  B. 

(11)  CAXB--Computes  the  complex  vector  cross  product  C  =  A  X  B. 

(10)  RECBS — Computes  boro sight  errors  of  antenna  enclosed  by 

tin  radonu  for  the  specified  orientation,  fineness  ratio, 
etc.  RECBS  calls  Subroutines  INCPW,  RECM,  and  AMPHS. 

(13)  RECPTN- -Computes  receiving  patterns  of  all  three  channels. 

RECPTN  calls  Subroutines  INCPW  and  RECM. 

(1-1)  OGI VI :- -Computes  point  of  intersection  of  ray  and  ogive  by 
solving  a  quart:  ic  < -(ju.it  ion.  OGIVE  calls  Subroutines  CBRT, 

SCK,  .tin!  XY . 

(!'■)  'BET- -i.’omt  'lit  es  cube-  root. 

(  ltd  SV'K — i‘om|  ul .  -  -.guari  root  with  test  for  negative  argument. 

I  1 '/' )  OGIVE?.-  -v’onj  ut«  tin  unit  inward  normal  vector  to  the  ogive 

urt.i  -  .  t  til-  pi  - 1 1. 1  P  ( x  ,  v  ,  :  )  . 

K  'R  K 

1  i  XY--"  .  i  b-,-  Siil-iout  1  in  OGIVE  to  compute  the  x  and  y  com- 

R  *  R 

•--I:  :.  t  !  t :  o|.‘|t  of  i  titersect  ion  of  a  ray  on  the 
. !  •  ■  i  i  'i  n-  a  i  t .  e  -i  . 

i  •  •  I .  K  —  .  : i,t  -  ; .  :  joint  .  -t  i  tit .  rsect  ion  of  a  ray  and 

1  e  .  •  i .  t ’  r  i  i •  k  t  -  |  ■ .  : .  •  nt  1 1  •  <  i  t  hi  bulkhead  inside  the 


U'l) 


cvi 


C 2  3) 


(.’4) 


(25) 


(2C) 


(2?) 


(-’.») 


(2')) 


BOX SKN--Comf  u t: cs  unit  normal  vector  to  bulkhead  (n  =  rx  ). 

K 

TDItiK — Comj >ute  s  the  point  of  intersection  of  a  ray  and 
tin'  base  if  tin-  metal  tip  on  the  radome. 

TDISXri — Com  |  ut.i  unit  normal  victor  to  metal  tit  (n  -  ~z  ) 

!< 

F’AK — Conj  -ut  es  tlie  ami  lit  u«.i»  of  tin  powr  oattein  from  t  In 

Corns  lex  :  lain  wave  set  ,1  ni  A  Ik  ,k  ),  A  (k  ,k  )  of  at. 

x  x  v  y  x  y 

ante  nna . 

AMI  lb' — Converts  a  comt-li  x  ■  umber  from  r>  ctanyul  ,ir  to 
'polar  form.  This  subteul  ine  utilizes  the  intrinsic 
function  AT  All  2 .  Tin  ami  I  i  t  ud.  •  ; reduced  is.  linear  (not 
decii.ils),  and  ttw  j'hasi  is  in  degree;;  on  the  ranue 
(-1.SO,  180). 

IIBPV — Converts  a  rial,  t wo-dim.  nsh.nal  array  from  linear 
to  logarithmic  values  in  dicibeis  on  tin  ramir  e  fi.-4')db. 
NORMH — Normalises  a  two-di inonsionni  real  arr.iy  t o  values 
between  P  and  1. 

CNPLTH--rlots  single  dimensional  far  field  (attorns  on 
axes  patterned  aft  or  standard  pattern  recorder  paper. 

CNPLTli  calls  Subtout  i  lie  PSi  in  addition  to  tins  usual 
raloom;  subrout  iiws. 

PS  I — lined  by  Subrout. ine.  CNPLT1I  to  compute  the  asiii.ut'nal 
angle  y. 

PLTHDtI~~Yic  Ids  three-dimensional  \  lots  of  tne  data  in  tin 
two-d imonsi onal  real  array  FIELD.  PLT.-.DH  cal  ls  Subroutines 
PLTT,  NORMH  as  Well  us  t  tie  usual  Co  1  Com)  subroutines. 

1‘LTT — Used  by  Subroutine  PLTtnil  to  eliminate  irovinc  the 
pen  for  hidden  lines. 


* 
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ist  Foui'ii  r  'i'r-.itisforin  of  t 


Linos  40-41’:  Declare  namolists  for  print inq  data.  These 

namelists  are  no  longer  used  except  for  occasional 
debuqq  inq  I  'urj  >oser . 

Lines  43-57:  Set  data  in  DATA  statements  as  described  above 

in  Section  2-3. 

Lines  61-62:  Set  SMAX  and  VMAX  to  unity  to  prevent  division 

by  core. 

Lines  (>3-64:  Read  ami  write  TITLE  aceordinq  to  18A1  format. 

Lines  b5-67:  Read  input  data  using  free-field  format. 

Line  68 :  Compute  sine  of  the  offset  anqle  0 

^  OS 

Line  69:  Set  TAR1J5=  FALSE  so  that  normalizing  factor  VAIRM 

ea  be  computed  (bines  t|q-.S2<>)  via  .)  call  to 
Subrou Lines  RECM  and  RXMIT.  In  the  latter,  TABLE” 
FALSE  causes  T  ,  T,,  to  be  set  to  unity  as  in  the 

ill  r 

case  of  no  r a dome . 

Lines  71-75;  Write  input  data. 

Li  nes  76-77:  Read  in]  lit  data  and  set  VAIRM  needlessly. 

Lines  78-104:  Comments  explaining  input  variables. 

Line  105:  Set  NN=N+1-  Number  of  wail  layers  plus  one. 

bine  1 06 :  Initialize  DINClb  total  thickness  of  radome  wall 

in  inches . 

biii'-s  107-100:  Read  wall  da  ca  anil  compute  total  thickness, 

tone  Mi);  Compute  DIAIN  inside  base  diameter  of  the  radome 

in  inches. 

bines  111-112:  Compute  indices  of  the  center  element  of  near-field 

ariavs  coni  spending  to  x  -y  -0. 

'  A  'A 


Liiii-s  113-114:  Write  array  dimensional  data. 

Lines  11  '-133:  Read  fineness  ratios,  scan  plain's.,  and  sran  an. ties. 

Lines  l.'3-13t>:  Compute  wavelength  in  inches  and  cent  inteters. 

Coinput  e  3  =  3  n /A 

cm 

Lines  L37-13S:  Call  RXMTT  and  compute  table  id  transnnssie.  coeffi¬ 
cient;;  versus  sine  of  incidence  angle.  The  first  call 
to  RXMIT  builds  the  table.  Subsequent  calls,  use  the 
table  if  TABLE=TRUE. 

l,i IK*  129:  Compute  DAPWL=  diameter  of  antenna  aperture  in 

wavelengths. 

Lines  130-139:  Convert  variables  in  inches  to  centimeters  for 

input  to  subroutines.  Some  variables  are  multiply 
defined  to  avoid  conflicts  in  labeled  common;  e . g .  , 

ZBOT  and  Zl.  Note  that  DIACM  is  the  inside  diameter 
of  the  radome  in  centimeters. 

Lines  1,40-144:  Convert  angles  from  degrees  to  radians  using  RAD=n/lB0. 

Lines  14'j— 1:>1:  Compute  near  fields  of  three  channel  monopulse  antenna 
using  Subroutine  IIACNF. 

Lines  1L2-1LH:  Let  KYMAX'KXMAX,  compute  magnified  folding  wavenumbers 
KXM,  KYM,  and  print  results. 

Lines  J  1  1 '  —  1  77:  Initial  i/.i  Ca  1  comp  plotter,  if  required.  The  commented 
i  n  i  .  t  i  a  1  i  za  l  i  <  >n  (him:;  I  til-  1  74 )  applies  to  the  IBM  3033 
system  at  J1IU/APL. 

Note:  hi  in's.  r/H-.'Stt  are  used  to  plot  the  near  fields  of  tile 

antenna  and/or  the  transmitting  principal  plane  |  ow.  r 

(.It  f  <  I  IIS  . 
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Lines  178-179:  Initialize  the  maximum  values  FMXEL,  FMXDAZ  of 

the  E-  and  H-plane  patterns  so  tliat  when  used  ini- 
tially  as  inputs  to  Subroutine  FAR,  the  resulting 
pattern  will  be  normalized  with  respect  to  its  own 
maximum  and  FMXEL  and  FMXDAZ  will  be  set  equal  to 
these  respective  maxima.  On  subsequent  calls  to  FAR, 
the  resulting  patterns  will  be  normalized  with  respect 
to  FMXEL  and  FMXDAZ.  Hence,  the  relative  gain  of  the 
difference  and  sum  patterns  will  be  correctly  dis¬ 
played  in  the  graphs. 

Line'  180:  Iterate  for  each  of  three  monopulse  antenna  channels. 

Lines  181-190:  Equate  complex  arrays  EXT,  EYT  to  the  selected  near 
field  and  compute  the  amplitude  NF  of  EXT. 

Line  191:  Assume  transmitting  near  fields  are  to  be  plotted 

( GRAFT R=T)  . 

Line  193:  Call  Subroutine  PLT3DH  to  plot  the  amplitude  of  EXT. 

The  inputs  XSIZE=6.,  YSIZE=2.5,  HEIGHT=2.5  yield  a 
3D  plot  that  will  fit  on  a  B'i"  x  11”  report  page. 

The  inputs  NF,  NX,  NY  specify  the  real  array  to  be 
plotted  and  its  dimensions.  The  input  NMZ=.TRUE. 
directs  the  subroutine  to  normalize  NF  so  that  its 
values  be  between  0  and  1.  The  input  LDR= . FALSE 
indicates  that  the  array  NF  contains  linear  values 
rather  than  logarithmic  values  (decibels) . 

Lines  194-201:  Compute  and  plot  phase  of  EXT  on  a  scale  of  -180 

degrees  to  +180  degrees.  Note  that  Line  19')  ensures 
that  the  real  array  NF  contains  these  phase  values 
scaled  to  the  required  0  to  1  range. 
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Line 


202-21  ' :  Repeat  amplitude  and  phase  3D  pl  at  s  far  EYT. 


Line  die: 


Line  2Jl>: 


L.ltie  2  2 : 


Line  .22  3: 


Line  224: 


Li  ip 


la  in 


Assume  GFAFSA-T  so  Liiat  principal  f  lane  patterns 
are  plotted. 

If  I P- 3 ,  go  to  Line  243  and  plot  H-plarie  patterns; 
otherwise,  [  lot  E-plane  patterns. 

Call  Subroutine  JOYFFT  to  calculate  the  inverse 
Fourier  transform  of  the  x  -component  of  near  field 
EXT  to  produce  the  plane  wave  spectrum  XEEL  from 
which  the  radiation  field  can  be  computed.  In  the 
process  of  computing  the  transform,  provide  increased 
resolution  from  NX  x  NY  points  to  NYE  x  NXE  points 
through  the  point  (NXC,NYC)  in  the  array  EXT.  In 
tile  k^  direction,  the  ['lane  wave'  spectrum  is  magnified 
bv  MY;  it.  is  magnified  by  MX  in  the  k  direction. 

y 

The  array  FFTXY  is  a  working  array. 

Repeat  for  EYT  to  produce  the  plane  wave  spectrum 
YEEL  for  the  y  -component  of  field. 

Call  Subroutine  FAR  to  calculate  the  E-plane  elevation 
( I  PWR=  3)  power  pattern  FFSKI,  of  the  near  field  at 
equal  samples  in  sinC  over  the  range  (-KXM,  KXM  -AK) . 
[f  FMXEL'1  (and  it  is  for  ir=l),  normalize  FFSEL  with 
respect  to  its  own  maximum. 

Call  Subroutine  PBPV  and  convert  the  power  pattern 
tii  dei -ib'  Is  on  .i  scale  of  0  t.o  -40  dB. 

Scale  the  values  in  FFSEL  to  the  range  of  0  to  1  for 
I  1 1  >1  t  i  1 1  • i . 


3b 


IrrtT  i  i  ill  >J^ii 


Line  231:  Call  Subroutine  CNPLTH  and  plot  the  power  pattern. 

If  KXM<1,  the  pattern  is  plotted  over  the  angular 

range  corresponding  to  sin  ^ (KXM) ;  if  KXM>1,  the 

angular  angle  is  (-90°,  90°).  Subroutine  CNPLTH 

actually  plots  conical  cuts  corresponding  to 

k  =  constant  or  k  =  constant  as  specified  by  inputs 
x  y 

KXC,  KYC.  In  the  call  here,  KXC=KYC=0  so  that  a 
principal  pattern  is  produced. 

Lines  232-236:  Write  a  figure  title  for  the  plot  and  establish  a 
new  origin  for  the  next  plot. 

Line  237:  If  IP=2,  the  E-plane  patterns  are  finished. 

Lines  238-242:  Since  JOYFFT  changes  the  input  arrays  EXT,EYT,  it 
is  necessary  to  recompute  them  so  that  increased 
resolution  can  be  obtained  in  the  plane  wave  spectra 
in  the  H-plane. 

Lines  243-258:  Repeat  computation  and  plotting  for  H-plane  power 
patterns . 

Line  260:  Iterate  the  radome  analysis  for  NFINE  fineness  ratios 

Line  261:  Set  FINE  =  outside  fineness  ratio. 

Lines  262-266:  Calculate  and  write  R  ,  B,  F  ,  F  as  defined  in 

OS  O SIS 

figure  2-1  for  the  radome  geometry. 

Line  267 :  Compute  RDML  =  distance  from  the  base  of  the  radome 

to  the  theoretical  tip  on  the  inside  of  the  radome. 
Lines  268-272:  If  ZTOPIN<RDML,  the  radome  has  a  metal  tip,  and  a 
message  is  written  to  that  effect. 

Lines  273-283:  Compute  parameters  needed  by  Subroutine  OGIVE  to 

describe  the  radome  shape.  R  and  B  are  in  centimeter 
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and  apply  to  the  inside  dimensions.  AP,  the  height 
of  the  cylinder  in  centimeters,  is  not  used.  RTSQ= 
square  of  the  radius  of  Lhc  top  disk.  RBSQ=  square 
of  the  radius  of  tile  bottom  disk  (bulkhead).  The 
other  variables,  BSQ,  RINV,  RSO I  ,  RP,  and  RP2 ,  are 
precalculated  hero  to  speed  later  computations  in 
OGIVE. 

Line  Compute  conversion  factor  DPMR  for  converting  milli- 

radians  to  degrees. 

Lines  286-288:  Initialize  the  "last"  values  of  bore sight  error  in 
azimuth  ( A2L)  and  elevation  (ELL)  and  the  "last" 
value  THL  of  scan  angle.  These  variables  are  used 
later  to  compute  boresight  error  slope  in  degrees  per 
degree  from  the  ;  resent  and  last  values  of  boresight 
error. 


Lines  289-290 : 
Lines  291-203: 
Lines  294-296: 
Lines  297-301: 

Line  309: 

Lines  310-312: 
Line  313: 

Lines  314-316 : 

Line  317: 


Write  title  for  analysis  results. 

Write  parameters  of  radome  wall. 

Write  heading  for  table  of  boresight  error  and  gain  data. 
Write  this  same  data  to  logical  unit  7  for  subsequent 
storage  as  a  disk  file,  if  desired. 

Iterate  the  radome  analysis  for  NPHI  scan  planes. 

Compute  4>r  in  radians  as  required  by  Subroutine  ORIENT. 
Iterate  the  analysis  for  NTHE  scan  angles  in  each 
scan  plane. 

Compute  t!  in  radians  as  required  by  Subroutine  ORIENT. 
Cali  Subroutine  ORIENT  and  compute  the  rotation  matrix 
ROTATE  and  translation  matrix  TRANSL  requi red  for  coor¬ 
dinate  transformations  using  Subroutines  TOINT  and  VECTOR. 
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Line  318:  On  the  first  iteration,  TABLE  is  false  so  that 

the  maximum  amplitude  of  the  received  voltage  on 
the  sum  channel  is  computed  without  the  radome. 

Line  319-322:  Set  the  direction  cosines  of  the  incident  plane 

wave  so  that  it  arrives  from  the  z  direction. 

A 

Line  323:  Call  Subroutine  INCPW  and  compute  the  rectangular 

components  PWI  of  the  incident  plane  w.  v  aving 
polarization  specified  by  IOPT. 

Lines  324-325:  Set  TSUP=T  and  TABLE=F  so  that  an  air  radome  wall 
be  used  and  so  that  printing  by  Subroutines  RXMIT 
and  RECM  will  be  suppressed. 

Lines  326-327:  Call  Subroutine  RECM  and  compute  the  complex  vol¬ 
tages  VR  received  on  the  sum,  difference  elevation, 
and  difference  azimuth  channels,  respectively, 
corresponding  to  VR(I) ,  1=1,3. 

Line  328:  Compute  VAIRM= |vR(l) | . 

Line  °?9:  Set  TABLE=T  so  that  on  subsequent  iterations 

VAIRM  will  not  be  recomputed,  and  so  that  the  table 
of  transmission  coefficients  will  be  utilized  when 
RXMIT  is  called. 

Lino  33i>:  If  SUPPRS-F,  compute  and  print  the  E-plane  and  H- 

l'lane  receiving  power  patterns  of  the  antenna  with 
the  radome  in  place. 

Lines  333-334:  Iterate  in  J  for  E-plane  (ICUT=1)  and  H-plane 
(ICUT=2)  patterns. 


Line  335: 


Set  the  desired  far  field  component. 


I .  i  III 'i 


Line 

Lin.': 

Li  lit': 

Lint.'.' 

Line 

Line 

Line; 
I.i  ne 
Line; 

Li  tie 


,  336-337:  Set  KMAX=sin  1  (0  )=.(il>6.  If  KXMAX,  as  comt  ut.  d 

max 

by  UACNF,  is  less  than  KMAX,  then  use  tile  sm.il  It  r 
as  the  maximum  angle  in  the  principal  plane  at 
which  to  compute  the  pattern. 

338:  Set  the  temporary  logical  variable  TSUP=T  so  that 

printing  will  be  suppressed. 

:  33.'- 14.':  Call  Subroutine  RECPTN  and  compute  the  com)  lex 

i  ■  e.  i ved  voltages  on  each  of  three  channels  at  NRLC 
points  over  the  range  (-KMAX,  KMAX  -  DK )  . 

:  54  1  -  544  :  Increase  the  resolution  and  print  results  for  all 
three  channels.  Do  not  print  results  that  are 
known  to  be  identically  zero. 

343-346:  Transfer  the  received  voltage  into  a  one-dimensional 
array  VREC. 

147:  If  NREC'-'NXE,  there  is  no  need  to  increase  the 

re so  1 ut ion . 

348:  Call  Subroutine  MAGI’ FT  to  increase  t  h<_  resolution  <  f 

VREC  from  NRLC  points  to  NXE  points.  The  re.su  1 1  is 
contained  in  complex  array  XYII'T  on  output  . 

3  3  4 '  t  —  3  ‘ '  3 :  Compute  linear  power  pattern. 

134:  Select  NXX-larger  of  NXE  and  NRLC. 

■?  133-336:  Write  heading  for  printed  result  •:  from  Subroutiiu 
NORM!  I . 

t  .7:  t.\  1  Subroutine  NORMII  to  normalize  the  NXX  values  in 

real  array  MVKEC  to  bil  between  zero  and  one.  The 
input  argument  LDB  .FALSE,  since  the  values  . 

in  dee i be  I s . 
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are  not 


Li.ru;  3bH: 


Call  Subroutine  UBl'V  ti 


'■invert  t  tie  |mwer  (  ti  t  1**1!: 


in  MVRKC  te  dec  l  be  1  s  . 

Lines;  3V-3f<0:  Write  eorrei  t  heading  f<>»  K-|  lap.  «*i  11- ;  . 

Lj  tie  301 :  Compute  the  increment  in  sintl  at  wimii  the  |  ■.  iw<  •  r 

pattern  I. is  been  computed  and  res<  Ived. 

Lines  30 30rt:  Scale  the  power  pattern  to  have  v  dues  between 

and  1.  If  Sl’PI,kS  =  F,  r-nmj  Uto  t  h*  all'll.  ••  Afjo  and 
the  ph  ise  ot  the  pattern,  and  [ rint  the  results  for 
i  very  fourth  anqle. 

Line  37J:  If  0RAFRV=T,  [lot  the  receivinq  powei  [attorns. 

Lines  373- 37H:  Call  Subroutine  CNPLTH  and  plot  the  receiving 
patterns  in  turn.  Write  an  appropriate  figure 
titb  following  each  pattern  plot.  Re-origin  the 
plotter  pen  for  subsequent  plots.  The  result  of 
Lines  130-383  is  four  principal  plane  patterns: 

t’.-plane  sum,  E-plano  A  ,  !'.  plane  sum,  H-plane  A  . 

E  L  AZ 

Lines  .18-1-  1H(. :  fall  Subroutine  PR  CBS  and  compute  the  boresight  errors 
AZT,  KLT  in  the  azimuth  and  elevation  planes  of  the 
antenna  as  caused  by  the  radome.  On  output,  the  real 
array  KA  contains  the  direction  cosines  of  the  last 
target  return  and,  hence,  gives  the  true  direction 
to  the  target  at  the  time  that  the  tracking  functions 
m  the  azimuth  and  elevation  planes  indicated  the 
electrical  boresight  direction. 

Line  3H7 :  If  this  is  the  first  iteration  in  scan  angle,  do  not 

attempt  to  compute  boresight  error  slope. 
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l'i  Jin?  Ut  .  I  ..  ’! 


I.i  i  i  ■  : 


L  1 ;  1 1  n  1  1  —  > 


Line  ?l>  > : 


■sight  •rr.ii  slope  (  di  qr  . -i s/di -qrei  )  in 


i.'H’iuti.  .in..]  i  1.  vat  inn  .  imniii  Is. 

■  j  >  t  tn.  "last"  values  .it  lh.ii  i i  i|.‘it_  errors  and  s.  mu 
angle  t  •  ■  thi'  i  ui  1 1  i  '  values  it.  ;  1 1  ■ ;  a  t  a  t  1 1  .r,  t  r 
n. ■  X t  it.  i  at  ion  . 

I’umi  ut.  Kiss  i  ii  maximum  Main  id  tin  ant.-nna  mm 
i  ti. min  I  due  t_o  tin-  laiU'in.  . 


Li  rii’S  34 1  : 


l.l  tti’S  J'i.i-4.  ,,  , ; 


Ll  la  40i: 


Writ.-  result  s  to  logical  units  «•  and  7 . 

Writ.'  maximum  amplitude  <>f  received  .  on  voltage 
VAX  KM  witl  mt  radome . 

Terminate  plotting  software. 


STOP 


END 


.'-t*.  Test  Cases 

Four  test  cases  are  presented  in  Appendices  A,  B,  C,  and  D  to 
demonstrate  correct  operation  of  the  radome  analysis  computer  program 


i-TTRACP. 


Appendices  A  and  B  |  resent  the  test  data  and  results  for  a 
circularly  (RHC)  polarized  antenna  and  five-layer  tangent  ogive  radome 
at  a  frequency  of  11.8028l>GHz  (A-=1.0  inch).  The  diameter  of  the  aperture 
is  ll.b-li .  The  outside  diameter  of  the  radome  is  16.267  inches.  The  fine¬ 
ness  ratio  is  3.00.  In  Appendix  A,  the  program  is  exercised  without 
plotting,  and  printing  is  minimized.  In  Appendix  B,  all  plotting  and 
printing  options  are  exercised. 

Appendices  C  and  D  present  the  test  data  and  results  for  a  vertically 
polarized  flat  plate  antenna  of  diameter  ri .  11 02 A .  All  other  parameters  of  the 
analysis  are  the  same  as  in  Appendices  A  and  B.  Appendix  C  contains  the 
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Chapter  3 


SUBROUTINE  HACNF 

Purpose:  To  compute  near-field  aperture  distributions  for  two 

types  of  three-channel  monopulse  antennas:  ( ] )  circular  aper¬ 
ture  with  uniform  amplitude  and  phase  distributions;;  (.  )  flat 
plate  antenna  with  a  programmed  amplitude  distribution  and  uni¬ 
form  phase.  Four  polarizations  can  be  selected  for  the  circular 

aperture.  The  flat  plate  antenna  is  vertically  (y  )  polarized 
only. 

Usage:  CALL  HACNF  (E,  NX,  NY,  ICHAN,  IPOL,  IXY,  DAPWL , 

DXWL,  KXMAX,  XCASE) 

Arguments 

E  -  Complex  array  of  NX  by  NY  elements  which,  on 

output,  contains  the  values  of  the  specified  (IXY) 
rectangular  component  (x^  or  y  )  of  the  electric 
field  distribution  over  the  specified  (ICASE) 
antenna  aperture  having  the  specified  (IPOL)  pola¬ 
rization  for  the  specified  (ICHAN)  channel  of  a 
throe-channel  monopulse  antenna. 

NX, NY  -  Ev.-n  integer  number  of  points  in  a  rectangular 

array  at  which  the  aperture  distribution  is  com¬ 
puted  in  the  x^  and  y  directions,  respect ivoly . 


The  ;  oint  I-NX/C  +  1,  J=NY/2  +  1  corresponds  to 


:  V.'  -  :  :  •  •  'I.  t  r. . !  Variable  wi  til  Values  J  ,  ,  <  'I  4 

A  :  i .  i  •  ,  •  t  )  i -  .  lii: i ,  ■  •  1 1  v.i  t  it ' : i  d  i  1  t  a  r » -  in  »  ,  -  l 

i  mu’  :  ;  i  r  !  i  r  ■  ii-  •  hanui  I  ,  t  :  i  •  •  ;  v<  1  y . 

;  i  .  ::.t  t.-l  v.u  i  >t  1  wiiii.  >  U  .1:;  *  in.  antenna 

1 . 1 1  i  . ".  i  *  II!  i :  t  ■  i  1  ■  'w  : 

:  V<  i  t  i  - 1  !  J  v  „  '  1  1 1  1 1  -I. 

A 

-  It  -r  i  :;oi,t  a  I  (xr  ) 

A 

■  -  in  ■  i  n  u  1  u  i  " 

i  -  1,.  f  t  -hand  '  i  t  '.il.tr 

1 XV  -  int.  g.  i  o  1 1 *  >•  1  v.in.i;  b  having  values  1  or  .1  t< 

:  i  *  t I s,  ■  i  y  eonpom  nt  t..f  ai  •  riurr  c  U-ctrio 
fu  iii. 

:  ini  -  :  >  i  -  tint  t'-i,  1 1  -  wav-  1.  ngths,  of  tin.-  antenna  ;i[-i  rturt  . 

:  XW!.  -  .t.-aiiiuj,  in  wave  I  •  ngt  I  is ,  1"  twe.  u  samples  in  aper- 

lui'i  in  x  and  v  directions  (output.)  . 

A  ’  A 

KMAX  -  Maxi  mini  Va  1  u-  of  normalised  wavonumlior  cor responding 

to  K MAX  -  1 .  /  (  d  . *DXWL)  (output). 

ICASi:  -  Inteij.-r  control  variable  having  values  1  or  2  t;o 

•  ci  fy  a  circular  aperture  antenna  with  uniform 
am;  1  i  t  ude  and  |hase.  If  T CASH -3,  a  flat  plate 
antenna  having  a  programmed  amplitude  distribution 
(see  Table  3— 11 )  with  Vertical  polarization  is  *•  leet 
:itn.  lit.  md  M.  t  h-  d 

Ii't  i  NX,  NY  must  -  arh  bi  equal  to  eaeh  ollu  r  and  to  an 
'..■r  .  tv;  .  .  \x  NY-  It..  In  addition,  when  I  CASK-'  (1.1  at  1  late 


b.  The  act  ual  shape  of  the  circular  aperture,  a:;  upj  roxi mated  by 
a  rectangular  array  of  sample  points,  is  shown  in  Figure  3-1  for  the  case 
of  NX=  JY  —  lh .  Row  l  and  Column  1  of  the  array  contain  null  element:..  Thi 
elements  inside  and  on  the  boundary  of  the  aperture  may  contain  non-zero 
values  as  shown  in  Table  3-1  for  the  various  cases  when  ICHAN=1  (sum 
channel) .  Note  that  specification  of  D  in  Figure  3-1  determines  the 
sample  s;  ucings  according  to 


=  Ay 


D  cosa  D  cosa 
Ap _ _  Ap _ 

(N  -2)  "  (N  -2) 

x  y 


(1) 


wnere  a  -  Tan  ^ (2/7) . 

The  aperture  distributions  for  three  monopulse  channels  are  formed 
by  phasing  the  elements  in  the  four  quadrants  of  the  aperture  appropriately . 
The  sum  channel  distribution  is  formed  by  assigning  equal  phases  to  all 
elements.  The  azimuth  difference  channel  is  formed  by  multiplying  all 
elements  in  Quadrants  II  and  III  of  the  sum  distribution  by  minus  one  and 
by  zeroing  all  elements  along  x  =0 .  For  the  elevation  difference  channel, 
Quadrants  III  and  IV  are  negated,  and  all  elements  along  the  line  y  =0 
i  r •  made  Zen  for  symmet  ry  reasons. 

]  hasing  chosen  models  a  tracking  antenna  and  provides  outputs 
in  two  orthogonal  channels  from  which  the  direction  of  arrival  of  a  target 
r>  turn  can  be  mathematically  determined.  Let  k  be  a  unit  vector  which 
joints  from  the  antenna  origin  toward  the  dire'etion  from  whence  the  plane 


wav.>  (target  return)  emanates;  i.e.. 


A 


OF  C 
R  GR 


FIGURE  3-1.  APPROXIMATION 
BY  RECTANGULA 


k  =  x  k  +  y  k  +  z  k 
Ax  Ay  A  ; 


I H  ■  f  i  Ik*  the  tracking  functions  for  this;  |  lane  wave  a:; 


A. (k  , k  ) 

r.(k  rk  )  -fyr--Xr'{ 

1  x  v  T,  ( k  ,  k  ) 
x  y 


(2) 


(  3) 


where  \  represents  the  output  of  tile  elevation  (i)  or  azimuth  (a)  difference 
channel  and  represents  the  sum  channel  output.  Then  for  small  k  >0,  the 
phase  of  f  is  +  ti/2;  for  small  k  <0,  the  phase  of  f  is  -n/2.  Similarly, 

1  <  X  t)! 

for  small  k _/0,  arg  (f  )  =  n/2;  for  small  k  <0,  ary  ( f  )  =•  -n/2.  Hence, 
the  change  in  phase  by  n  in  either  channel  represents  the  bores i gilt  direc¬ 
tion  of  the  antenna,  and  tracking  is  doin'  using  the  imaginary  parts  of 
the  trucking  functions  rather  than  their  real  parts. 

c.  Tile  shape  and  sump]  ing  grid  used  to  model  the  flat  plate  antenna 
are  shown  in  l-’igur.  3-2.  In  Subroutine  H.’.CNF,  the  integers  NX  and  NY  must 
both  equal  le,  and  only  linear  polarization  (y  )  is  applicable  to  the  flat 
j  late  antenna  (ICASK-3) .  The  phasing  of  the  four  quadrants  is  done  as 
ics-  r i hi  d  above  to  model  the  three  monopulse  channels  so  that  tracking  can 
b' ■  simulated.  Note  that  specificat  ion  of  l>^  determines  the  sample  spacing 
m  c<  u  ■  i  i  iiy  to 
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FIGURE  3-2.  GEOMETRY  OF  FLAT  PLATE  ANTENNA 
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Symmetrical  Amplitude  Distribution  for  Flat  Plate  Antenna 


Sample  No. 

X 

A 

Amplitude 

^A 

Am]  'lit  ud> 

9 

0 

1.0280 

0 

1.0280 

10 

Ax 

1.0280 

Ay 

1.0280 

11 

2  Ax 

.9120 

2Ay 

.917; 

12 

3  Ax 

.  7959 

3Ay 

.8060 

13 

4  Ax 

.6077 

4Ay 

.  o  1 53 

14 

3  Ax 

.4194 

5Ay 

.  4250 

1:> 

6  Ax 

.2097 

6Ay 

.2125 

16 

7  Ax 

0.0 

7Ay 

0.0 
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field  value  zero  (Line  40) .  Multiply  the  non-zero 


elements  by  CFAC(IPOL)  to  generate  the  correct  pola¬ 
rization  (Line  38) . 

Lines  42-43:  Compute  sample  spacing  Ax  /A  and  go  to  statement  60. 

Lines  44-46:  Error  message  and  STOP. 

Line  47-48:  Flat  plate  antenna —  if  NX^1<>,  write  error  message 

and  STOP  (Lines  109-111). 

Line  49:  Compute  sample  spacing  Ax  /A. 

Line  30:  Ensure  NX=NY 

Lines  51-84:  Zero  all  elements  in  the  aperture.  If  I XY=1 
(x^-component ) ,  go  to  statement  60. 

Lines  58-62:  Assign  tapered  amplitude  values  to  eight  "even" 
elements  in  Quadrant  III. 

Lines  63-71:  Compute  amplitude  values  for  the  "odd"  elements 
in  Quadrant  III. 

Lines  72-75:  Compute  amplitude  values  for  elements  3-9  along 

y  =0  line  and  along  x  =0  line. 

A  A 

Lines  76-79:  Generate  symmetrical  amplitude  values  in  Quadrant  IV. 

Lines  80-83:  Generate  symmetrical  amplitude  values  in  Quadrants 
I  and  II. 

Line  84:  Compute  k 

xmax 

Lines  85-89:  Test  to  determine  if  the  sum  channel  data  generated 

should  be  phased  to  produce  the  aperture  distribution 
for  a  specified  difference  channel  (ICHAN) . 

Lines  '<0-98 :  Form  aperture  distribution  for  difference  elevation 
channel  by  zeroing  all  elements  along  y^=0  and 
negating  all  elements  for  y^<0.  RETURN. 


Lines  99-107:  Form  aperture  distribution  for  difference  azimuth 


channel  by  zeroing  all  elements  along  x  =0  and 

A 

negating  all  elements  for  XA<C|-  RETURN  . 

Lines  108-112 :  Error  message  for  ICASE=3  and  NX^lft.  Comment  of 

DAPWL=5.047  applies  to  the  test  described  in  Chapter 

END 

Test  Case:  See  discussion  in  Chapter  2. 

References 

1.  D.  R.  Rhodes,  Introduction  to  Monopulse,  McGraw  Hill,  New 
York,  1959. 

Program  Listing:  See  following  pages. 


3U3°CJT  I  Nr  ^ArNF«£,NK,NY,ICHAri,IPCL«IXV,OAPHL,OXWL,KXHAX,ICASE» 
SUBS  HACNS  COMPUTES  ELECTRIC  FIELUCOHFONENTS  CVFR  a  CIRCULAR  APERTURE 
JF  SACILC  F.MAX=  (N*/ ?-l» /COS  <  ATAM  2/ 7|  l  AN  3  ®E  TURNS  SAME  IN  ECNX.NY). 
NX  "UST  EOUAL  NV  ANC  MUST  iF  c-VFN. 

IC  m  A  K=  1  FOR  SUM  CHANNEL  IPCL=1  E  OR  tfERT-Y  PC  L  .  IXY=1  FOR  X-COMP. 


X  (\J  rO  ^UNvCKflC0'OHN^^^vDN®0'O^fVjW4^^N(C 


a 

X 

o 

o 

I 

► 

Or 

o 


«\j 

ii 


o 

a 


z> 

a 


a 

z  x< 
<x  xi 


LU 
UJ  l/) 
</>  <X 
-a  o 

o  ^ 


o  • 

tjj  o 
O  Uj 

a 

o 

(T  o 
CT 

V 

CD  > 
Hi 

> 


X 

_> 

-J 

z 

J 

H 

Ul 

X 

1 

o 

o 

X 

u 

o 

l 

uj 

M 

a 

a 

w 

or 

a 

cr 

cr 

to 

O 

X 

z> 

• 

o 

to 

X 

O' 

o 

o 

tO 

X 

-J 

CLI 

o 

o 

<3 

o 

o 

i 

X 

X 

i- 

CL 

<X 

X 

X 

Ld 

d 

a 

T 

a 

-J 

X 

Id 

(\J 

-J 

X 

UJ 

a 

CL 

X 

— * 

II 

O 

o 

or 

O 

ft 

z 

V— 

<r 

• 

V 

o 

o 

X 

X 

CO 

a 

o 

o 

Id 

or 

H 

X 

Uj 

■d 

o 

u. 

u- 

u 

-J 

o 

UJ 

4- 

-* 

UJ 

w 

<* 

:> 

• 

• 

— » 

• 

• 

rvj 

x 

of 

■> 

*■» 

-1 

• 

VjD 

A, 

X 

fM 

• 

• 

a 

ii 

li 

II 

Ul 

H 

=> 

*■ 

rj 

V 

-*• 

in 

• 

• 

UJ 

UJ 

IVJ 

X 

O' 

D 

o 

<3 

W 

• 

H 

f\J 

o 

CT 

• 

• 

• 

1) 

a 

y: 

Z 

• 

■S’ 

O 

• 

u 

UJ 

X 

X 

CT 

z 

u 

X4 

z 

w 

o 

o 

• 

• 

U-1 

X 

CL 

Z> 

<-> 

Cu 

• 

z 

• 

o 

X 

cn 

X 

V 

• 

UJ 

o 

X 

.t 

O 

< 

4 

• 

CM 

X 

X 

> 

u 

a 

z 

ir» 

X 

o 

a 

X 

X 

»4 

■X 

X 

rv 

o 

X 

«r 

o 

• 

O 

»~4 

*j 

X 

o 

o 

1 

X 

— j 

CL 

Q 

<3 

H 

II 

z 

II 

II 

P-H 

X 

— 

o 

• 

• 

C* 

♦ 

d 

Z 

U 

O 

UJ 

u 

w 

C3 

d 

-J 

-J 

o 

w 

A 

X 

4 

X 

a 

cr* 

• 

u 

or 

X 

CD 

u 

x 

u> 

z 

LJ 

a 

u 

• 

c  -» 

o 

1 

1 

z 

/ 

II 

U 

Li. 

»u 

X 

X 

<x 

» 

a 

<y 

LO 

U 

a 

CL 

• 

o 

X 

"0 

— 

o 

<3 

X 

II 

o 

X • 

U- 

a 

o 

Z> 

z 

_J 

— 

• 

CJ 

►H 

►— « 

o 

O 

UJ 

♦ 

4 

CM 

o 

• 

« 

UJ 

X 

• 

Q 

u_ 

<x 

*r 

2 

3 

a 

> 

CJ 

o 

o 

• 

2L 

A 

4*. 

4 

U 

A 

X 

X 

UJ 

i 

« r 

• 

ho 

V 

• 

fO 

(VJ 

CM 

4 

X 

CM 

• 

l»J 

"> 

o 

U 

o 

> 

CL 

X 

♦ 

• 

• 

A 

•-K 

•H 

• 

V 

X 

> 

X 

• 

• 

• 

• 

U 

O 

<J 

O 

<1 

> 

H 

\ 

C? 

• 

• 

>- 

• 

Ou 

CJ 

X 

X 

> 

> 

4 

d 

O 

UJ 

u 

— 

•-* 

X 

-J 

u 

X 

u 

-J 

z 

• 

UJ 

» 

h~ 

t— 

z 

X 

» 

UJ 

z 

z 

z 

z 

CM 

X 

UJ 

UJ 

«S 

<  J 

— J 

•  i  1 

<x 

Y 

f  o 

U 

w 

X 

\ 

<\j 

• 

OJ 

o 

-j 

• 

X 

• 

*■ 

w 

*> 

4 

Of 

• 

• 

X 

► 

• 

• 

UJ 

r 

f\J 

U 

<r 

O 

Ld 

X 

• 

• 

id 

• 

Z 

UJ 

▼4 

1 

▼4 

1 

4 

• 

-J 

-J 

X 

X 

X 

a 

or 

X 

•d 

3 

Or 

X 

<1 

z 

to 

X 

-J 

-j 

z 

> 

w 

to 

II 

w 

II 

X 

X 

o 

a 

M 

— 

o 

a 

O 

UJ 

X 

a 

o 

X 

X 

U 

<j 

a 

o 

a 

• 

X 

Q 

<* 

X 

X 

X 

o 

a 

a 

> 

II 

II 

a. 

u 

u 

X 

a 

x 

o 

u 

UJ 

V 

UJ 

v- 

o 

u 

a 

V 

H 

o  o 

«x 

< 

X 

• 

X 

X 

U 

V 

X 

d 

X 

-J 

H 

II 

►-« 

z 

X 

X 

u 

o 

o 

o 

ur 

cr 

w 

w 

'd 

f\i 

X 

x 

< 

H 

m 

a 

-J 

<3 

»» 

•*— 

X 

'M- 

»*- 

H 

—I 

X 

-J 

C? 

** 

•• 

♦ 

h 

II 

II 

C 

-o 

X 

II 

n 

X 

LT 

■a 

u 

u 

to 

u 

o 

X 

II 

ii 

M 

Uj 

o 

Uj 

< 

z 

u 

X 

u 

u 

u 

u 

U 

u 

o 

II 

o 

II 

ii 

u 

u. 

u 

T 

J 

*•' 

U 

-J 

-i 

X 

to 

o 

or 

o 

<T 

(Y 

dM 

(-4 

h-4 

XI 

XI 

X 

o 

X 

o 

V 

or 

X 

X 

X 

cr 

-J 

Ui 

X 

1  I  J 
CL  X  X 
X 
V 


< 

o 


a 


u  u. 


)uoou  a  o  '  j  i '  •  5  « )  o 


bH 


F  <1  ,  J>  =F  (I,  J»  *U  F  AC  ( I  POL  I 


GO  TC  1C 
0  E  1 1  *  J>  =  C 
10  CONTINUE 


< J) 


N_c0<7'O-»HC\J*"),J,lf\>_0hwC0CT'o*H(\Jfk0.^tPv£)N.00<7NO 

h^N.f^C30Q0<3C<0a0®a000a0c00^(7»CT>(T‘<7'(T'<7'<T'CT'Cf'O 


-HCVJfO^-lf»vX)P^aO<r 

ooooooooo 


2T 

Or  2  Z 

3  rv  rv 
♦-  r>  Z) 

IU  ♦—  H- 

ft  UJ  UJ 

a-'  o' 


oj  ^ 

»  ii_i 
ii 

n  ~  u 
")  "5  3 
♦  Z’ 
» rN  o  *-* 

IfN  r  K* 
hH  I 

O  *- 
Q  U  O 


J  J  lit 
O  ON. 

a  a 


<\i  • 

•  •  C? 

U  O  UJ 
<iJ  I J  • 

•  •  21 

>-  X  >>1 
X  X  I 
*“«  K-»  O 


1  H  Q  (\J  H  H 
H  II  ^  N  II  II 

♦  tn  H  V  “5 

(M  -*  Z 


UJ  I  I.)  CVJ  ▼H 

|  ►"  ^  II  w  \  H  II 

II  D  ♦  T  II  X  h  1 

-  ^  r  (\j  —  2 

->  ii  *r\  irv 


v  ut  ")  n  fjc  la;  h  \  ii  -)  ii  irvirv 

>  vC  r  X  r^-  N*  ►  D  N  x  T  *  cc  <C 

ZL  >-♦  >“  «4  I  ►-<  *  t 

II  O  w  T.  o  o  ui  nowroo 

"5  O  lu  1  O  O  liJ  ft  n  M  O  lu  ^  p  o 

<x 


96  F0“Mir  <//"***cqcos  EXIT!  NX  NOT  “EQUAL  TO  16  IN  SU8R  H A CnF* ♦***// )  110 

STOP  111 

■E  NO  112 


*+ 


p 


I 


I 


Chapter  4 
SUBROUTINE  ORIENT 


4-1  . 


4-4  . 


rurjoste:  To  com) -uto  the  rotational  matrix  of  direct  ion  co:  i  j..  : 

ROTATE  and  tin-  translational  matrix  TRANSL  required  t<  carry  <  ait 

coordinate  and  vector  transformations  between  antenna  coord! i.ntn 

svstem  (x  ,  v  ,  z  )  and  radome  coordinate  system  (x  ,  v  ,  . 

A  /■  A  R  "  R  P 

1’sayo :  CALL  0RT1-NT  (RA,  THETA,  PIIIA,  RR,  THETAR,  PHIR,  AGAK3A, 

ROTATE,  TRANSL) 

Arquir.t,  r.t  s 


RA, 

Spherical 

coordinates  (cm 

,  radians) 

of  t  ill 

THETA 

oriqin  of 

the  antenna  coo 

ruinate  s’. 

•stem  with 

I'll  I A 

respect  ti 

-  the  reference 

ooord  i  r el  t  e 

■  system 

(x,  y,:-.)  a; 

:  indicated  in  1- 

i  quri  4-  1  . 

Not i  t  licit 

the  c  't:  i  it  i : 

1  of  tin-  refi  reuse  system 

coincides  wi 

Lin  si  pii  .  i  ; 

joint,  wlii  eh  i 

1  era ted 

on  t  lie  axi  s 

:  yriiin  ■  t  ry  . 

<  f  tiie  radomi 

• 

id-:, 

-  s!  her  1  I  .1  1 

coordinates  (oni 

,  radians) 

of  tiie  or  is 

Till  TAP, 

'f  tin  ran 

:<  Tii  coordinate 

svsUr:  wit 

il  T-I  'SI  i  Ot  t  f 

lap  i  he  ref.  r*  no*-  . cyst  > tti . 

AoA.’lA  -  AnciL  ( radians)  bctw«  cn  the  :•  and  s  ex-,  s. 

P'  TATI.  -  Real  attay  of  <  x  '  •  1  ■  rv  I .  t  wi  ;■  !  i  i .  r 

<>ut  )  ut  tin  ma  t  r  i  x  ■  !  « :  s  t  >  ■  •  t  i  <  ■:  -os  i  : .•  1 1- 

•  xi  I  a  i  in  4  I  "  low. 


TK, 


i-'i  ,  i  !  i  i  >  .  i  of  t  ;  a  ■  •  •  I  •  Tn  t,t  ; .  wi .  i  i  ■ :  i  t  ■:  1 1  a  l  ; , , 

■  ■  i i *  :  .if  *  : i ■  l  rails  I  a  t  i ,  a,  pint  r  i  x  T  as  •  x:  1  a  i  m  ■  i 

i 

:  .  1  -  w. 

.■  i 


4-4  . 


Comments  and  Method 


a.  The  coordinate  systems  of  Figure  2-2  are  obtained  from  those 


shown  in  Figure  4-1  by  setting  THETAA  =  0,  PHIA  =  ir/2,  and  AGAM3A  =  0. 


When  this  is  done,  the  z  and  z  axes  coincide,  the  x  and  x  axes  are 

A  A 


parallel,  and  the  y  and  y^  axes  are  parallel.  The  angles  <£_  and  0^  in 


Figure  2-2  are  related  to  0  and  d>  as 

r  r 


r 


■J>  +  TT 

P 


(1) 


6  =  TT 

r 


0 

L 


(2) 


The  unit  vectors  x  ,  y  ,  z  were1  chosen  to  coincide  with  the 
R  R  R 

spherical  coordinate  unit  vectors  r  ,  0  ,  <i>  associated  with  the  point 

r  r  r 

0  ( r  ,  0  ,  4>  )  ;  hence,  x  always  lies  in  the  olane  of  scan  as  indicated 

R  r  r  r  R 

in  Figure  4-2.  This  observation  is  important  if  the  properties  of  the 
radome  wall  are  not  symmetric  with  respect  to  rotation  about  tile  z  -axis, 
such  as  in  the  case  of  circumferential  variations  in  wall  thickness, 
nonuniform  heating,  etc. 

b.  Tiie  details  of  the  coordinate  system  transformations  are  de¬ 
scribed  in  Reference  1  and  summarized  below.  The  transformation  of  the 

point  P  in  antenna  coordinates  (x  ,  y  ,  z  )  to  radome  coordinates  (x  ,  y  , 

AAA  R  R 

z  )  is  given  by 
R 


(3) 


73 


T!u  t 1 ansforina t ion  of  a  vector 


F  -  x  F  +  y  F  +  z  F  =  x  F  +  y  F  +  z  F 

A  xA  vA  A  zA  R  xR  R  yK  R  zR 


(4) 


is  given  by 


xR 


vR 


zK. 


1  ! 


xA 


yA 


F 

L  zAJ 


(  :>) 


In  tlu:  above,  |R  )  is  the  matrix  of  direction  cosines  which 
i  .1 

ic scribes  the  rotation  of  the  radome  coordinate  system  with  respect.  to  the 
antenna  svstem;  IT  ]  describes  t'ne  Tanslation  of  the  radonv-'  oriqin  0  wit 

l  l 

r.'siwt  t.o  t.lu.  oriqin  O  of  the  antenna  svstem.  In  fact,  setting  (x  i  , 

a  ‘  A 

v  =  i.i,  z  -  •»)  in  liquation  (3)  shows  that  (T  ,  T  ,  T  )  represents  the 
A  A  x  y 

local  i ,  i.i ,  in  radome  coor.ii  nates,  of  the  aiitenna  origin. 

The  matrix  !  t  I  cat.  i  <  expanded  and  written  explicitly  as 


where  IR. . ]  denotes  the  transpose  of  [R.  ];  i.e.,  |R. .]  =  [P . . )  since 

i.l  i]  11  1i 

rows  and  columns  are  interchanged.  Also,  since  [R.  }  is  a  unitary  matrix, 

its  inverse  is  equal  to  its  transpose. 

To  facilitate  the  specification  of  a  particular  antonna/radome 
orientation,  the  reference  coordinate  system  (x,  y,  z)  is  used.  Trans¬ 
formations  from  the  reference  system  to  the  antenna  system  are  described  by 


x  -  r  sinO  cosdi 
a  a  a 


y  ,  v  -  r  sinO  sin f 

A  i  |  '  a  a  a 


z  -  r  cost* 
a  .i 


while  t.ransfonnat  ions  from  re  feronce  system  to  radome  system  are  described 


x  -  r  sinO  cosdi 
r  r  r 


r  sinO  si  in)' 
r  r  r 


/  -  r  cost! 
r  r 


wiieti  I  r  I  and  |,  I  rei-resi  nt  tin-  rotations  of  t  lie  two  systems,  with 
i  l  i  ! 


SI  •  e  t  (  e  the-  re  I  i  -r. 


;y;;tem.  When  these  two  separate  trails  forma  t  i  <  :s 


are  combined,  there  results 


(12) 


(13) 


where  X  ,  Y  ,  etc.,  .ire  d-fino  i  in  Equations  (  1  1  and  (11);  •  .g .  , 

&  cl 

X  =  r  sint'  cosi  . 
a  a  a  a 

4-r'.  Program  Flow:  See  listing  bi  low  and  compare  directly  to  method 
described  above.  Note  that  in  GA.‘-:(!,J)  ,  the  index  I  re;  re.  ■  nts  th.c  row 

number  in  [y  1,  and  index  a  •  •  resents  the  column  numi <  r . 

1 .1 

4-1’.  Test  Case:  Sec  discussion  in  Cii.g  ter  2. 

4-7.  Reference 

L.  E.  P.  Joy  and  G.  K.  Huidleston,  "iiadome  Effects  on  lie 
Performance  of  Ground  Mapping  V  idnr , "  Technical  Report, 
Contract  DAAHO 1 -J.’-C-On  :  •> ,  V.  Army  Missile  Command, 

March  1**73  . 

4-ft .  Program  Listing:  G<>e  fol  lowir.  :  j  i  :•  . 
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G  j  M  P' J  T  f  THf  PGTATt  APFAY  b  Y  MULTIPLYING  THE  PHG  APPAY 
A  N  C  Tmt  TPANSPOSf  OF  IMP  GA~I  APPAY. 


Chapter  S 


SUBROUTINE  POINT 

Purpose:  To  transform  a  roint  P  in  antenna  coordinates  {x  ,  y  ,  z 

AAA 

to  radome  coord  i  na  t  os  (x  ,  y  ,  z  )  ,  and  vice  versa. 

R  'R  R 

isayo:  CA.bh  'POINT  (i  ,  PT,  ATOR,  T,  PO> 

Armmeu  t: 

■’  Ke.ii.  i  .input)  array  of  t  itr-*  *-•>  o  .1  •.  nts  is  •  n  e-  a:.,  i  no 

the  Cartesian  coordinate:;  (on)  of  tie  :  oir.l  ' 
transformed;  e .  g . ,  I’(l)  =  xt  ,  i’(-)  "  V,.  I'd) 

NT  -  Heal  (output)  array  of  tlus  o  .  J  ru  nt  p  r-  i  i  •.*  ir.u 

t  hi  Car'esian  coordinator.  (cm)  ■  t.  tin  to  i  « .•  t- 

lir.  ite  s\  ■  t -  m ;  •  . a.,  PT(1)  --  x  ,  to. 

K 

ATOK  -  Log;  li  input  variable  wise;,  conti  •>;  ns.  i,  f 

t  i  ansforma  I  ion  :  ATOK  =.TRbb.  tier  ant*  <  -r,id<iin 

( ••  isnaticn  (■!-',));  AT'  T  .  EALSE .  f<>r  radon*  -*■  > 
..isc  ••niia  idii.a!..  trails  format  ion  jfi.-iuat  ion  (■;  —  ■-■)  ) 

Ri  . .  i  v  i :  i  ;•  1 1 '  )  array  of.  ?  x  '  i  i-  noids  i>  ;  n  e  in.: 

tin  ROTATE  array  •••oim  ut.«.d  nv  Bui  rout,  i  :i-  ORIENT . 

Ri  si  (ini nt)  array  of  three  elements  re;  resenting 
t  !:-■  TRANS  I.  array  computed  by  Subroutine  ORIENT. 

Conme.’: '  s  and  .’A  ti.ed 

i .  dm  rout  in.  s  i  ■  •  :u i i  ■  d :  Subroutine  ORIENT  must  be  cailed  prior 

first  .all  to  : oINT  bo  that  T  and  PO  are  available . 

i-  bi  r  nn  ‘hoi,  se  Suin',  ill;  i  lie  ORIENT  in  Chapter  ■!  . 

I  !  uia;  '  •  :  .'i:  . .  i  •.  !  i  :  I  i  tin  i  ■■  1  ov  t  i  r-  id  1  v  to  P  sis  I  '  or  ; .  i  i  -  <  V 

m  !  t  1-0  • 

,  ,  i i. iji'i.c 


r 

i 

f 


f'-6.  Tost  Case:  Sot;  Chapter  C. 

Reference:  Sot:  Ciiaptor  4. 

r>-s.  Program  Listing:  Sot:  following  pa  go  • 
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Chapter  (> 


SUBROUTINE  VECTOR 

6-1.  Purpose:  To  transform  a  vector  F  in  antenna  coordinates  to  radorru 
-  .xirdina tes ,  and  vice  versa. 

6-i.  Usage:  CALL  VECTOR  (V,  VT,  ATOR,  T } 

6-1.  Arguments 

V  -  Real  (input)  array  of  three  elements  representing 

the  rectangular  components  of  the  vector  to  be 

transformed;  e.g.,  V(l)  =  F  ,  V(2)  =  F  ,  V ( 3)  = 

xA  yA 

F  . 

zA 

77  -  Real  (output)  array  of  three  elements  representing 

t : : '  rectangular  conp  >:u  r.ts  of  tine  vector  i  n  t 

t  n."  s  a  -lain-p.  a.'stem;  ■  .g.,  VT  (!  )  -  F  .  ■  to. 

x  R 

.VI--  -  I, i u!  r  i.ai-.i  •  wh  i  eh  ■•••.w  tro]  s  t  !:■  d  i  r-  •  t.  i  e 

of  '.ran. 77  rm.it  i  n:  ATOR  -  .TRUE,  for  ant'  i.ua- 

t. .  ■  —  r ,  ■  lom*  (Far  i "  i  •  n  ■)  )  :  A'A.-l  .  FAIjSE  .  1  or 
:  ]  .i:,>  -  t  -.ui ;  •  nn.a  tEguut.  t.  h  •)  )  . 

>  Ci!  :  (x  ki  \  AT'. .  a-  i ; '  •  :  in  ;  a  n  ’  i  .in  .'uirij  uh  : 

i .y  nui  » i  -nt  ii,.  PiR [ l  NT. 

:  .  ;  n.il  »-»  -j*  i  s..  :■  rn  ■  ;u i  !  •  d  :  :Jui  a  u t  ;  n.  ORIENT  mi::  t  :  •  •  n  .  M  •  . ;  :  r  i  r 

•  t  i  i:l  .  71  :  .  VRi  TOR  no  Ei.it  T  will  availai  1 .  . 

:  .  For  mi  lin'd,  ;•  •  ;;u!. rout  i  ne  Ol  ■  I  ENT  in  Cii<<j  t  •  r 
•  -  .  Ir'-ii.ir.  il'-w:  -'I'R]  a  Vi  listing  belov:  dir-  tin  to  K  i-.i.i  t  i  •  mi.-  (  i  -  ) 

■  no  ( 
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Chapter  7 


SUBROUTINE  INCPW 

Purpose:  To  compute  the  rectangular  vector  components  of  tin 

electric  field  of  a  ;  lane  electromagnetic  wave  propagat  ing  in  t 

direction  -k  and  incident  on  the  z  =0  plane  of  the  antenna  coor- 
A  A 

dinate  system,  where  ( x  =0 ,  y  -0,  z  -0)  is  used  as  the  phase  urigi 
Usage:  CALL  INCPW  (KA,  El,  IOPT) 

A  r  pi  union  t  s 

KA  -  Real  input  array  of  three  elements  containing 

the  direction  cosines  of  the1  direction  k  from 

A 

whence:  tin  plane  wave  emanates;  e.g.  ,  KA  ( 1 )  =  k 

xA 

KA( J)  -  k  ,  KA(  3)  --  k  where  k  =  x.  k  +  v  k 
yA  zA  A  A  xA  A 

:i  k  .  . 
i’/v 

::i  output  array  of  thr<  •<::  elementr  font  aininn 

t  f  .  '  COPlpOV  of  the  elec!  C  i 

•  '  e  i  d  r:  ,  u  ;r.a  i::  ■■  .  h  *  hat  |  . ! -  1  . 

;  .  s;  it.  \  ■.  ire  .  v.mch  .a  -c.  •  oir 

a  '  :  '  {  '■■■•  ■  '  n-c 

10  ^  vet  -  a,  • 

I(  >P'  1  V-  mu  '  .  Ilpi  :  'll  I  :  - 

'  OPT  jilt  li.un.  i  la  ■  ■  ' 

•;  ;  ■■  •  :  i  • 

i  i :  [si*  ;:,i’  P  !:  is  ;  ,  i  lied  tl  .. 

in...;  •  '1 


.all.  t  i 


£  =  elevation  component 

A 

Ot  =  azimuth  component 


— ^  Qa  >  0 

— ►  e€  >  o 


e  .Coordinate  System  for  Far  Field  Patterns 


r 


will  see  the  tip  of  the  E  vector  tree-  out  a  circle  in  a  plane  of  equal 

phase,  with  the  direction  of  rotation  being  clockwise  for  right-hand 

circular  and  counterclockwise  for  left-hand  circular. 

b.  Spherical  to  rectangular  coordinate  transformations  11]  are 

used  to  define  the  rectangular  vector  components  E  ,  E  ,  E  in  terms  of 

x  y  z 

the  transverse  spherical  components  E  ,  E  .  Let  r=k=xk  +yk  + 

t  a  x  y 

z  k__  represent  the  direction  from  whence  the  plane  wave  emanates,  where 

k  ,  k  ,  k  are  the  direction  cosines  of  r  =  k.  Then,  with,  reference  to 
x  y  z 

Figure  7-1,  there  results 


-k  k 

x  v 


1  -  k 


.  , - -  .  -k  k 

,  ,  Z  V- 

+  v  v  1  -  k  +  z  -  i 


7*= 


i  -  k 


(1) 


1  -  k 


=■•  +  y('O')  + 


d  -  k 


t  - ; 


except  at  k  =  '  1,  where  those  equations  reduce  to 


a  -  x 


(  ■■) 


l 1 ' 


Tile  expressions  for  the  f i ■  Id  components  for  tin  !  sir  on  <  s  of 
inter-St.  are  simnarized  In  Table  7-1 .  The  corrospondi  tic  r  ignetie  field, 
can  -  obtained  from 


H  -  (11  x  k) 


ireiir.im  Fiiw:  i'ompare  expression:-  in  T  ii  b  -  1  .  i .  i  e  t  ; 


r  i- (t  arn  1  i  st  i  tie 


below. 


r1 — - 

AO-A099  162 


GEORGIA  INST  OF  TECH  ATLANTA  F/S  t7/9 

parametric  investigation  OF  RADOME  ANALYSIS  METHODS.  VOLUME  II.— ETC (U) 
FEB  81  6  K  HUOOLESTON.  H  L  BASSETT  AF0SR-77-3469 

AFOSR-TR-81-0460  NL 


UNCLASSIFIED 


SUBROUTINE  INCPW(KA,£I.IOPT > 

KA^NEGATTVE  OF  DIP  OF  PR OP**N  OF  INCIDENT  PLANE  WAVE  (ANT  COORD) 
El*  ElcctpiC  f  IELC  veCTC°  OF  IfiCIOENT  PLANE  WAVE  (OUTPUT) 
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Chapter  8 


SUBROUTINE  RE CM 

8-1.  Purpose:  To  compute  the  complex  voltages  produced  at  the  terminals 
of  the  three  channels  of  a  radome  enclosed  monopulse  antenna  by 
a  plane  wave  of  specified  polarization  and  direction  of  arrival. 
8-2.  Usage:  CALL  RECM  (PWI,  KA,  NX,  NY,  KXMAX,  KYMAX,  FGHZ ,  ROTATE, 

TRANSL,  SUMX,  SUMY,  DELX,  DELY,  VR,  TABLE,  SUPPRS,  RSQMAX) 
8-3.  Arguments 
PWI 

KA 


NX, NY 

KXMAX, KYMAX-  Real  variables  which  represent  the  normalized 

folding  wavenumbers  corresponding  to  the  sample 
distances  Ax  ,  Ay  according  to  Ax  =A/ (2*KXMAX) , 

A  A  A 

Ay  =1/ ( 2*KYMAX) ,  where  A  is  the  free  space 
wave  length. 

FGHZ  -  Frequency  in  gigahertz  of  the  monochromatic  plane 

wave. 


A  complex  array  of  three  elements  containing  E^, 

E  ,  E  of  the  incident  plane  wave.  See  Subroutine 
y  z 

INCPW. 

-  A  real  array  of  three  elements  containing  the 

direction  cosines  k  „ ,  k  „ ,  k  .  of  the  unit 
xA  yA  zA 

vector  kft  which  points  from  the  antenna  origin  in 
the  direction  from  whence  the  plane  wave  emanates. 
The  even  integer  number  of  sample  points  in  xft 
and  y^  directions  used  to  represent  the  antenna 
aperture  fields. 


k* 


05 


ROTATE .TRANSL- 


SUMX,SUMY 


where 


DELX , DELY 


DAZX , DAZY 


VR 


Real  matrices  of  direction  cosines  and  translation 
distances  used  to  carry  out  coordinate  transfor¬ 
mations  of  points  and  vectors  from  antenna  to 
radome  coordinate  systems,  and  vice  versa.  See 
Subroutine  ORIENT. 

Two  dimensional  (NX  X  NY)  complex  arrays  of  the 
x  and  y  vector  components  of  the  antenna  aperture 
fields  for  the  sum  channel  of  a  three-channel 
monopulse  antenna.  The  element  at  I=NX/2+l,  J=NY/ 
2+1,  corresponds  to  that  at  x^=0,  yA=0  -in  the 
aperture.  The  general  correspondence  is  giv^.r.  by 

x  —  —  x  +  (1-1 ) *Ax  =  (I-MIDNX)*Ax 

A  max  A  A 

y  =  -  y  +  ( J-l ) *Ay  =  (J-MIDNY) *Ay 

A  max  A  A 

x  =  Ax  *NX/2  and  y  =  Ay,,  *NY/2. 
max  A  max  A 

Also  see  Subroutine  HACNF. 

Antenna  aperture  fields  for  the  difference  ele¬ 
vation  channel. 

Antenna  aperture  fields  for  the  difference 
azimuth  channel. 

Complex  array  of  three  elements  which  on  output 
contains  the  complex  terminal  voltage  of  the 
antenna  for  the  sum,  elevation  difference,  and 
azimuth  difference  channels,  respectively. 
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TABLE 


Logical  variable  required  by  Subroutine  RXMIT: 


if  TRUE,  a  look-up  table  is  used  to  calculate  the 

transmission  coefficients  of  the  radome  wall;  if 

FALSE,  these  coefficients  are  calculated  exactly 

for  each  angle  of  incidence  specified. 

SUPPRS  -  Logical  variable  used  to  control  the  printing  of 

results  from  Subroutine  RXMIT:  if  FALSE,  a  table 

of  power  transmission  and  reflection  coefficients 

for  equal  increments  in  the  sine  of  the  incidence 

angle  is  printed.  The  phases  of  the  complex 

voltage  transmission  and  reflection  coefficients 

of  the  radome  wall  are  also  printed. 

RSQMAX  -  Real  variable  denoting  the  maximum  radius  of  the 

2  2 

antenna  aperture  such  that  any  point  (x  +y  ) >RSQMAX 

Pi  Pi 

is  omitted  from  the  ray  tracing  and  summation  pro¬ 
cedure  used  to  compute  the  received  voltages  VR. 

8-4.  Comments  and  Method 

a.  Subroutines  Required:  TRACE,  VECTOR,  POINT,  RXMIT,  CAXB. 

b.  Method:  The  voltage  VR  induced  at  the  terminals  of  a  linear 
antenna  by  a  "received"  electromagnetic  plane  wave  E  ,  H  is  given  by 

K  K 

the  Lorentz  reciprocity  theorem  as  [1] 


VR(k) 


C  <j>  (Em  x  H  -  E  x  H  )  *  nda 
J  -T  -R  -R  -T 

S 


(1) 


wheiL-  k  is  the  unit  vector  which  points  in  the  direction  from  whence  the 
plane  wave  emanates  and  wh''ro  E  ,  H  are  the  electromagnetic  fields  of 
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the  antenna  as  produced  on  the  closed  surface  S  which  surrounds  the 
antenna  when  it  is  transmitting.  The  unit  vector  n  is  the  normal  to  S 
pointing  into  the  region  not  containing  the  antenna,  and  C  is  a  com¬ 
plex  constant. 

When  the  closed  surface  S  is  the  z  =o  plane,  n=z  ,  da=dx  dy  ~&x  Ay 

A  A  A  A  A  A 

and  the  integral  in  (1)  can  be  approximated  by 

V  (k)  =  CAx  Ay.  V  V  (E  H  -  E  H  -  E  H  +  E  H  )  (2) 

R  A  yA  ,  '  Tx  Ry  Ty  Rx  Rx  Ty  Ry  Tx 

1  m  ‘ 

where  Ax.,  Ay.  are  equal  sample  distances  in  x  and  y  and  where  the  roc- 
A  A  A  A 

tangular  vector  components  of  the  fields  on  the  z^-0  plane  are  given 
gonerically  by 


E  =  X  E  +  y  Em  +  Z.  E 
— T  A  Tx  A  Ty  A  Tx 


(3) 


It  is  assumed  that  the  fields  E  ,  H  on  S  with  the  radome  in  place 

-T  -T 

are  unperturbed  by  the  radome.  Also,  E^  is  specified  according  to  the 
aperture  distribution  and  polarization  desired  as  is  usually  done  in 
antenna  analysis.  The  corresponding  magnetic  field  H^,  however,  presents 
something  of  a  vexation  in  that  a  non-Maxwel ] ian  aperture  field  can 
result  if  11^  is  specified  independently  of  E^  and  Maxwell’s  equation 
HJr=VxE/r/- jwp .  On  the  other  hand,  specification  of  independently  of 
is  tantamount  to  specifying  magnetic  and  electric  current  sheets  in 
the  antenna  aperture  which  produce  two  independent  solutions  to  Maxwell's 
equations  whose  sum  yields  the  total  fields.  This  latter  approach  is 
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taken  here  when  the  geometrical  optics  approximation 


— T 


z  x  E 
A  — T 


(4) 


is  utilized,  where  t\—v  yi/e  ~  377  ohms.  Also,  the  magnetic  field  Hr  is 
given  by  a  similar  formula 


-R  = 


-k  ^ 


(5) 


where  -k  is  the  direction  of  propagation  of  the  incident  plane  wave. 

Combining  the  results  of  Equations  (4)  and  (5)  into  (2),  and  desig¬ 
nating  the  origin  xA=yA=zft=0  as  the  phase  reference  for  the  complex 
fields,  there  results 


VR(k, 


C‘  l  It  KErp  +  Em  E  }  (1  +  k  )  - 

,  u  Tx  Rx  Ty  Ry  zA 

i  m 


(k  Em  +  k  E  )  E  1  . 
xA  Tx  yA  Ty  Kz 


j  2tt 


<kxAXA+kyA  yA>)(6l 


where 


k 


x  k 
A  XA 


yA  kyA 


ZA  kzA 


(7) 


1 


(8) 


99 


and  where  the  exponential  factor  accounts  for  the  phase  of  the  incident 


wave.  It  is  noted  that  k  ,  k  ,  k  are  direction  cosines  of  k;  hence, 

xA  yA  zA 

kz^=cos6,  where  0  is  the  polar  angle  measured  from  the  z^-axis  in  the  usual 
spherical  coordinate  system.  The  (l+cos0)  term  in  Equation  (6)  is  charac¬ 
teristic  of  the  geometrical  optics  approximation  of  Equation  (4)  [2], 

The  other  factors  have  been  absorbed  into  complex  constant  C'. 

The  effects  of  the  radome  on  the  received  voltage  given  by  Equation 
(6)  are  accounted  for  by  tracing  a  ray  from  each  aperture  element  Ax^Ay^ 
in  the  direction  k  and  weighting  the  field  Er  associated  with  the  ray  by 
the  complex  insertion  transmission  coefficients  T^,  T^  of  the  radome  wall 
as  shown  in  Figure  8-1.  These  coefficients  depend  on  the  incidence  angle 

0.  and  the  plane  of  incidence  defined  by  k  and  the  unit  inward  normal  n„ 

1  R 

to  the  radome  wall  at  each  point  of  incidence  for  each  ray  as  illustrated 
in  Figure  8-2.  The  ray  tracing  is  carried  out  in  the  direction  k,  and 
the  direction  of  propagation  of  each  ray  is  assumed  to  be  the  same  on 
both  sides  of  the  wall,  an  assumption  that  mandates  use  of  the  insertion 
transmission  coefficients  defined  for  an  infinite  sheet  by 


E±  (P) 

Ti  =  E  (P) 
li 


(9) 


E,,  (P) 


En  >  m 


(101 


where  the  numerator  in  each  case  is  the  field  at  point  P  with  the  sheet 
in  place  and  the  denominator  is  the  field  at  the  same  point  with  the 
sheet  removed. 
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Figure  8-2. 


Plane  Wave  Propagation  Through  an 
Infinite  Plane  Sheet 
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The  ray  tracing  is  carried  out  in  radome  coordinates  (x  ,  y  ,  z  ) , 

K  R  R 

and  transformations  of  points  and  vectors  from  antenna  coordinates  (x  ,  yft, 
z  )  to  radome  coordinates,  and  vice  versa,  are  required.  (These  transfor¬ 
mations  are  described  in  detail  in  Subroutines  ORIENT,  POINT,  and  VECTOR.) 
Let  (x  ,  y  ,  0)  be  the  point  in  the  aperture  from  which  the  ray  (line) 

A 

emanates  in  the  direction  k.  Convert  this  point  and  unit  vector  to  the 

radome  coordinate  system.  Find  the  intersection  (x  ,  y  ,  z  )  of  this 

RI  Ri  RI 

/  2  2 

ray  with  the  inner  radome  surface  as  described  by  f(/x  +  y  z)  -constant 

R  R 

since  it  is  a  surface  of  revolution  (Subroutines  TRACE  and  OGIVE).  Compute 

the  unit  inward  normal  n  to  the  surface 

R 


Vf 

IW1 


x  n_+y_n_+z_n_ 
R  xR  R  yR  R  zR 


and  convert  it  to  antenna  coordinates 


(ID 


x,  n  *  + 
A  xA 


y  n 

A  yA 


+  z . 


zA 


(12) 


Use  n  and  k  to  determine  the  plane  of  incidence,  angle  of 

the  transmitted  plane  wave  E '  H‘  (see  Subroutine  RXMIT) 

— R  '  ~R 

Substitute  into  Equation  (f>)  and  sum  the  results  to  obtain 
expression  for  the  received  voltage 


incidence,  and 
for  this  ray. 
the  following 


VR(k) 


c"  l  l  I-U 

1  m 


+  kzA,(ERx 


hTx  +  hRv 


E  ) 
Tv 


(k  ,  Em  +  k  „  Em  ) E 1 
xA  Tx  yA  Ty  Rz 


(k 


xA 


+  K 


yA  '  A 


(13) 


1(13 


where  a  sign  change  and  n  have  been  absorbed  into  C". 


Equation  (13)  with  C"=l  is  used  in  Subroutine  RECM  to  compute  the 

received  voltage  on  each  of  the  three  monopalse  channels.  Note  that  the 

received  field  Ej  ,  H'  at  each  point  (x  ,  y  ,0)  is  tile  same  for  all  three 
— R  — R  A  A 

channels  so  that  three  summations  can  be  carried  simultaneously  to  maximize 
computational  speed.  In  each  summation,  only  the  data  corresponding  to 


Em  for  the  sum, 
Ty 


elevation  difference. 


and  azimuth  difference  channels 


need  to  be  changed.  Note  also  that  Equation  (13)  can  be  rewritten  as 


V  fk)  =  I  X  (h  H'  -  E  '  )  -  E  (n  H’  +  E  '  )]ej  ~  (kxA  XA  +  KyA  VA) 

r  L  L  <■  Tr.  Ry  Rx  Ty  Rx  Ry  1  J 


where  n  ,  n  are  given  by  Equation  (5)  . 

8-5.  Program  Flow 

Line  12:  Initialize  the  ray  counter  NRAY. 

Lines  13-18:  Compute  A  (cm),  k  =2n/A,  Ax  ,  Ay  ,  and  the  midpoint 

o  A  A 

of  the  NX  X  NY  data  arrays  corresponding  to  yft=0 

Lines  21-24:  Set  z  =0=PA(3)  and  precalculate  k  k  and  k  k  . 

A  o  xA  o  yA 

Transform  k  to  radome  coordinates  k  =x„  k  +  y  k  + 
A  R  R  xR  YR 

z  k  in  (^reparation  for  the  ray  tracing. 

K  ZK 

Lines  26-28:  Initialize  the  summations  VR(1),  VR(2),  VR(3)  for  the 


received  voltages  on  the  sum,  difference  elevation, 
and  difference  azimuth  channels,  respectively . 


Lines  30-33:  Iterate  for  each  aperture  point  x^=PA(l);  precalculate 
2 

x  and  k  k  x  outside  the  subsequent  loop  for  y  . 

A  O  XA  A  A 


Lines  34-40:  Iterate  for  each  aperture  point  y^=PA(2) .  Compute 

2  2 

x  +y  =RSQ:  if  point  is  outside  RSQMAX,  omit  from 

A  A 

computation. 

Lines  41-47:  Transform  (x^,yA,0)  to  radome  coordinates  and  trace 
ray  to  radome  inner  surface  to  find  unit  inward 
normal  nR.  If  metal  tip  or  bulkhead  is  encountered 
by  ray,  omit  this  ray  from  computation  of  received 
voltages. 

Lines  48-52:  Transform  nR  to  antenna  coordinates  and  compute  the 

transmitted  plane  wave  PWT=(E‘  ,  E’  ,  E’  ). 

Rx  Ry  Rz 

j  2 71  (k  x  +  k  y  ) 

Lines  53-57:  Compute  phase  PC=e  —  xA  A  yA  A  and  apply 

A 

to  E*  ,  E’  ,  E’  . 

Rx  yK  Rz 

Lines  58-71:  Form  n  H*=E'  x  k  and  use  Equation  (14)  to  add  the 

R  R 

contribution  of  this  ray  to  the  received  voltage  on 
each  of  the  three  channels. 

Lines  72-73:  Increment  ray  counter  and  continue  the  iteration  until 
all  aperture  points  have  been  used.  Upon  completion, 
NRAY  equals  the  number  of  rays  used  in  the  summation 
for  each  received  voltage. 

Lines  74-75:  If  SUPERS  is  false,  write  NRAY. 

RETURN 


Test  Case:  See  Chapter  2. 
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Chapter  9 


SUBROUTINE  TRACE 

9-1.  Purpose:  To  direct  the  ray  tracing  to  find  the  intersection  of  a 
ray  emanating  from  a  point  inside  the  radome  and  the  inner  radome 
surface.  All  dimensions  are  in  centimeters.  Radome  coordinates 
are  implied. 

9-2.  Usage:  CALL  TRACE  (PO,  K,  P,  N ,  METAL) 

COMMON/TRACC/Z  2 ,  Zl 
9-3.  Arguments 

PO  -  Real  input  array  of  three  elements  containing  the 

point  P0( x  ,  y  ,  z  )  from  which  the  ray  emanates. 

o  o  o 

K  -  Real  input  array  of  three  elements  containing  the 

direction  cosines  of  the  ray;  i.e.,  K(l)  =  k^, 

K  ( 2)  =  k  ,  K  ( 3)  =  k  . 

y  z 

P  -  Real  output  array  of  three  elements  containing 

the  point  of  intersection  P(x,  y,  z)  of  the 

ray  and  the  inner  radome  surface. 

N  -  Real  output  array  containing  the  direction 

cosines  of  the  unit  inward  normal  vector  to  the 

radome  inner  surface  at  P(x,  y,  z) ;  i.e., 

N ( 1 )  =  n^,  N  ( 2)  =  n  ,  N (  3)  =  n ^  where 

n  =  xn  +  yn  +  zn  . 
x  y  z 

METAL  -  Logical  output  variable  which  indicates  any 

opaque  surfaces  encountered,  such  as  a  metal 
tip  or  bulkhead:  METAL  =  .TRUE.  indicates  that 
P(x,  y,  z)  lies  on  such  an  opaque  surface. 

l-rttOEux.w  bi^-NOi  blU-iu 
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Z2  -  Real  input  variable  which  designates  the  zr 

coordinate  (cm)  of  the  intersection  of  the  ogive 
section  of  the  radome,  and  the  metal  tip  (if  any); 
must  be  set  in  main  program  prior  to  the  first  call 
to  TRACE. 

Z1  -  Real  input  variable  which  designates  the  zR  coor¬ 

dinate  (cm)  of  the  intersection  of  the  ogive 
section  and  the  bulkhead  of  the  air  frame;  must 
be  set  in  main  program  prior  to  the  first  call  to 
TRACE. 

y-4.  Comments  and  Method 

a.  Subroutines  required:  OGIVE,  TDISK,  BDISK,  OGIVEN,  TDISKN, 

BDISKN 

b.  The  inner  surface  of  the  radome  is  represented  by  three  distinct 
surfaces  as  indicated  in  Figure  9-1:  a  planar  bottom  disk  (bulkhead),  a 
tangent  ogive,  and  a  planar  top  disk  (base  of  a  metal  tip) .  The  ray  is 
traced  to  the  ogive  surface  first  to  find  a  point  of  interso.'Ction  P(x,  y,  z) 

(1)  If  z  <z<z^»  then  the  ogive  section  of  the  radome  was  struck, 
the  unit  normal  is  computed  (OGIVEN) ,  METAL  is  set  to  .FALSE, 
and  the  program  returns. 

(2)  If  z>z^,  it  is  assumed  that  the  ray  encountered  the  top  disk 
before  impincmg  on  the  ogive  surface  (which  actually  extends 
beyond  the  z^  coordinate).  The  ray  is  then  traced  to  find  its 
intersection  with  the  plant;  z  =  r.,( .  If  the  toj  disk  is  indeed 
struck,  then  METAL  is  set  .TRUE,  and  n  =  -z  is  returned. 
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(3)  If  z<z^  from  (1)  above,  it  is  assumed  that  the  bottom  disk 

was  struck  by  the  ray  before  it  encountered  the  ogive  surface. 
The  ray  is  traced  again  to  find  its  intersection  with  the 
plane  z  =  z^.  If  this  bottom  disk  is  indeed  struck,  then 
METAL  =  .TRUE,  and  n  =  z  is  returned. 

The  steps  in  (2)  and  (3)  above  appear  to  be  unnecessary;  however, 
they  are  included  to  ensure  that  the  ray  tracing  procedure  works  correctly 
and  to  alert  the  user  if  it  does  not.  For  example,  if  incorrect  variable 
values  are  passed  to  the  supporting  subroutines,  there  is  a  good  chance 
that  no  intersection  will  be  found  with  any  one  of  the  three  surfaces,  in 
which  case  the  following  error  message  is  outputted: 

"THERE  IS  A  HOLE  IN  THIS  RADOME" 

The  message  is  continued  with  the  values  of  (x  ,  y  ,  z  )  and  (k  ,  k  ,  k  ) . 

o  ooxyz 

Incorrect  values  of  geometry  variables  passed  to  the  supporting  subroutines, 
and  attempts  to  trace  a  ray  from  a  point  exterior  to  the  inner  radome 
surface,  will  prompt  the  error  message  and  alert  the  user  of  his  mistake. 
9-5.  Program  Flow:  See  listing  below. 

9-6.  See  Chapter  2. 

9-7.  References:  None 

9-8.  Prognm  Listing:  See  following  pages. 
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Chapter  10 


10-1 


10-2 


10-3 


SUBROUTINE  RXMIT 

Purpose:  To  compute  the  complex  rectangular  vector  components  of 

the  electric  field  E  transmitted  through  the  radome  wall,  where 
-T 

it  is  assumed  that  the  incident  field  E.,  H.  is  locally  a  plane 
wave  and  that  the  radome  wall  behaves  as  an  infinite  plant-  dielec¬ 
tric  sheet.  The  direction  of  propagation  of  the  plane  wave  -k 
and  the  unit  inward  normal  n  at  the  point  P^(x,  y,  z)  are  used  to 
determine  the  angle  of  incidence  and  the  plane  of  incidence  of 
the  plane  wave.  All  dimensions  are  in  centimeters.  Radome  coor¬ 
dinates  are  implied. 

Usage:  CALL  RXMIT  (PWI,  PWT,  K,  NORM,  PI,  TABLE,  SUPERS,  BETA) 

COMMON/TRANSC/DIN  ((•)  ,  ER(6),  TL>(h),  TZ ,  WALTOL,  N ,  NN , 

D ( b) ,  ZB,  TK 
Arguments 

PWI  -  Com|  lex  input  array  containing  the  vector  com¬ 

ponents  of  tin  incident  electric  fit-id;  i.e., 

PWI  (T  .  ,  E  .  ,  E  .  )  . 

xi  y  i 

PWT  -  Com|  lex  out | ut  array  containing  the  vector  com¬ 

ponents  of  tin  transmitted  electric  field;  i.e., 

PWT  ( E  .  E  ,  E  ) . 
x  t  y  t  z  t 

K  -  Real  input  array  containing  the  direction  cosines 

of  tin  direction  from  whence  the  plane  wave 

emanate;;;  i.e.,  K(k  ,k,k)=xk  +  v  k  +  k 
x  y  z  x  y 
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N,  the  exit  side  of  the  wall;  layer  N  is  the  first 

NN  layer  on  the  incident  side.  ER(NN) ,  TD(NN) 

specify  e^_,  tanfi  of  the  medium  in  which  the  N-layer 
structure  is  immersed  (normally,  free  space  so  that 
ER(NN)  =  1.0,  TD ( NN)  -  0.0).  The  real  array  D 
contains,  after  the  first  call  to  RXMIT,  the 
thickness  in  centimeters  of  each  layer. 

TZ,  -  Real  variables  used  previously  to  specify  longi- 

WALTOL,  tudinal  and  circumferential  variations  in  wall 

ZB,  TK  thickness  and  in  the  tolerance  on  thickness.  These 

variables  are  not  active  in  this  version  of  RXMIT. 

10-4.  Comment  and  Method 

a.  Subroutines  required:  WALL,  AMPHS,  AXB 

b.  The  transmission  of  an  incident  plane  wave  through  a  plane 

-12 

dielectric  sheet  immersed  in  free  space  (e  =  8.854  X  10  farads/m,  ii  = 

o  o 

-7 

4tt  X  10  henries/m)  may  be  described  in  terms  of  the  insertion  voltage 

transmission  coefficients 


T 

l 


E.  (P') 
11 


(1) 


where  E^ ,  are  the  transmitted  fields  at  P'  with  the  sheet  in  place, 

and  E.  ,  E. are  the  incident  fields  at  the  same  point  in  the  absence  of 
ii  ill 

the  sheet.  The  point  P'  lies  on  the'  colinear  extension  of  the  incident 
ray  and  is  located  on  the  exit  side  of  the  sheet. 
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'I',,  ari-  di  f  ferent,  it  is 


Simv  thf  transmission  coefficients  T  , 
necessary  to  resolve  the  incident  electric  field  E.  into  jh 

—  l 

anil  |  ar  il  h  1  comi  onents;  i.e.,  vector  components  which  are 
to  and  I  aiiil  l'-l  to  tin’  [  lane  of  incidence  (POI)  defiiud  1  >y 
illustrated  in  Figure  lu-1.  The  unit  vector  per] end icu 1  a r 

Cl  Veil  by 


k  x  n  k  x  n 

R  _  _  R__ 

|k  x  nR|  sin|k,nR 


A  unit  vector  [at  all'.  1.  to  the  Pol  is  qiveti  by 


k ,  |  a  x  k 

II  1 


Th  incident,  electric  field  may  be  written  as 


IF  -•  x  K  .  +  y  I!  .  + 
■  i  xi  yt 


11.  *  k  E.  +  k 

iz  t  1  J 


whe  r  i 


i: .  -=  k  •  E.  =  k  E  .  +  k  E  +  k  E  . 

it  1  ~i  xl.  xi  yt  yi  zi  l 


h  ii "  k 


F  -  k  M  F  +  k  I,  K  .  +  k  , .  E 

-J  x  I i  xi  y  II  yi  z||  z 


and  win.  t'  k  ,  k  ,,,  etc.  are  the  vector  com]  >onent  s  of  k  ,  k 
x  L  x  1 1  1 


>f  the  coordin.it. e  system  (x,  y,  z)  , 


rpendicu  1  ar 

I erpendi cul ar 

k  and  n  as 
R 

to  the  POI  is 


( i) 


(4) 


(h) 


(6) 

(7) 

In  terms 


E 

— 1 


X  ( E  .  +  E  . 

Xll  XI 


+  y (E  .  +  E  . 

yii  y  i 


+  z(E  .  +  E  . 

7.11  71  II 


wile  re ,  for  example 


(8) 


E  . 
Xll 


k 

xi 


E 

li 


The  transmitted  plane  wave  is  then  given  by 


(9) 


E  =  k 
-T  1 


E.  +  k 
il 


(10) 


r: 

~T 


x  ( T  E  . 

1  Xll 


+  T„  E 


ll>  +  V(T, 


E  .  +  T. 

yil  I 


yi 


,)  +  z(T 


Zll 


+  T, 


(ID 


10-'). 


Program  Plow 
Lilies 
Line  9: 

Lines  10-12: 


Lines  15-16: 
Lines  17-59: 


Comments 

Set  NANGLE  =  number  of  entries  used  in  the  )ook-u)> 
tables  for  T  ,  T... 

.1  II 

NDO  causes  initialization  of  variables  and  the  com¬ 
putation  of  the  look-up  tables  on  the  first  call  to 
RXMIT  (lines  11-59). 

Convert  layer  thicknesses  from  inches  to  centimeters. 
Compute  look-up'  tables  for  ,  T  j  |  at  NANGLE  points 
spaced  equally  in  sint).  over  the  range  (0,  1)  .  If 
SUPPRS  -  .FALSE.,  print  a  table  of  transmission  co¬ 
efficients  (every  fifth  point  only).  If  ER(1)<  1.05, 

set  AIR  -  .TRUE,  and  compute  unity  transmission  coeffi¬ 


cient:;  for  the  "air"  railome  (for  testing). 


Lines  60-78:  Compute  sin0^. 

Lines  79-86:  Interpolate  in  table  to  compute  T^,  T| | a t  sin© . . 
Lines  87-100:  Normalize  the  vector  k  . 


1 


Lines  101-112 :  Compute  E 
Lines  113-124 :  Compute  E 
Lines  125-129 :  Compute  E 
Lines  130-136:  If  sin0 . 

i 

message , 
Test  Case:  See  Chapter 
References:  None 


.  ,  E  .  ,  E  .  . 

xri  yil  zii 

xiir  Eyiir  Eziir 

xt'  Eyt'  Ezt  and  return- 

is  out  of  range  of  the  table,  write  error 
set  ,  T j |  to  unity,  and  continue. 

2. 


Program  Listing:  See  following  pages. 
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Chapter  11 


SUP  ROUTINE  WALL 


11-1.  Purpose:  To  compute  the  transmission  and  reflection  coefficients 

of  a  N-layer  dielectric  sheet  having  thicknesses  d  ,  dielectric 

n 

constants  e  ,  and  loss  tangents  tan6  for  each  layer  when  a  plane 
rn  n 

wave  is  incident  at  angle  0  . 

11-2.  Usage:  CALL  WALL  (BETA,  SINE,  D,  ER,  TD,  N,  NN,  TPER ,  TPAR, 

RPER,  RPAR) 


11-3.  Arguments 
BETA 


SINE 

D, 


ER, 


TD 

N 


NN 

TPER, TPAR 


RPER, RPAR 


Real  input  variable  =  2n/A,  where  A  is  the  free 
space  wavelength. 

Real  input  variable  =  sin  0. . 

l 

Real  input  arrays  containing  the  thickness  (cm) , 
dielectric  constant  e  ,  and  loss  tangent  tan6  of 
each  layer. 

Integer  input  variable  equal  to  the  number  of 
layers . 

Integer  input  =  N+l . 


Complex  output  variables  equal  to  the 
voltage  transmission  coefficients  for 
[onenfs  of  the  incident  electric  field 
eular  to  and  parallel  to  the  plane  of 
respect i ve ly . 

Coin]  lex  output  variables  equal  to  the 


coefficients  R^ ,  Ry 


i nsertion 
the  oom- 
perpendi- 
i ncidenco , 

reflection 


r ru  e r. .  i  fiLe-ii; 


139 


11-4. 


Comment  and  Method 


a.  Layer  1  is  the  first  layer  on  the  exit  side  of  the  panel;  layer 
N  is  the  first  layer  on  the  incident  side.  T  ,  T„  have  the  same  value  for 

1  II 

either  side  of  the  panel  being  the  incident  side;  however,  ,  Rj,  are 
different  (in  phase)  for  the  two  cases. 

b.  The  details  of  the  method  are  presented  in  Reference  1  and  are 


reproduced  in  Appendix  E. 

11-’.  Program  Flow:  See  Reference  1. 

11-b.  Test  Case:  See  Chapter  2. 

11-7.  References 

1.  K.  B.  .Toy  and  G.  K.  Huddleston,  "Radome  Effects  on  the 
Performance  of  Ground  Mapping  Radar,"  Technical  Report, 
Contract  DAAHol -72-C-Ubdtf ,  U.  S.  Army  Missile  Command, 
Mo.  el.  1  '7  1. 

i  1  -••••.  Program  List  ir.g:  following  g  ag<  s. 
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Chapter  12 


SUBROUTINE  AXB 

12-1.  Purpose:  To  compute  the  real  vector  cross  product  C  =  AX  B,  whore 
A,  B,  C  aro  expressed  in  rectangular  congxments . 

12-2.  Usage:  CALL  AXB (A,  B,  C) 

12-3.  Arguments 

A,  B  -  Real  input  arrays  containing  the  rectangular  com¬ 

ponents  of  A  -  x  A  +  y  A  +  z  A  and  B;  i.e-, 

—  x  J  y  z  — 

A(A  ,  A  ,  A  ) ,  B(B  ,  B  ,  B  ) . 
x  y  z  x  y  z 

C  -  Real  output  array  containing  the  rectangular  com¬ 

ponents  of  the  vector  C=AXB;  i.e.,  C(C  ,  C  ,  C) 

—  —  —  x  y  z 

12-4.  Comment  and  Method 

a.  Both  input  vectors  must  be  real. 

b.  The  computation  of  C  =  A  X  B  is  elementary. 

12-5.  Program  Flow:  See  listing  below. 

12-6.  Test  Case:  None 
12-7.  References:  None 

12-8.  Program  Listing:  Sec  following  page. 
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Chapter  13 


SUBROUTINE  CAXB 


13-1.  Purpose:  To  compute  the  complex  vector  cross  product  C  =  A  X  B, 
where  A  is  a  complex  vector  and  B  is  a  real  vector  expressed  in 
rectangular  coordinates. 

13-2.  Usage:  CALL  CAXB  (A,  B,  C) 

13-3.  Arguments 


13-4. 


13-5. 


13-6  . 
13-7. 
13-8. 


A 


B 


Complex  input  array  containing  the  rectangular 

components  of  the  vector  A=xA  +yA  +  z  A  ; 

—  x  y  z 

i  .e .  ,  A  (A  ,  A  ,  A  )  . 

x  y  z 

Real  input  array  B  (B  ,  B  ,  B  )  representing 

x  y  z 


the  vector  B. 

C  -  Complex  output  array  C  (C  ,  C  ,  C  )  representing 

x  y  z 

the  vector  C  =  A  X  B. 

Comment  and  Method:  None 
Program  Flow:  See  listing  below. 

Test  Case:  None 
References:  None 

Program  Listing:  See  following  page. 
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Chapter  14 


SUB  ROUT  INI'  RECBS 

14-1.  Purpose:  To  compute’  the  anglt  of  arrival  k  of  a  plane  wave  on  a 
monopulse  antenna  which  yields  an  electrical  boresight  indication 
which,  due  to  the  radome ,  may  be  different  from  the  mechanical 
boresight  along  t  ho  2  axis  of  the.’  antenna.  The  antenna  aperture 
lies  in  the  x y  plane.  All  dimensions  are  in  centimeters.  Antenna 
coordinates  are  implied. 

14-4.  Usage:  C'ALl-  RECBS  (SUMX,  SUMY,  DEL>  ,  DELY,  DAZX,  DAZY ,  NX,  NY, 
UMAX,  NS,  I OPT ,  VR,  DMRAD ,  ROTATE,  TRANSL ,  FGHZ,  KXMAX , 
KYMAX,  TABLE,  SINGS,  K,  AZTM,  ELTM,  RSQMAX,  VMAX,  SMAX, 
SUPPRS) 

14-3.  Arguments 

SUMX,  SUMY,-  Complex  input  arrays  of  NX  by  NY  elements  each 

DE1.X,  DELY,  containing  the  aperture  distributions  of  the 

DAZX,  DAZY  monopulse  antenna.  See  Subroutine  HACNI-’. 

LMAX  -  Integer  input  variable  which  controls  the  maximum 

number  of  iterations  that  will  be  done  to  find 
the  electrical  boresight  within  the  tolerance 
specified  by  DMRAD. 

NS  -  Inactive  integer  variable. 

I OPT  -  Integei  input  variable  which  selects  the  polari¬ 

zation  of  thi  incident  plane  wave.  See  Subroutine 
INi'.PW  . 

VK  -  Complex  array  of  three  elements  representing  the 

received  voltage  on  the  sum,  elevation  difference, 


13D 


and  azimuth  difference  channels  of  the  antenna. 


respectively  (output) . 

DMRAD  -  Real  input  variable  equal  to  the  tolerance  to  which 

the  electrical  bore-sight  in  milliradians  is  to  be 
computed;  i.e.,  0.1  miiliradian. 

ROTATE,  -  Variables  required  by  Subroutine  RECM. 

TRANS L,  See  Chapter  8. 

FGHZ , 

KXMAX, 

KYMAX, 

TABLE 

SINOS  -  Real  variable  equal  to  the  sine  of  the  angle  0 

os 

(measured  from  the  z-axis)  in  the  <f>  =  45°  plane 
(<{>  measured  from  +x  toward  +y)  at  which  the  first 
target  return  arrives;  e.g.,  0  =  3  degrees. 

K  -  Real  output  array  containing  the  direction  of 

arrival  of  the  final  target  return;  i.e., 

K(k  ,  k  ,  k  )  . 
x  y  z 

AZTM,  -  Real  output  variables  equal  to  angles  (mrad)  in 

ELTM  the  azimuth  and  elevation  planes  of  the  antenna 

which  specify  the  direction  of  arrival  of  the  final 
target  return.  If  k  is  the  unit  vector  fiointing 
from  the  origin  in  the  direction  of  the  final 
return,  then  the  orthographic  projection  of  this 
vector  onto  the  xz-plane  makes  an  angle  AZTM  with 
the  z-axis;  its  projection  onto  the  yz-plune 
makes  an  angle  ELTM  with  the  z-axis. 
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RSQMAX  -  Real  variable  needed  by  Subroutine  RECM.  See 

Chapter  8. 

VMAX  -  Unused. 

SMAX  -  Real  output  variable  equal  to  the  magnitude  of 

the  received  sum  voltage  for  the  final  return; 
used  to  compute  loss  in  antenna  gain. 

SUPPRS  -  Logical  input  variable  which  controls  the  compu¬ 

tation  and  printing  of  additional  antenna  outputs 
around  the  boresight  direction.  If  TRUE,  the 
complex  voltage  outputs  of  the  difference  channels 
will  be  computed  at  one  milliradian  increments  over 
the  range  13.0  mrad,  centered  on  the  direction  of 
the  final  target  return. 

14-4 .  Comments  and  Method 

a.  Subroutines  required:  AMPHS,  RECM,  INCPW. 

b.  Subroutine  REC3S  uses  a  linear  tracking  model  to  determine  the 

direction  of  arrival  k-xk  +vk  +xk  of  a  plane  wave  which  will 

x  '  y  z 

produce  null  indications  in  the  elevation  and  azimuth  difference  channels 
of  the  monopulse  antenna  inside  the  radome.  Subroutine  RECM  is  used  to 
compute  the  received  voltage  on  each  channel  for  the  specified  polarization 
(IOPT)  and  direction  of  arrival. 

The  first  target  return  is  made  to  arrive  from  the  direction 

'  r~  2 

k  =  x  sinO  +  y  sin0  +  z/l-2sin  0  (1) 

1  os  os  os 

to  produce  out.  juts 
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U  .  -  I  in  [^l 
AZl  I  Vy  \ 


'  ael  f 

E  ,  =  In.  , - ' 

ELl  I  Vj.  | 


wiK  re  V..,  V  ,  V  an-  t  he  complex  voltage  outputs  of  t!u  three- 
;vEL  a A/, 

channel  v»ut|  ut«  oi  t  ho  three-channel  anti  nna.  The  second  return  is  made 
to  .irrivi  from 


k  =■  x(-sinO  )  +  y(-sint‘  )  +  zv  1  -dsin^t) 

A  os  os  os; 


to  j roduco  out|  uts  U  ,  U 

AZe  EL2 

Construct  a  linear  model  for  each  channel  independent! y  using 
tlu  slopo/i nter. 'opt  equation  for  a  line;  i . e . , 


U  =  M  k  +  b 
AZ  Ani  x  AZ 


U  =  M  k  +  b 
El,  EL  v  EL 


■M  -  (I)  -  U  ,)/  (k  -  k  ) 

AL  Al.  1  AZ2  xl  x2 


M  “  (U  ,  -  U  )/(k  ,  -  k  ) 
EL  ELI  EI.2  yl  y2 


b  =  U  ,  -  M  k 
AZ  AZ  1  AZ  xl 


1),  ,  ■  UT„  ,  -  M  k  , 

EL  ELI  EL  yl 


X,  *• 


p 


Use  tins  nodi' 1  to  estimate  thi1  valm  s  of  k  , 

x 

-  I.';  i  .i-.  , 


k  that  will  make  U 
v  AZ 


U 


EL 


k 

x 


b  /M 
AZ  AZ 


(7a) 


k 

V 


b  /M 
EL  EL 


(7b) 


where  the  value  of  k_  follows  from 


k“  +  k^+k~  =  l 
x  y  z 


(8) 


The  third  target  return  is  made  to  arrive  from  this  direction  and 
the  values  of  U  ,  U  are  computed  via  Subroutine  RECM.  Now,  according 

JL 

to  the  last  linear  model,  a  value  of  U  within  the  range 


AZ  1 


to  1 


(9) 


would  indicate  that  the  null  has  been  found  within  the  tolerance  0 

to  1 

(=DMRAD)  specified.  rf  this  tolerance  is  not  satisfied  for  both 
channels,  then  the  process  is  repeated  until  it  is  or  until  LMAX  is  ex¬ 
ceeded.  In  continuing  the  iterations,  kl  becomes  k,  and  the  estimated 

Z. 

point  becomes  k.  . 

On  the  last  return,  tin  direction  of  arrival  in  specified  by  k. 
The  angles  in  the  azimuth  and  elevation  planes  are  given  in  milliradians 

by 


AZTM 


y 


]  iino 


(10) 


14  3 


1 


Li ru  s  lt.-  ji):  Coin:  ute  first  two  target  returns  to  construct 

linear  tracking  model. 

Lines  11-3H;  Ooitij  uto  sioivs  M  -  SLI’AZ,  M  ~  SLPEL  from  first 

AZ  EL 

two  returns. 

Line  3'»;  Iterate  on  linear  model  up  to  LMAX  times. 

Ianet;  41-44:  If  tlu  increment  in  Ak  is  larger  than  sin  (DMRAli/ 1000 ) 

then  use  it  to  compute  slopes;  if  not,  use  the  last 
conj  uted  values  of  slopes  to  avoid  division  by  too 
small  a  number. 

bines  4  -4<- :  Compute  mt  ercet  ts  b  ,  L 

AZ  EL 

Lines  47 — is;  Corn  ut>  accuracy  criteria  base  d  on  current  slopes 

and  intercepts. 

Lines  4  '-cl:  Confute  direction  k  that,  the  men  If  i  indicated  will 

;  reduce  nul  ls  P  ri~n,  i\  ~0  in  both  planes . 

AL  1,1 1 

bines  comiute  r  ,  f  for  tins  direction  k . 

EL  AZ 

Ppdats  t  in  linear  tracking  model  bv  storing  the  last 
two  ints  in  each  channel  as  U(l),  tJ  ( 2 )  ;  e  .  g  .  , 
e  t  I  i  i;  (2)  and  U  (2)  -  ,  Kj  (1)  =  K.,(l)  and 

K  ( l )  K ( 1 ) ,  etc . 

Lira  '1:  At  b  ,i/.t  t  nr.--  iterations  are  always  used . 

Lines  c, :  it  p  ,  v  are  wit  hin  error  bounds,  exit  the  loop; 

AZ  Li. 

if  not,  continue  to  iterate. 

Lines  “-<7:  If  LMAX  exceeded,  inform  the  user. 

Li:.»  Comiute  am|  I  i  tude  LMAX  on  sum  channel  for  final 

t  .i re"!  return. 

Lines  1  ■- I  :  Ci  input  i  slopes  for  final  n  turn. 

Lm-  s  I  Corn;  at  >  bore;;  j  ,;lit.  error  AZTM,  ELTN. 
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Chapter  15 


SUBROUTINE  RECPTN 

15-1.  Purpose:  To  compute  the  receiving  patterns  of  a  monopulse  antenna 
at  NREC  points  in  a  specified  principal  plane.  A  plane  wave  of 
specified  polarization  (ICOMP)  is  made  to  be  incident  on  the 
antenna  at  equal  increments  in  sine  over  the  range  (-KMAX,  KMAX  -DK) 
in  either  the  elevation  plane  (ICUT  =  1)  or  azimuth  plane.  The 
received  voltage  in  each  channel  is  computed  in  the  presence  of 
the  radome  and  stored  for  return  to  the  calliny  program. 

15-2.  Usage:  CALL  RECPTN  (SUMX,  SUMY,  DE  ',  DELY,  DAZX,  DAZY, 

NX,  NY,  ICUT,  ICOMP,  KMAX,  ’C,  VREC,  KXMAX,  KYMAX, 

FGHZ,  ROTATE,  TRANSL,  TABLE,  SUP'i'RS,  RSQMAX) 

15-3.  Arguments 

SUMX,  SUMY,-  Complex  input  a.  lays  of  NX  bv  NY  elements  con- 
DELX,  DELY,  taininq  the  aperture  field  uistributions  of  the 

DAZX,  DAZY,  monopulse  antenna.  See  Subroutine  HACNF. 

ICUT  -  Integer  input  variable  which  specifies  the  prin¬ 

cipal  plane  in  which  the  pattern  is  computed: 
elevation  (ICUT  =  1)  or  azimuth  (ICUT  =  2)  . 

ICOMP  -  Integer  input  variable  which  specifies  the  linear 

polarization  of  the  incident  plane  wave:  elevation 
component  e  only  (ICOMP  =1)  or  azimuth  component 
a  only  (ICOMP  =  2).  See  Figure  15-1  for  further 
clari f ication . 
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elevation  component 
azimuth  component 


KMAX 


He a  1  input  variable  equal  to  sinO 


,  where1  the 


max 

pattern  is  computed  over  the  angular  range  (-0  , 

max 

i'  ),  but  in  equal  increments  in  si  mi  so  that 
ma  x 

Fourier  interpolation  can  be  applied  directly  in 
the  wavenumber  domain  using  the  Fast  Fourier 
Transform . 

NREC  -  Integer  input  variable  equal  to  the  number  of 

points  at  which  the  pattern  is  computed. 

VREC  -  Complex  output  array  of  NREC  by  3  elements  con¬ 

taining  the  computed  receiving  patterns  for  the 
sum,  elevation  difference,  and  azimuth  difference 
patterns  of  the  monopulse  antenna. 

KXMAX,KYMAX, -  Input  variables  required  by  Subroutine  RECM.  See 
FGHZ, ROTATE  Chapter  «. 

TRANS L, TABLE, 

SUPPRS , RSQMAX 
15-4.  Comments  and  Method 

Subroutines  INCPW  and  RECM  are  used  to  compute  the  incident  plane 
wave  and  the  received  voltage  in  each  channel  for  each  direction  of  arrival 
in  the  specified  plane.  For  Liu  elevation  plane,  the  direction  of  arrival 
is  given  by 


k  =  x  ( 0 )  +  y 


*  r~~ t 

!  si  nO  +  z  r 1 —sin 


:i) 


where  ■■  is  the  angle  measured  from  the  z-axis.  For  the  azimuth  plane 


X  si  111:  +  V  (0  ) 


+  z/l  -sin^O 


1 3  3 


k 


(2) 


^HCo^J,LT>\£)fsw0C(JNC^*r*4fvJ^D  ,f  lJ>  vD  rs.  cC  <X  O^CNj^O  ^LTisDN-ODCTCDT^OJ^O^LA^N*® 

»-|HHHH^Hr4^THC\jMf\J^C\JN(\jMf\Jf\Jn^P0r0K)r0^f0^ 


Chapter  16 


SUBROUTINE  OGIVE 

16-1.  Purpose:  To  solve  for  the  intersection  PH(x,  y,  z)  of  a  line  (ray) 

and  a  tanqent  ogive  surface.  The  ray  starts  at  point  po(x  ,  y  ,  z  ) 

o  o  o 

and  travels  in  the  direction  K(k  ,  k  ,  k  )  =  k  =  xk  +•  yk  +  zk  . 

x  y  z  x  y  z 

Dimensions  are  in  centimeters.  Radome  coordinates  are  implied. 

16-2.  Usage:  CALL  OGIVE  (PO,  K,  PH,  HIT) 

COMMON/OGIVC/RF ,  BSQ,  AP ,  RINV,  B,  RSQ1,  RP2 
16-3.  Arguments 

PO  -  Real  i nput  array  containing  the  point  of  origin 

of  the  ray  PO(x  ,  y  ,  z  ) . 

o  o  o 

K  -  Real  input  array  containing  the  direction  cosines 

of  the  lay  K (k  ,  k  ,  k  ) . 

x  y  z 

PH  -  Real  output  array  containing  the  point  of  intersec¬ 

tion  i'll  (x,  y,  z)  ,  if  HIT  =  .TRUE. 

HIT  -  Logical  output  variable  which  indicates  if  an 

intersection  solution  was  found  (TRUE). 


The  following  variables  are 

common  witii  the 

main 

program  and  are  pirecal- 

eulated  to  speed  up 

the  ray 

1 1  ae  > 

iki  com;  i  it  at  ions. 

Refer  to  Figure  16-1  of 

the  radome  geometry 

for  the 

■  let  init  ions  of 

k  and 

B. 

KP 

-  Real 

11}  'Ut 

v.u  l  abl  i  • 

> 

-  R“  - 

2 

•  B  . 

Brio 

-  Real 

i  n|  ’Ut 

var l able 

*-  B". 

AP 

-  Rea  1 

input 

vat i able 

=■  0. 

See  APIN  in  Section  2-4. 

RINV 

-  Rea  1 

input 

variable 

=  1/H. 

B 

-  Real 

i  nput 

variable. 

See 

Figure  16-1 . 

T.7 


INCIDENT 
PLANE  WAVE 


Figure  1.G-]  .  Tangent  Ogive  Radome  Geometry. 


KS'jl  -  kt--.il  in;  ut  variable  =  R~ . 

RP2  -  Real  input  variable  -  k‘  +  B“. 

6-4.  bomment  and  Metno.l 

a.  The  common  variables  must  be  commuted  in  the  main  program 
in,  i  to  t.he  fiist  call  to  o  'I  Vi: . 

b.  .  ubrout  irit" :  to-jun  od:  CbKT  ,  SyR,  XY  .  Heal  '"ui.cIki.  „'bRT 
am;  ut.es  cuot'  root?  Svk  computes  square  root,  with  test,  for  negative 

t  ■  j union t . 

The  intersection  of  a  ray  and  ogive  sui race  requires  the 
olution  of  a  quartic  equation  in  the  parametei  z  as  follows  [1J 


'hei  e 


2  t  1  iH)  ( JA+V) 

•1  "  2 

(1+U) 

;  ')  1 

+  A  +W)  +  ( 2A+V)  -  413  “li 

i _ _ _ 

2  ’ "  . 

U+U) 

tV)  t-k  1  A"  +W  )  -  4H  *"V 


(  3) 


2  <  1  ♦ i: )  t 


>  K 


an.i  t;  ,iro  defined  or.  F-'i-Turo  16- J  and  by  the 


ax  1  ;  o*  revolut  ion.  !.>:  t'ui  cki  i  Vo  tshuj.o.  T 
i  ;ir.  o:'t  sol  .1  Lon<j  of  tin-  .  00;  ■!  i  nutor  for 
•*'  ;.j  1 1  when  tiio  lol  .lowing  equation.-;  for 

the  point  no  ox  ,  v  ,  )  in  tin-  direct  ion  k 


;  n  t  •  :  ;.rt  u  a.-  (]■’): 


Cot.:  tOIlt 


Me-  mil  tic  i-.tu  it  iiiti  itv.iv  ho  Found  t  rom  th 


v''l 

Cd) 

(d) 

O' :  ive 

U  - 1 ) 

t.he  O'ti’.’i 

n  r  ay 

x  k  4 

x 

d;' 


t  -fr* 


(29) 


win  re  the 
to  (x,  y, 
fol lowing 

t  ”  +  2 

The  quaii r. 


Tin  '.in:*  ’ 
1  -  .  :  1 

L  i .  .  ,  • 

l  -  . 


y  =  y  +  k  t 

o  y 


parameter  l  is  the  distance  along  the  line  from  (x  ,  y  ,  /.  ) 


z) .  Substituting  Equations  (27)  -  (29)  into  (10)  yields  the 


quadratic  equation  in  t 


>2  /  2  2  2  2 

Ik  x  t-  k  v  )t  +  (:<  +  v  )  (  /R  -  (z  -a  )  -  B)  =  R 

X  i  t)  •<>  o  \  s 

(30) 


tic  f.'iTjnJ.i  yi<  Ids  tiic  following  solutions  to  the  above  equation 


2  2  2  U 

-.X  *  k  v  )  -  (x  +  y  )  +  K  ] 

x  o  o  o  s 


:  :  i.i'  ions  (24)  through  (26). 


x  uui  compare-  directly  to  the  above 


i  .  .  a  .  .i,,  "Kadome  Effects  on  the 

rb  :  •:  •  :  Mai  j  l-  Radar,"  Technical  Report , 

Cent  ra  '  A/-.!!  1-7.  -  -  •  ,  p.  5.  Army  Missile  Command, 

Mar cn  1  •?  '•  . 

Steiju:.,  and  Abromowitz ,  Handbook  of  Mathematical  Functions, 
National  Hiii.mu  of  .standards,  June  1964,  p.  17. 

See  following  pages. 
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10-a.  Program  Fisting: 


*  * 


C:ui|  t  <  1  17 


17-..  . 

17-1. 


17—  i . 


win  re 
l  nw.irii 


.SUBROUTINE  (H 1 1  VKN 


Pur:  n'Si' :  To  compute  the  unit  inward  normal  vector  n  -  x  n  + 

x 

v  r:  +  x  n  to  the  tangent  ogive  surface  at  the  point  Pl(x,  y,  z)  . 
Dimensions  are  in  centimeters  and  radome  coordinate:,  are  im; lied. 
Usage:  CALL  OGIVEN  (PI,  N) 

COMKON/OGIVC/RP,  BSQ,  AP,  RINV,  B,  RS<p 1 ,  RP7 
(Sue  Chaj-ter  lb  for  common  variable:;.) 

Arguments 


PI  -  Real  input  array  containing  the  point.  PI(x,  y,  z) 

on  the  tangent  ogive  surface  at  which  the  unit 

normal  is  desired,  as  computed  by  Subroutine  OGIVE. 

N  -  Real  output  array  containing  the  direction  cosines 

(n  ,  n  ,  n  )  of  the  unit  inward  normal  vector, 
x  y 

Comments  and  Method 

Tin  tangent  ocsivo  surface  is  descri  bed  by 


f  (  r , ) 


+  B  =  0 


(1) 


rs~ '. 

\  x  fy 


and  wher<  k  and.  R  ar-  def  i  ned  in  Figure  lo-l. 


Tiie  unit 


nnrm.i  1  t'  tin  s  surface  is  given  by 


n 


df  dr 
x  -  -  -  -  + 

dr  d x 


df  dr 
dr  dy 


df 

dr. 


l(v> 


17-7 


17-tf 


Reference 

1.  Smail,  L.  L.,  Analytic  Geometry  and  Calculus,  Appleton- 
Con  tury-Crof  ts  ,  Inc.,  New  York,  11J'j3. 

Program  Listing:  See  following  page. 
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.sihwh'ttni:  xv 

'ur post.’ :  To  compute  tin.-  x  an<i  y  c.  nir.ii  imi  t  ..  <  •:'  i : j t .  rsect  ion 

i-int  1*1  (x,  y,  /.)  of  a  lis;.  (ray)  haVii:>>  -i:  n  :  ■  •.  i  m  • 

yk  ,  k  ,  k  )  with  a  siurfao-  of  rovoluti ;  n  when  .•  ; ..  ki..  v.  i . .  Tin 
x  y  :: 

me  ;  .issis  tnri'iMi,  t  in  know:.  moint  p(x  ,  v  ,  ).  All  .  i  1  non:;  i  <  <ns 

<  :  (  ; 

:v  in  r.t  i  n.  1 1.  rn  . 

'n  .in-  ■ :  .'Ail,  MY  (•',  K,  1 1 ,  PI) 

,r>iu:ni  nls 

i ■  -  Ken  i  i  n(  ut.  array  coiitaini  ikj  the  known  point 

t  nrou'jh  wn  i  oh  tin  ray  pusses;  i  ,  r  (x  ,  . 

K  1-  a  1  input  array  of  t  :>•  ■  di  tvi:t  ion  aaia-':;  of  tin 

r.  v ;  i  .  -  .  ,  K  ( k  ,  k  ,  k  )  . 

V  \* 

Ivai  li.jui  variable  equal  to  t.iu  known  r.  eoordirinti 
of  fir  i  i.t .  rs-'ft  ion  a:;  found ,  for  -  xamt  In,  from 
M al  -n  mt  i  i...  OGI  VK  . 

II  -  Ki  a  i  out  j  array  uontuininq  tin  •  i.  si  ivd  point. 

• ;  i  tit  or  soot  ion  PI  ( x  ,  y ,  . 

i mime: . :  :  ■■  am  i  Mt  t  in- .  i 

'in-  i.iraiT.t  trio  -juat  ions  for  u  I  irt  in  spnoo  pass  l  rut  ti-.rouqli  the 

x  ,  v  ,  r  )  and  liavinq  diroot  ion  cosines  (k  ,  k  ,  k  )  atm  qivon  bv 

O  o  (.  XV  Z 

(la) 

t  tb) 


+  k  t 


(ki 


Chapter  19 


SUBROUTINES  BDISK,  BDISKN ,  TDISK,  TDISKN 

I"-J.  Purpose:  To  compute  the-  intersection  PI(x,  y,  z)  of  a  1  iru-  (ray) 

einanut-  ing  From  the  point  P(x  ,  y  ,  z  )  Ibviih;  oireot  ion  cosines 

o  -  a 

K(k  ,  k  ,  k  )  with  a  planar  disk  at  z  =  z  or  at.  z  -  z 
x  y  z  hot  tot 

Subroutine  BDISKN  is  used  to  compute  tiie  unit  inward  normal  n  =  z. 

Subroutine  TDISKN  is  used  to  compute  the  normal  n  =  -z . 

19-2.  Usage: 

CAUL  BDISK  (I-,  K,  PI,  filT)  CALL  TDISK  (P,  K,  PI,  HIT) 

C  OMMON /BDISKC/ZB  OT ,  RBSQ  COMMON/TDISKC/ZTOP,  RTSQ 

CALI,  BDISKN  (N)  CALL  TDISKN  (N) 

1 1  ■  —  3 .  Arguments 

P  -  Real  input  array  containing  the  point  P(x  , 

o 

y  .  •  )  from  which  tno  ray  emanates. 

’  ) 

K  -  Fa  u i  incut  array  of  direction  cosines  K(k  ,  k  ,  k  ). 

x  y  z 

PI  -  Real  input  array  containing  the  desired  point 

of  intersection  PI(x,  y,  z)  .  , 

iil'i  -  Logical  output  variable  which  is  TRUE  if  an  1 

i  ntersoction  is  found. 

lit ’T  -  Real  input  variable  equal  to  the  z  c  ■ordinate  of 

the  : lunar  di sk . 

RBSp  -  Real  inj  ut  variable  equal  to  tin  square  of  tin 


nii ]  i  us  of  tin  |  ! , mar  d j  : :k  . 

Real  (  up  ;  ,  i  .inuy  con  t  a  l  ii  i  ng  I  is  d  i  i  •  c  M  on  cos  i  lies 
t  t  i  i<  unit  mold  nc  n  mu  I  Vrcln  ;  vis.,  N  ( ; ' ,  ,  1  )  . 


’<> 


DISK  *C9IZ0NTAt 


DISK  HORIZONTAL 


1 
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Chapter  20 
SUBROUTINE  FAR 

20-1.  Purpose:  To  compute  the  far  field  pattern  in  wavenumber  coordinates 

(k  ,  k  )  of  an  antenna  who.se  radiating  characteristics  are  specified 
x  y 

by  the  com  I  lex  j  lane  wav-,  spectra  A  (k  ,  k  )  ,  A  (k  ,  k  )  .  Thu  an- 

x  x  y  y  x 

tenna  is.  Locat'd  n,  a  i  1-an  perpendicular  to  the  z  (polar)  axis. 

20-2.  Usage:  CALL  PAR  (FU.U),  X FI ELD,  YFIELD,  NX,  NY,  FGHZ,  KXMAX,  KYMAX, 
RADIUS,  i  i’Wk,  *1- v X ) 

20-3.  Arguments 

FIELD  -  Real  >  at  put  irrav  of  MX  by  NY  elements  containing 

tile  far  field  power  pattern  at  discrete  wavenumbers 
k  -  sinOc-'si  ,  k  _  =  sinOsirnt,  where  0  and  ip  are  tile 
usual  polar  and  azimuthal  angles. 

XJ'IELD,  -  Con;  lex  input  arrays  of  NX  by  NY  •>leraenfs  con- 

YFIELD  taming  t  in  t  lane  wave  spectra  A  ,  A  at  discrete 

x 

wavenumbers  k  ,  k  . 

x  y 

NX,  NY  -  Integer  input  variables  equal  to  the  array  sixes. 

Kiln 2  -  Real  input  variable  equal  to  the  frequency  Ln 

gigahertz  . 

KXMAX, KYMAX-  Real  input  variables  equal  to  the  maximum  wave- 
number  associated  with  the  elements  of  the  arrays 
FIELD,  XFIEi.D,  and  YFIELD.  The  element  I  1,  J=1 
in  these  arrays  corresponds  to  the  wavenumber 
eoordinatu  (-KXMAX,  -KYMAX) .  For  any  (I,d),  the 
wavenumber  coordinates  are  given  1  >y 


1  u 


NX 

KX  =  (I  -  — -  -  1)*  KXINC 

NY 

KY  ^  (J  -  - )  )  *  KYINC 

where 

KXINC  2 * KXMA X/N X 

KYINC  =  2*KYMAX/NY 

RADIUS  -  Real  input  variable  equal  to  the  radius  r  in 

c..  n  time  tors  of  the  sphere  on  which  the  far  field 
pattern  is  computed.  This  variable  effects  only 
tlie  term  e  jKr/r,  and  r  is  set  to  unity  in  the 
calling  program  for  normal  use. 

IPWR  -  Integer  input  variable  which  selects  the  vector 

components  to  be  used  in  coin]  ut:  i.ng  t  iie  power 
i  .1 t:  !  ern : 

1  -  K leva i  ion  component  only 

1  Azimuth  component  only 

I  -  Total  power 

•5  Right  hand  circular  polarisation 
be ft  hand  c  reuiar  polarization 
l 'MAX  -  i  i  i  input  and  output  variable.  t>n  input,  if 

! 'MAX  -  0,  the  program  will  normalize  the  array 
bUlLD  fro  zero  to  one  and  output  the  normalisin' 
factor  as  KMAX .  Tf  FMAX  >  O  ■  mi  intul,  it  will  b 
•i:-,..-d  as  tiii  norma)  is.  inn  factor;  on  out :  ut  it  wil 


in e 'hanqi  1 1  . 


20-4  . 


Comments  and  Method 


Let  K  (x,  y,  0)  ,  K  (x,  v,  0)  be 
x  y 

a  rectangular  antenna  aperture  located 
tlie  oriqin  of  the  coordinate  system. 
Cure  fields  are  defined  by 


the  tangential  electric 
in  the  z  plane  and  cen 
The  {.'lane  wave  spectra  of 


A  ( k  ,  k  ) 
x  x 


1  (•  f  +i  (k  x  +  k ^  y ) 

- —  ;  i  E  (x,  v,  0)  o  X  ■  dxdv 

(271)  2  J  i  X  ' 


A  (k  ,  k  ) 

v  *  v 


( 2t ) ' 


E  ( x , 

S' 


y> 


o) 


+j(kx  +  k  y) 
e  '  dxdv 


A  ( k  ,  k  ) 

Z  X  V 


-k  A  -  k  A 

2L2L _ JiLJb 


t  ■  • 


k  f-  k  +  k 
x  y 


-  k“--  ( 


Th -■  electric  field  ..  (  «'  x,  y,  or  .;)  ..t  any  point  (>:, 

given  by 


i  A  ( k  f  k  )  e  dk  ok 

j  ■  -  v  v 


1  xk  * 
x 


Holds  of 
tv-red  at 
the  aper- 

(1) 

(2) 

(3) 


(4) 

:•  01 


) 


tui  for  tin'  s:  ivia  1  ima'  of  large  r ,  the  rectangular  field  component 
a;  :  roach  t  heir  as\!:a  t.otie  values  [1] 


ir,  k  ,  k  )  ~  ll'iik - cost'  A  (k  ,  k  ) 

xo  vn  •  r  ?  xo  vo 


.  i !  L.  !:at".  :  Jet  ■ 


I oi nt s  ar>  given  by 


k  -  k  sin-  cosy 

xo 


k  -  k  sinf-  si:  ;■ 
VO 


k  k  fost 


•it-  a!  ov  i  .mat  ions,  is  Me-  ;ela:  ancle  measured  from  the  r.  axi 
;  .is  t  tv  •  tr  iiautha  !  ..nai,-  nv-.isur-b  from  +x  toward  +y  in  tin  convi 
.  i  .  e !  n>  r  i  •  \i.l  co  iina!  •  riairs  :  . 

Consiiit  r  i  anio.  iii.a  a  s::  rone  1. 1  oo  ■rain-ite  system  '.a  !■’’  muV'  1 


v.v>V'  nnmi  •<  r  ..  ,  k  ,  ■■  ;:i. 


is  ■ :  <  r-  e  .  Jo-:  ct  s  j  u;  !  :»  oiri  et 


■  :u  i  va  k  iitiv  :  .v  ( ■  ,  .  i 


sin  '.  v  v  .  •  >r:;  •  ,  <  may  ! wv l  :  ■  eii  as 


nor;:. a  !  ,i  ::i.m  v.'.i 


.  ;x.r 


•  r 

v'. 


t  !i-  (a  fob 


(  -  elevation  component 

A 

a  =  azimuth  component 


,»iui  is  ■  receiving  ;  .a  !.•  -rn  win.’ it  tin  probe  is  ;  <.'lariz.it ion  mat.  eh*,  .-a  at 
every  ;  i  mt.  it..  Liu-  t  antenna . 

lu  tin-  of  circularly  polarized  fields,  t!.e  | -robe  n  can  be 

ex;  tossed  as 


KiiC:  n 


C 


+ 


LUC :  n 


( 20 ) 


(21) 


The  receiving  ; -attorns  in  tin-  two  cases  are  then  given  by 


where  tile  ai  crot  r iati-  n.  is  used. 

n 

Subroutine-  FAR  implements  t  in.-  above  equations  and  computes  the  power 


at  terns  for  an 

aperture  in  an 

infinite  ground 

plane;  i.e.,  t 

'.ho 

USe  of 

mlv  A  and  A 

tantamount-  tc 

-i  t  hi  assuim  t  ion 

that  E  out t 

;.i  dt. 

■  the  finit. 

X 

— t  a  n 

a  r  t  uire  a  r-  -a  j  : 

-  /.tiro.  For  ti.t 

■  -  X  t  -  - 1  lied  case  o 

if  a  finite  apr 

■r  l  .. 

.re  .it!  fri  i 

t 

t  ainu  tit  ia  1 

-  -t  i 

o  f  1  e  1  d  H  1 1  so 

-•!  an 

coni  r  ibut.es 

to  t.ii-  radia 

•  U'i>  i 

t -  far  fi- 

:ids  ,m- 

aiv-  by  liquation 

,s  (l-4b)  - 

(3-49)  of 

■> 

Tr.  fact 

,  it  i 

.  ail  y  t  y  i.m  ■ !  ud  i  r.q 

the  effects 

of  ii  that 

-tan 

H  i  .‘To  i 

t  t  tn-i  and 

■o*  •  i  v  i  r  . 

•a  form;  i  at  i  ons  foi 

the  finite 

aperturi  can 

to  : »  •  qui  va  1  out  [-1  ]  . 

Tin  our  r-  nt.  veroi  - -r,  of  Pul -rout  i.’it-  FAR  listed  below  could  bo  easily 
m  :•  i  i  ?'  l  •  -d  f  :  i  tie  1  ud«'  tie  add  i  t  i(ii,il  terms.  Tf  the  aoomo  tr  ica  1  optics 
r.  :  >i  Aim.it  for  *  in-  .it-’ftuie  f  i  ■  Ids  i:  mad",-  vis.. 
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Chapter  21 


SUBROUTINE  AMPHS 


21-1.  Purpose:  To  convert  a  complex  number  c  =  x  +  jy  from 

i  j  j  $ 

rectangular  to  polar  form  c  — |c|e 
21-2.  Usage:  CALL  AMPHS  (C,  AMP,  PHS) 

21-3.  Arguments 

C  -  Complex  input  variable  containing  the  rectangular 

components  of  the  complex  number  to  be  converted; 


AMP 

PHS 


i .e. ,  C  =  CMPLX(X,Y) . 

J~2  2 

Real  output  variable  equal  to  /x  +y  . 

Real  output  variable  equal  to  the  phase  angle  <j> 
in  degrees. 


21-4.  Comment 

The  intrinsic  Fortran  function  A TAN 2  is  used  to  compute  PHS. 
21-5.  Program  Flow:  See  listing  below. 

21-6.  Test  Case:  None 
21-7.  References:  None 

2i-8.  Program  Listing:  See  following  page. 
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Chapter  22 


SUBROUTINE  DBPV 

22-1.  Purpose:  To  convert  a  real  array  of  linear  values,  normalized  to 
lie  between  zero  and  unity,  to  decibels. 

22-2.  Usage:  CALL  DBPV  (FIELD,  NX,  NY,  IPV) 

22-3.  Arguments 

FIELD  -  Real  input/output  array  of  NX  by  NY  elements:  on 

input,  it  contains  the  values  to  be  converted;  on 

output,  it  contains  the  corresponding  decibel  values 

on  the  range  (-40,  0) .  All  input  values  less  than 
-2 

10  are  set  to  -40  dB  on  output. 

NX,  NY  -  Integer  input  variables  which  specify  the  size  of 
the  array  FIELD. 

IPV  -  Integer  input  variable  which  specifies  whether  the 

input  values  in  FIELD  represent  power  (IPV=1)  or 
voltage  (IPV=2).  If  IPV=1,  F(I,  J)  =  10  log1Q 
F(I,  J)  is  returned;  if  IPV=2,  F(I,  J)  =  20  log^ 
F(I,  J)  is  returned. 

22-4.  Comments 

It  is  intended  that  the  xrout  array  FIELD  be  normalized  prior  to  the 
call  to  Subroutine  DBPV. 

22-5.  Program  Flow:  See  listing  below. 

22-6.  Test  Case:  None 
22-7.  References:  None 

22-8.  Program  Listing:  See  following  page. 
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Chapter  23 


SUBROUTINE  NORMH 

23-1.  Purpose:  To  normalize  a  two-dimensional  real  array  of  field  values 
so  that  all  values  in  the  array  lie  between  zero  and  unity. 

23-2.  Usage:  CALL  NORMH  (FIELD,  IMAX,  JMAX,  LDB) 

23-3.  Arguments 

FIELD  -  Real  array  of  IMAX  by  JMAX  elements.  On  input, 

it  contains  the  field  values  expressed  as  non¬ 
negative  real  linear  amplitude  or  as  amplitude 
in  decibels.  On  output,  the  linear  amplitudes 
are  replaced  by  their  scaled  values  FIELD(I,J)/ 
FMAX,  where  FMAX  is  the  maximum  amplitude  value 
in  the  array;  the  logarithmic  amplitude  values 
are  replaced  by  (FIELD(I,J)+40.)/40. ,  where  -40 
decibels  is  assumed  to  be  the  lower  bound  on  the 
original  data. 

IMAX, JMAX  -  The  number  of  elements  in  FIELD. 

LDB  -  A  logical  variable  set  TRUE  if  the  values  in 

FIELD  are  in  decibels. 

23-4.  Comments  and  Method 

A  function  f(x,y)  of  two  variables  having  minimum  value  f  .  and 

min 

maximum  value  f  may  be  normalized  to  O^f  (x,y)<l  according  to 
max  n 


provided  that  the  denominator  is  not  zero.  In  this  procedure,  the  f  =0 

corresponds  to  f=f  .  ,  and  f  =1  corresponds  to  f=f 

min  n  max 

When  f(x,y)  represents  a  linear  (vice  logarithmic)  variable,  it 

is  desirable  to  force  f  .  to  be  zero  if  the  minimum  value  of  f  is  actually 

min 

greater  than  zero.  In  this  special  case,  f  becomes 

n 


f  (x,y)  -  pal- 
n  f 

max 


Equation  (2)  is  also  used  to  treat  the  special  case  of  f  -  f  .  -  0; 

max  min 

however  if  If  I <1 ,  f  is  set  equal  to  ±1.,  where  the  sign  used  is 
'  max '  max 

that  of  f  .  This  refinement  has  the  effect  of  producing  a  constant 
max 

function  whose  value  lies  between  zero  and  unity;  without  it,  f^  would 
be  simply  set  to  unity  or  division  by  zero  may  result. 

When  f (x,y)  represents  a  logarithmic  variable,  such  as  the  ampli¬ 
tude  in  decibels  of  an  electromagnetic  field,  all  of  the  foregoing  dis¬ 
cussion  applies;  however,  a  minimum  value  f  .  must  be  imposed.  If 

min 

f  .  <-40,  f  .  is  set  equal  to  -80  (decibels) ;  otherwise,  a  -40  decibel 
min  min 

level  is  assumed.  A  value  of  f  equal  to  zero  decibel  is  also  assumed. 

max 

23-5.  Program  Flow 

Lines  9-16:  Find  minimum  MN  and  maximum  MX  values  of  data  in 

FIELD;  form  their  difference  DR=MX-MN. 

Line  17:  If  array  values  are  in  decibels,  go  to  50. 

Line  18:  If  all  values  in  the  array  are  the  same,  go  to  25 

and  scale  the  data  to  lie  between  zero  and  unity 
(Lines  28-37) . 
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Lines  19-27:  If  all  linear  amplitude  values  in  FIELD  are  not 

identical,  scale  the  data  according  to  FIELD (I, J)  = 
(FIELD(I , J)  -  Min.  Value) /(Maximum  Value  -  Minimum 
Value) . 

Line  38:  If  values  in  FIELD  are  in  decibels,  and  the  minimum 

value  is  less  than  -41dB,  then  assume  a  -80dB  lower 
bound,  go  to  60  (Lines  47-52) ,  and  scale  the  data 
according  to  (FIELD(I,J)  +  80.)/80. 

Lines  39-46:  Scale  the  data  according  to  a  -40dB  lower  bound;  i.e. 
(FIELD (I , J)  +  40.)/40. 

Lines  53-54:  Write  MN  and  MX. 

23-6.  Test  Case:  See  Chapter  2. 

23-7.  References:  None. 

23-8.  Program  Listing:  See  following  pages. 
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Chapter  24 


SUBROUTINE  CNPLTH  AND  FUNCTION  PS I 


4-1.  Purpose:  To  plot  (Calcomp)  single  dimensional  far  field  patterns 

at  constant  wavenumber  k,.  . 

fix 

4-2.  Usage:  CALL  CNPLTH  (FIELD,  N,  KMAX,  KCNTR,  KFIX) 


PSI  =  AT  AN  2  (K//l.  -  K2  -  K^.  ) 

fix 


4-3.  Arguments 


FIELD 


Real  input  array  of  N  elements  containing  the 


field  values  in  decibels  but  normalized  so  that 
-40  dB  corresponds  to  0  and  0  dB  corresponds  to 
unity  on  the  normalized  scale. 

Integer  input  variable  which  specifies  the  number 
of  elements  in  FIELD. 

Real  input  variable  equal  to  the  half  width  of 
the  wavenumber  range  corresponding  to  the  array 
elements  1  through  N  of  the  array  FIELD;  i.e.,  the 
increment  in  wavenumber  corresponding  to  the  dis¬ 
tance  between  the  Ith  and  (I+l)st  element  is 


2  KMAX/N . 


KCNTR 


Real  input  variable  equal  to  the  wavenumber  coor¬ 


dinate  of  the  (N/2  +  l)st  element  of  the  array 
FIELD.  FIELD(l)  has  wavenumber  coordinate 


KCNTR  -  KMAX. 


Real  input  variable  equal  to  the  fixed  value  of 
the  other  wavenumber  coordinate.  For  example,  if 


k  varies,  then  k  =  KFIX. 
x  y 


.  UE  biuU*K-N0i  blU'iL 


24-4.  Comment  and  Method 

Let  F(k^,  k  )  represent  the  far-field  power  pattern  of  an  antenna 

where  k  and  k  are  normalized  wavenumbers  as  defined  in  Chapter  20.  A 
x  y 

pattern  cut  at  constant  wavenumber  is  a  conical  cut  about  the  real  axis; 

e.q.,  a  k  =  constant  cut  is  a  conical  cut  about  the  x  axis  of  the  coor- 
x 

dinate  system. 

For  principal  plane  cuts,  =  0  yields  an  E-plane  pattern  as 
defined  in  Figure  2-3;  k^  =  0  yields  an  H-plane  pattern.  For  principal 
plane  cuts,  KCNTR  =  0  and  KFIX  =  0. 

The  plotting  commands  are  set  up  to  produce  a  4"  X  8"  rectangular 
pattern  plot  on  a  standard  pattern  scale.  The  plot  is  positioned  on  the 
paper  to  give  margins  of  2"  on  the  left,  1"  on  the  right,  and  2.25"  from 
the  bottom  and,  hence,  is  suitable  for  direct  use  as  a  figure  in  a  tech¬ 
nical  report. 

24-5.  Program  Flow:  See  listing  below. 

24-6.  Test  Case:  See  Chapter  2  and  pattern  plots  in  Appendices  B  and  D. 
24-7.  References:  None 

24-8.  Program  Listing:  See  following  pages. 
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Chapter  25 


SUBROUTINES  PLT3DH  AND  PLTT 


25-1.  Purpose:  To  plot  (Calcomp)  the  two-dimensional  array  FIELD  (I,  J) . 
25-2.  Usage:  CALL  PLT3DH  (XSIZE,  YSIZE,  HEIGHT,  FIELD,  IMAX,  JMAX,  NMZ, 
LDB) 


25-3.  Arguments 
XSIZE, 
YSIZE, 
HEIGHT 
FIELD, 
IMAX, 
JMAX 

NMZ 


LDB 


Real  input  variables  in  inches  defined  on  Figure  1. 


Real  input  array  of  IMAX  by  JMAX  elements  con¬ 
taining  the  values  to  be  plotted.  These  values 
must  be  normalized  to  the  range  (0,  1)  before 
plotting. 

Logical  input  variable.  If  NMZ  =  .TRUE.,  the 
array  FIELD  will  be  normalized  with  respect  to 
its  own  maximum  value;  if  NMZ  =  .FALSE.,  no 
normalization  will  be  done. 

Logical  input  variable  required  by  Subroutine 
NORMH  (Chapter  23)  . 


25-4.  Comments 

In  Figure  25-1,  the  axes  and  labels  shown  are  not  produced  by  the 
subroutine;  these  axes  are  presented  to  demonstrate  the  perspective  of  the 
plot  and  to  identify  its  dimensions.  Report  size  plots  will  be  produced 
suitable  for  one  8  1/2"  X  11"  page  when  FACTOR  =  1.0  and 
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Figure  25—1.  Dimensions  of  Three-Dimensional 
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r  X. 


XSIZE  =  6.0" 

YSIZE  =  2.6" 

HEIGHT  =  2.5" 

Margins  in  this  case  will  be  1.5"  on  the  left,  1"  on  the  right,  and  4.25" 
from  the  bottom  of  the  plot  paper.  Margin  lines  are  provided  on  the  plot 
paper  to  outline  the  8  1/2"  X  11"  page.  Also,  the  plot  itself  can  be 
carefully  cut  from  the  plot  paper  and  cemented  onto  a  set  of  axes  as  has 
evidently  been  done  in  Appendices  B  and  D. 

25-5.  Program  Flow:  See  listing  below. 

25-6.  Test  Case:  See  Chapter  2  and  Appendices  B  and  D. 

25-7.  References:  None. 

25-8.  Program  Listing:  See  following  pages. 
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Chapter  26 


SUBROUTINE  FFTA 

26-1.  Purpose:  To  compute  the  Discrete  Fourier  Transform  (DFT)  or  its 

N 

inverse  of  a  sequence  of  complex  numbers  consisting  of  2  elements, 
where  N  is  an  integer.  The  Cooley-Tukey  algorithm  is  used  to  per¬ 
form  computations  in  place  to  speed  up  the  computations  and  to 
return  the  transformed  values  in  the  input  array. 

26-2.  Usage:  CALL  FFTA  (FIELD,  NEXP,  IBMISN) 

26-3.  Arguments 

FIELD  -  Complex  array  of  2**  NEXP  elements:  on  input  it 

contains  the  sample  data  to  be  transformed;  on 
output  it  contains  the  transformed  data.  See 
below  for  ordering  of  data. 

NEXP  -  Integer  exponent)  e.g.,  for  64  elements  in  FIELD, 

NEXP=6 . 

IBMISN  -  Integer  parameter  which  controls  operation: 

IBMISN  =  3  performs  the  inverse  DFT 
IBMISN  f  3  performs  the  DFT  as  defined  in 
4  below. 

26-4.  Comments  and  Method 

a.  Subroutine  FFTA  is  machine-dependent  in  that  the  bit  reversed 
number,  IFLIP,  must  be  generated  using  Fortran  instructions  which  are 
peculiar  to  a  particular  machine.  Also,  the  word  length  must  be  taken 
into  account.  Lines  38-42  of  the  attached  program  listing  are  used  to 
effect  che  desired  operation  for  the  CDC  Cyber  70  (60-bit  word,  numbered 
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0  through  59  from  right  to  left  with  Bit  0  being  the  least  significant) : 

IFLIP=0 

DO  4  11=1,  IEXP ,  1 

J=60  -  II 

IFLIP=2*IFLIP  +  AND  (SHI  FT ( 1 , 1+J) ,  IB) 

4  CONTINUE 

The  SHIFT(I,1+J)  operation  shifts  the  bits  of  the  integer  I  to 
the  left  by  1+J  bit  positions.  The  AND  operation  strips  off  the  right 
most  bit  of  the  shifted  result.  E.G. ,  when  11=1,  the  right  most  bit  of 
I  (Bit  0)  is  extracted  from  I  by  the  AND(SHIFT)  operation.  The  current 
value  of  IFLIP  is  then  shifted  one  bit  to  the  left  by  the  2*IFLIP  opera¬ 
tion.  The  two  results  are  then  added  together.  A  total  of  NEXP  bits 
are  extracted,  starting  with  Bit  0,  followed  by  Bits  1,  2, . . . (NEXP-1) . 

The  net  result  of  these  operations  is  to  take  the  NEXP-bit  binary 
representation  of  the  array  element  number  I,  reverse  the  order  of  the 
bits,  and  right  justify  the  result.  Array  elements  in  FIELD  numbered  I 
and  IFLIP  are  then  interchanged  if  I>IFLIP.  The  first  and  last  elements 
of  FIELD  always  remain  in  place.  The  array  elements  are  rearranged  in 
this  manner  so  that  they  will  be  ordered  after  transforming  [1] . 

b.  To  explain  the  ordering  of  the  data  in  the  complex  array  FIELD, 
it  is  convenient  to  consider  the  specific  example  of  using  FFTA  to  com¬ 
pute  the  Fourier  transform  G(f)  of  a  time  function  g(t)  as  defined  by 


G(f) 


T 

max 

!  g(t) 


- j2nft 
e 


dt 


-T 

max 


(1) 
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and  as  approximated  by 


G(f)  =  E  g(ti)  e_j2  fti  At 


(2) 


where  t  are  the  equally  spaced  points  along  the  t  axis  when  g  is  sampled 

over  the  interval  -T  ^  t  <  T 

max  max 

NEXP 

There  are  N=2  samples  in  the  input  array  FIELD (X)  corresponding 

to  I=1,N.  The  first  sample  (1=1)  corresponds  to  g(-T  ).  The  last 

max 

sample  (I=N)  corresponds  to  g(T  -At).  The  I=(N/2+l)th  sample  corresponds 

max 

to  g(0) ;  i.e.,  the  value  of  g  at  t=0.  The  DFT  assumes  periodicity  of  the 

sampled  data  so  that  the  value  at  t=T  is  identical  to  that  at  t=-T 

max  max 


The  sample  spacing  is 


At  =  2 


T  /N 
max 


(3) 


and  corresponds  to  a  folding  frequency  f  of 

max 


f  =  l/2At 
max 


(4) 


On  output,  the  array  FIELD  contains  the  frequency  components  G(f) 

at  N  equally  spaced  frequencies  Af  over  the  band  -f  <  f  ^  f  ,  where 

max  max 

1=1  corresponds  to  f=-f  ,  I=(N/2+l)  to  f=0,  and  I=N  to  f=f  -Af,  where 

max  max 


Af  =  2  f  /N  (5) 

max 
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and  where 


T 

max 


1 

2Af 


Also,  by  the  inversion  integral  [2], 


(6) 


g(t) 


f 

max 

G ( f )  e+j27rft  df  *  Af  Z  G(f  )  ej2lTfpt 

f  P  P 

max 


(7) 


This  version  of  Subroutine  FFTA  is  written  so  that  division  by  N 

“  j  27Tf  t 

is  done  when  the  Fourier  transform  (kernel  =  e  )  is  computed.  When 

the  expression  in  Equation  (3)  for  At  is  used  in  (2) ,  there  results 


G(f  )  =  2  T  J  I  g(t.)  e"j27Tfpti 
p  max  N  3  i 


(8) 


Transposing  2  T  and  using  Equation  (6)  yields 
mcix 


Af  G ( f  )  =  jjj  £  g(tj  e  ^"7Tfpti 
P  i 


(9) 


where  the  righthand  side  is  the  definition  of  the  Discrete  Fourier  Trans¬ 
form  as  computed  by  FFTA.  Inversely, 


g(t.)  =  E  Af  G(f  ) 
1  P 


e 


+j2irf  t. 
pi 


(10) 


224 


] 


which  is  the  Inverse  DFT  as  computed  by  FFTA. 

Conversely,  if  the  original  data  in  the  input  array  FIELD  are 
samples  of  a  frequency  spectrum  G(f) ,  a  similar  analysis  shows  that  FFTA 
computes  At  g(tj  as  the  inverse  transform  (IBMISN=3) ;  i.e.,  the  time 
function  is  modified  in  amplitude  by  At.  Of  course,  when  the  forward 
transform  (IBMISN7O)  is  performed  on  this  result,  tne  original  sampled 
data  G(f J  are  obtained  in  FIELD  on  output. 

From  the  above  considerations,  the  following  conclusions  can  be 
drawn  concerning  the  use  of  FFTA  to  compute  the  Fourier  transform  G(f) 
of  a  windowed  time  function  g(t) : 

G(f  )  =  2  T  •  FFTA{ g ( t . ) }  (11) 

p  max  x 

g(t.)  =  — ^ -  •  IFFTA{G(f  )}  =  IFFTA{ Af  G(f  ))  (12) 

1  2  T  p  p 

max 


As  an  example,  let  g(t)  be  the  rectangular  pulse  function  which 

has  constant  amplitude  V  for  I t I <t  and  which  is  windowed  in  the  larger 

o  '  1  o 

time  interval  |t|<T  .  The  Fourier  transform  G(f)  is  given  by  [3] 

max 


G  ( f ) 


sin  2nft 

2  t  V  — — — - 2 

o  o  2nft 

o 


(13) 


Let  g(t)  be  sampled  at  N=2NE*P  points  over  the  interval  I t | sT  ,  and  let 

11  max 

these  sampled  poir,  s  be  placed  in  the  array  FIELD.  Then  the  si>ectrum 
G ( f )  will  be  closely  approximated  at  discrete  frequencies  f  by 
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*  FIELD ( I ) 


G(f  )  ~  2  T 
p  max 


where 


f  =  -f  +  (1-1)  *  Af 
p  max 

and  where  FIELD  is  the  output  of  FFTA  according  to  CALL  FFTA  (FIELD,  N,  0) . 

Proper  consideration  should  be  given  to  the  sampling  of  the  time 

function  so  that  the  DFT  produces  a  good  estimate  of  the  actual  integral 

transform.  For  example,  if  t  =T  ,  and  all  samples  are  constant,  then 

o  max 

the  DFT  will  produce  a  single  nonzero  frequency  component  at  f=0  (corres¬ 
ponding  to  the  (N/2+l)th  element  of  FIELD);  i.e.,  a  delta  function.  Such 
a  result  follows  from  the  facts  that  the  Fourier  transform  of  a  constant 

g(t)=V  is  G(f)=V  6(f)  and  that  the  DFT  assumes  a  periodicity  of  the 
o  o 

sequence  of  samples  provided  to  it. 

Consider  the  other  extreme.  Let  the  pulse  g(t)  be  represented  by 

only  one  sample  at  t=0  in  the  window  I t I <T  .  The  Fourier  transform  of 

1  max 

g(t)=V  5(t)  is  G(f)=V  ,  a  constant, 
o  o 

It  is  clear  from  the  above  considerations  that  the  time  function 
must  be  properly  windowed  and  properly  sampled  to  produce  a  good  estimate 
of  its  transform  via  the  DFT.  Simply  stated,  the  time  function  should  be 
sampled  at  a  rate  At  which  is  twice  the  highest  frequency  contained  in 
the  function  as  interpreted  by  the  DFT. 


*• 


~  *i; 
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26-5.  Program  Flow 


NEXP 

Lines  22-24:  Compute  N=2  and  set  the  sign  ISN  of  the  exponent 
in  the  Fourier  kernel. 

Lines  26-29:  Compute  IEXP=NEXP  from  N.  This  is  a  redundant  com¬ 
putation  made  when  the  original  FFT  subroutine  was 
modified  to  conform  to  the  call  to  a  library  version 
on  another  computer  system. 

Lines  30-35:  Rearrange  the  order  of  the  input  data  so  that  samples 
for  t>0  are  placed  in  the  lower  half  of  the  array, 
and  those  for  t<0  are  placed  in  the  upper  half.  For 
a  frequency  function,  the  data  are  rearranged  so 
that  the  first  N/2  points  give  the  components  for 
non-negative  frequencies  (1=1  corresponds  to  f=0) , 
and  the  last  N/2  points  contain  the  data  for  the 
negative  frequencies. 

Lines  36-49:  Rearrange  the  data  in  FIELD  so  that  it  will  be  ordered 
after  transforming  as  described  for  Lines  30-35  above. 

Lines  50-73:  Perform  the  summation  using  the  Cooley-Tukey  algorithm 

[1]. 

Lines  74-79:  If  forward  transform  is  being  done,  divide  all  values 
in  FIELD  by  N. 

Lines  80-85:  Rearrange  the  output  data  in  FIELD  so  that  it  conforms 

to  that  used  on  input;  i.e.,  f.=f  +  (I-l)Af  or 

l  max 

t . =-T  +  (I-l)At  as  appropriate. 

i  max 

26-6.  Test  Case 

A  rectangular  pulse  function  with  amplitude  Vq=100  was  chosen  for 

g(t)  with  t  =.10  second  T  =1.60  seconds,  and  N=2048=21'*' .  The  resulting 
o  max 
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sample  increment  At  and  folding  frequency  were  0.116  second  320.0 
Hertz,  respectively.  The  comparison  of  the  central  nine  points  of  the 
computed  and  true  frequency  spectra  were  as  follows  (CDC  Cyber  70) : 


True  6(f)  Computed  G(f) 


I 

f  (Hz) 

Am}'. 

Phase ( °) 

Amp . 

Phase ( °) 

1021 

-1.250 

18.006 

0.0" 

18.006 

0.35 

1022 

-0.938 

16.863 

0 . 00 

18.863 

0.26 

1023 

-0.623 

1". 490 

0.00 

19.490 

0.18 

1024 

-0. 313 

19.872 

0.00 

19.872 

0.09 

102  5 

0.000 

20.000 

0.00 

20.000 

0.00 

1026 

0.  313 

19.872 

0 . 00 

19.872 

-0 . 09 

1027 

0.625 

19.490 

0.00 

19.490 

-0.  26 

1028 

0.938 

18.863 

0.00 

18.863 

-.033 
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26-8.  Program  Listing:  See  following  pages.  The  second  listing.  Subroutine 
FFT,  is  for  use  on  the  IBM  3033  at  JHU/APL.  It  employs  the  subroutine  FFTA 
available  on  that  system  library.  Use  of  this  subroutine  requires  the  calls 
in  Subroutines  JOYFFT  and  MAGFFT  to  be  changed  from  CALL  FFTA  to  CALL  FFT. 
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RETURN 


Chapter  27 


SUBROUTINE  MAG F FT 

27-1.  Purpose:  To  increase  the  resolution  of  a  complex  array  of  data 

points  using  Fourier  interpolation  and  the  Fast  Fourier  Transform. 
The  number  of  points  in  each  array  must  be  an  integer  power  of  two. 
~J- 2.  Usage:  CALL  MAGFFT  (A,  NA,  B,  NB) 

27-3.  Arguments 

M 

A, NA  -  Complex  input  array  of  NA  =  2  <  NB  data  points. 

N 

B ,  NB  -  Complex  output  array  of  NB  =  2  data  joints. 

27-4.  Comment  and  Method 

a.  Subroutines  required:  FFTA ,  PWRTWO 

b.  By  Shannon's  sampling  theorem,  a  band-limited  function  is 
represented  by  its  samples,  and  it  can  be  reconstructed  at  any  joint  from 
them.  The  computation  of  the  value  of  the  function  at  a  joint  other  than 
a  sample  point  is  called  Fourier  interpolation.  Such  inter jolut ion  can  be 
used  to  increase  the  resolution  of  a  function. 

The  Fast  Fourier  Transform  ( F FT )  can  be  used  to  facilitate  Fourier 

interpolation.  Briefly,  the  original  function  A(k  ),  know*:,  at  NA  joints 

x 

on  the  range  (-K  ,  +  K  )  ,  is  transformed  to  yield  F.tx)  FiAik  1  }  at  NA 
MM  x 

sample  points.  These  NA  values  of  K(x)  are  then  j  laced  in  tn*.  center  of 

N  M 

an  array  containing  NB  -  2  *  NA  joints  to  fi  rm  the  function  1. 1  (x)  . 


TIi  i  s 

f  urict  ion 

is  then 

1  liver 

transformed  to 

produce  Al 

:k  )  at  NR 
X 

joint  s 

over 

the  same 

range  ( 

-K  ,  > 

M 

K  )  .  (Actually, 
M 

the  range 

is  ( -K  , 

M 

+  K  - 
M 

Ak; 

since  the 

FFT  con 

sider  s 

the  sum:  led  function  to  be 

periodic 

outside 

the  known  rang*.'  so  that  the-  (NB  +  l)st  joint  would  bo  the  same  as  the 
first  point  in  the  array.) 

2  3‘j 

l  rt.  ..  .  ■  »  .  ...it 


27-5.  Program  Flow:  See  listing  below. 


27-6 

27-7 

27-8 


Test  Case:  See  Chapter  2. 

References:  See  Chapter  26. 

Program  Listing:  See  following  page. 
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Chapter  28 


SUBROUTINE  JOYFFT 


28-1.  Purpose:  To  compute  the  two-dimensional  Fast  Fourier  Transform 

of  a  complex  array  of  NXI  by  NYI  points  and  to  provide  magnification 
of  a  specified  portion  of  the  transformed  data. 

28-2.  Usage:  CALL  JOYFFT  (INPUT,  NXI,  NYI,  MX,  MY,  NXC,  NYC,  OUTPUT, 

NXO,  NYO,  XYFFT,  NXY ,  ISN) 


28-3.  Arguments 
INPUT, 
NXI,  NYI 
MX,  MY 

NXC,  NYC 

OUTPUT, 
NXO,  NYO 

XYFFT, 

NXY 

ISN 

28-4 .  Comment 


Complex  input  array  of  NXI  by  NYI  points. 

Integer  input  variables,  equal  to  an  integer  power 
of  two,  which  specify  the  magnification  in  the  I 
and  J  directions,  respectively. 

Integer  input  variables  which  specify  the  center 
coordinate  I  =  NXC,  J  =  NYC  of  the  sector  to  be 
magnified . 

Complex  output  array  of  NXO  by  NYO  points  con¬ 
taining  the  transformed  points  of  the  magnified 
sector . 

Complex  working  array  of  NXY  points. 

Integer  input  variable  which  specifies  the  direction 
of  the  FFT :  ISN  =  3  for  inverse  FFT;  ISN  =  1  for 
FFT.  See  Chapter  26. 


a.  Subroutines  required:  FFTA ,  PWRTWO. 
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b.  All  integer  input  variables  must  be  integer  powers  of  2  and 
must  satisfy  the  following  restrictions: 

(1)  NXO*NYO  <  NXI*NYI 

(2)  NXO  <  NXI  or  NYO  <  NY I 

(3)  MX*NXI  <  NXY  and  MY*NYI  <  NXY 
28-5.  Program  Flow:  See  listing  below. 

28-6.  Test  Case:  See  Chapter  2. 

28-7.  References:  None 

28-8.  Program  Listing:  See  following  pages. 
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Figure  D-3.  Phase  of  E  of  Flat  Plate  Antenna. 


Figure  D-8.  Transmitting  E-Plane  A  Pattern  of  Flat  Plate  Antenna 
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Figure  D-10.  Phase  of  E, 


VAA7. 


of  Flat  Plate  Antenna 


Figure  D-14.  Receiving  H-Plane  Sum  Pattern  of  Flat  Plate  Antenna  With  Radome 


Figure  D-15.  Receiving  H-Plane  A  Pattern  of  Flat  Plate  Antenna  With  Radonte 


Appendix  E 


Plane  Wave  Transmission  Through  Multilayered  Radome  Wall 
(Excerpted  from  Reference  1  cited  in  Chapter  11.) 

The  derivation  below  and  the  computer  program  implementation 

listed  in  Appendix  F  are  based  on  work  done  by  Richmond  at  Ohio  State 
University.  Although  Richmond's  matrix  formulation  for  the  analysis 
of  plane  multilayers  has  been  previously  documented  [3],  an  outline  of 
the  theory  is  repeated  here  to  provide  a  convenient  reference  in 
defining  the  quantities  described  in  the  computer  program  of  Appendix  F. 

Consider  a  plane  electromagnetic  wave  incident  on  the  surface  of 
a  stack  of  plane,  homogeneous,  dielectric  slabs  of  finite  thickness 
and  infinite  width  surrounded  by  free  space  as  shown  in  Figure  7(a). 

The  wave  illustrated  has  perpendicular  polarization  (electric  field 
intensity  vector  perpendicular  to  the  plane  of  incidence)  and  the 
symbols  E_^  and  represent  the  electric  field  intensities  of  the 
incident  and  reflected  waves  at  the  "incident  point  "  P,  and  repre¬ 
sents  the  electric  field  intensity  of  the  transmitted  wave  at  the 
"normal  exit  point  "  Q.  The  reflection  coefficient  R  and  the  "normal 
transmission  coefficient  "  T^  of  the  multilayer  are  defined  by 

VP) 

R  =  z’.'C?)  (perPendicular  polarization)  (171) 

and 

E  (Q) 

Tn  =  E* '(?")'  {perPen<*icular  polarization)  (172) 
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The  "insertion  transmission  coefficient"  T  is  defined  as  follows 

E  (Q) 

T  =  —  (perpendicular  polarization)  (172) 

_  T  ejkdcos6 
n 

where  d  is  the  total  multilayer  thickness,  8  is  the  angle  of  incidence 

measured  from  the  normal,  and  k  is  the  free-space  phase  constant 

uj/u  e  =  2tj/A  . 
o  o  o 

The  resultant  field  in  each  layer  consists  of  an  outgoing  wave 
and  a  reflected  wave.  In  Figure  l  the  complex  constants  and 
represent  the  electric  field  intensity  E^  of  the  outgoing  wave  in  layer 
n,  evaluated  at  its  two  boundaries,  and  and  represent  the  reflected 
field  intensity  at  the  two  boundaries. 

The  field  intensity  in  layer  n  can  be  written  as 

_  ,  n  .  ,  n  n  -jkysinB 

=  (ae  +  be  )e  J  J  (173) 

The  propagation  constant  yn  is  expressed  in  terms  of  the  attenuation 

constant  a  and  the  phase  constant  6  as 
~  r.  n 

y  =  a  t  jg  (174) 

n  n  J  n 

It  is  assumed  that  the  permeability  of  each  layer  is  real  and  the  com¬ 
plex  permittivity  is  expressed  as 
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e  =  e ' (l-jtan6  > 


(17b) 


Using  the  wave  equations  and  Equations  (173),  (174),  and  (175),  it  can 
be  found  that 


=  (k//2  ) 


r\  /■ 

-  sin  0)  +  (y  e  'tan6)' 

r  r 


<Vr  "  sin  0) 


6  =  (k//f  ) 


(ye'-  sin^S)^  +  (y  e'tan6)^  +  (ye'  -  sin20) 
r  r  r  r  r  r 


where  y^  and  e^  are  the  relative  permeability  and  permittivity: 


yr  =  y/yQ  (178) 

and 

«'  =  e'/e  (179) 

r  o 


Ev.ilu.ii  i.nf.  E^  in  Equation  (173)  at  the  left  and  right  boundaries 
of  layer  n ,  it  can  be  shown  that 


A  = 
n 


C  e 
n 


-Y  d 
n  n 


(180) 


and 


B 

n 


D  e 
n 


Y  d 
n  n 


(181) 


where  is  the  thickness  of  layer  n.  Equations  (180)  and  (181)  can 
be  expressed  by  the  following  matrix  equation: 


(17G) 


(177) 
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-y  d 
n  n 


0 


0 

Y  d 
n  n 

C 

n 

e 

D 

. 

(182) 


Let  t  and  r  t ,  denote  the  interface  transmission  and 

n+l,n  n+l,n 

reflection  coefficients  for  a  wave  in  layer  n+1  incident  on  the  boundary 
of  layer  n.  Further,  let  t  and  r  represent  the  interface 

coefficients  for  a  wave  in  layer  n  on  the  boundary  of  layer  n+1.  In 
terms  of  these  coefficients,  the  electric  field  intensities,  evaluated 
at  both  sides  of  the  boundary  between  layers  n  and  n+1,  are  related 
linearly  as  follows: 


and 


C  =  t  A  ,  +  r  D 
n  n+l,n  n+1  n,n+l  n 


(183) 


Bn+1  ”  ^n,n+lDn  +  rn+l,nAn+l 


(184) 


The  relations  follow  from  the  superposition  theorem  and  the  definitions 
of  the  interface  coefficients . 

It  can  be  shown  that 


rn,n+l  “  ”rn+l,n 

t  s  l  +  r 

xn+l,n  x  n+1 ,n 

Si, n+1  s  1  +  rn,n+l=  1-rn+l,n’ 


(18b) 


(186) 


(187) 
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iliU 


t 

:  i+i 


,n  n,n+l 


r  r 

n+l,n  n,n+l 


1 


by  using  Equations  (.185)  through  (188),  Equations  (183) 
>-ur  L>i  ariunged  as 


C  -  (A  ,  +  r  ,)/t 

n  n+1  n  ,n+l  n+1  n,ntj 


Lin  d 


( b  -  r  .  A  .  )/t 
n+1  n+i,n  n+1  n,n+. 


These  oar.  ~e  <_ xoroSaeU  in  :..uti  ix  1  orin  as 


f 


1 


-r 


n+1  ,n 


n  ,n+l 


l 


-r 


n+ 


1 


,3 

i  l  nt  1 


The  matrix  equations  ( i.  e>2 )  and  (191)  can  bo  combined  to  obtain 
lowing : 


f  a 

y  d  )  i 

n  n 

n  :  i 

1 

’•■1 

.Dn  *1. 

.  ‘ 
t  . 
n  - 1 ,  r. 

n 

-y  d 
n  n 

-r  ,  e 

n,n-l 

-r  .  e 

n  ,n-l 

y  d 
n  n 
e 

1 

I 

1 

(188) 


and  (188) 


( 1 89  ) 


(  19U) 


(191) 


the  fol- 


(192) 


Let  the  two-by-two  matrix  in  Equation  (192)  be  denoted  by  Mn : 


M 


n 


-r  ,e 
n,n-l 


■Y  d 
'n  n 


Y  d 
n  n 


(19  3) 


Repeated  application  of  Equation  (193)  yields  the  following  matrix  rela¬ 
tionship  between  the  electric  field  intensities  at  the  incidence  and 
exit  surfaces : 


o 

=  (l/t)  M1  .  M2  .  M3  ...  ^  .  S  . 

ft 

an+i 

D 

o 

l  J 

bn+i 

.  J 

(194) 


where  the  dots  denote  matrix  multiplication,  N  represents  the  total 
number  of  layers,  S  denotes  the  matrix 


1 

^N+l.N 


"rN+l,N 

1 


(.  19b ) 


and 


1  '  t0,l’tl,2,t2,3* '  ^N.N+l* 


(  19b) 


In  the  situation  used  to  define  the  transmission  and  reflection 


coefficients  of  the  structure,  a  wave  of  unit  amplitude  is  assumed  to 


be  incident  on  one  outer  surface,  so  that 


097) 


B 


N+l 


R 


(198) 


C  -  T  (199) 

o  n 

and 

D  =0  (200) 

o 


Thus  liquation  (li4)  becomes 


-  ( i  / 1 )  .  k  .  m  . .  m,. 

1  dl  J 


(201) 


l 


T.'.e  o  1  ut  ion  :  o  r  "pai  alie  1  polarisation"  (electric  l.eid  citt'ic  i 
: y  ;  u  i :  ••  i  i  •  '.be  ,l.i!iv  :  incidence )  is  obti.ne*:  .-v  apply-:.,'  ' i  • 
theorem  •  duality  tv  the  above  solution.  Thus,  the  retie  t  lor.  ..nu 
*  r  i  r  s  i  1  t  i  .  *  ■  1. 1  s  •  1 1  ,  •  lift  ineci  oy 


i  i  ) 
~r'1 


\ : arul lei  polarisation ) 


i  <) 


-'■^y  ^  j  urdiitl  polar  a  ita  t  ion  t 


1*3. 


The  matrix  equations  given  above  apply  also  for  parallel  polarization, 

in  which  case  the  complex  constants  A  ,  B  ,  C  and  D  represent  the 

r  n  n  n  n 

amplitudes  of  the  magnetic  field  intensities  of  the  traveling  waves 
in  layer  n.  Equations  (185)  through  (188)  also  apply  for  parallel 
polarization,  in  which  case  the  interface  reflection  and  transmission 
coefficients  are  defined  by  the  ratio  of  the  magnetic  field  intensities 
H^.  The  interface  reflection  coefficients  are  given  by 


.  Vn+1  -  Wn 

rn+l,n  -  Vn+1  +  VlYn 


(perpendicular  polarization)  (204) 


and 


‘  n+1 ,n 


Vn+1  ~  Wn 
Vn+1  +  Wn 


(parallel  polarization) 


(205) 


where  y  is  given  by  Equations  (174),  (176),  and  (177)  if  the  permea¬ 
bility  u  of  each  layer  is  real. 

After  the  indicated  matrix  multiplications  of  Equation  (44)  are 
perlorrned,  and  the  division  by  t,  the  equation  has  the  form 


(206) 


Thus  , 


and 


T 

n 


a  +  bR 


(207) 


R  =  -c/e. 


(208) 
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